Résolution des problëmes aux limites de Laplace en cylindriques 3D par les transformations de Hankel et les 
transformations de Kontorovich-Lebedev- p 212 


Méthode de résolution des problèmes aux limites de l'équation de Laplace en coordonnées 
cylindriques 3D sur des domaines non bornés par les transformations intégrales de Hankel 


On aborde en quelques pages les problèmes aux limites cylindriques sur des domaines non bornés 
radialement et partiellement bornés sur l'axe z. Les conditions aux limites sur l'axe z sont 
inhomogènes. Les solutions de base obtenues par séparation sont les suivantes dans le cas ой les 
conditions aux limites présentent une dépendance d'angle polaire : 
T(r,0,z) = (А Cos(n0) + B Sin(n0))(C Cosh(Az) + D Sinh(Az)XE Ј (Ar) + F Y,(Ar)) 
Si les conditions aux limites sont axi-symétriques (pas de dépendance d'angle polaire), alors les 
solutions se restreignent au cas n=0 : 
T (r, z) = (C Cosh(Az) + D Sinh(Az)XE J (Ar) + F Y,(Ar)) 
Si le domaine radiale est non borné, par exemple r€[0,4 ee], dans ce cas les solutions séparables ne 
forment pas en général un système à partir duquel on peut développer une série discrète de 
fonctions propres notamment parce que les conditions de finitude des solutions ne peuvent être 
respecter en même temps en О et += pour les fonctions radiales. Si le domaine est non borné 
supérieurement, r€[r0,+2°], alors toute fonction de la forme : 
Y (Ar, )J (Ar) - Y (Ar) (Ar) 
respecte manifestement une condition aux limites homogène de Dirichlet en r=r0, mais pas la 
condition de finitude en r=+œ. 


Devant cette impossibilité de construction, on peut envisager des formes intégrales comme 
solutions d'un problème aux limites axi-symétrique non borné en r€[0,+=] et borné en z avec des 
conditions aux limites inhomogène en z€[0,z0], soit proposer une solution de la forme : 


T(r,z) = | dA J (Ar)(C(2) Cosh(Az) + D(A)Sinh(Az)) 


Les coefficients C et D sont alors devenus des fonctions du paramètre d'intégration que l'on peut 
considérer variant dans spectre continu de valeur propre À. Pour cela on est aidé par l'existence 
d'une transformation de Hankel qui garantie le développement de toute fonction ayant les 
propriétés adéquates en un couple de formes intégrales dénommées couramment intégrales de 
Fourier-Bessel. 


La plupart des exemples donnés ci-dessous sont tirés de deux ouvrages d'exemples et d'exercices de 
Physique Mathématique : N.N.Lebedev, Special Functions and their applications, 1965 ainsi que 
N.N.Lebedev, I.P.Skalskaya et Y.S.Ufliand, 1965, < Problems of Mathematical Physics », section 2 
« The Hankel transform » pages 160 et suivantes. Les énoncés ont été parfois modifiés pour ne 
faire apparaître que la formalisation mathématique du problème aux limites de Laplace, sans se 
préoccuper de la forme exacte des grandeurs physiques en jeu. On utilisera également la 
transformation de Hankel dans la section dédiée aux systèmes de coordonnées paraboloïdal ou 
parabolique de révolution, en s'inspirant également des exercices proposés dans cet ouvrage. 
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Le théorëme suivant a lieu : 


Théorëme de Hankel avec les fonctions de Bessel de premiëre espëce 


Pour toute fonction f(x) sur l'intervalle [0,+ce) dont la variation est bornée sur tout intervalle [0,x0], 


et telle que l'intégrale est finie : je | fo) x <œ 
0 


Cette dernière admet la représentation intégrale suivante, dès lors que le paramètre v > -1/2, aux 
points de continuité de la fonction : f(x) = far АЈ, (Ах) [а tJ (AD f (b) 
0 0 


Au point de discontinuité de la fonction, on écrit : 


2 


Cette formule intégrale peut se décomposer afin d'introduire la transformée de Hankel de la 
fonction f(x) : 


f(x-0)-f(x-0) _ je AJ, (Ax) [dt tJ, (At) f (t) 


F (2) = [at tJ (At) f (t) transformée de Hankel 
0 


f(x) = [da AJ (Ax)F (A) transformée inverse de Hankel 
0 


Une forme équivalente est obtenue par un changement de fonctions : 


g(a) = fa VAtJ, (At)g(t) transformée de Hankel 
g(x)=Vxf(x) et #(А)=/АЕ(А)= 0 


+оо 
g(x) = [а VAJ, (Ax)S(2) transformée inverse de Hankel 
0 


Transformation de Hankel avec les fonctions de Bessel de deuxiëme espëce 


Un développement intégral existe utilisant les fonctions de Bessel de deuxiëme espëce et les 
fonctions de Struve pour l'inversion, avec les mêmes conditions imposées aux fonctions f(x) (la 


+œ 
variation est bornée sur tout intervalle [0,x0], et telle que l'intégrale est finie : Í dx| f Vx <œ): 
0 


fŒ) = | dA МАХН, (Ax) | dt JAY, (at) f (b) ee jo dene 
0 


Y,(At) fonction de Bessel de deuxième espèce => К 
H,(Ax) fonction de Struve f(x) = [а [АхН, (Ax)F(2) 
0 
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Une forme équivalente est obtenue par un changement de fonctions : 


g(1)= Í dt t gY, (А) 
Кс. а и 0? 


JA 


g(x)= faa H, ахл 8) 
0 


Transformation sur l'intervalle [x;+°e) appelée transformation de Weber-Orr ou Weber- 
Dirichlet : condition homogéne de Dirichlet sur le rayon intérieur 


Un développement intégral existe également sur l'intervalle [Xœ ^9) utilisant les fonctions de Bessel 
de premiére et deuxiéme espéce, avec les mémes conditions imposées aux fonctions f(x) (la 


+œ 
variation est bornée sur tout intervalle [x0,x1], et telle que l'intégrale est finie : [ dx| f olx < œ» 
0 


pour des fonctions f(x) s'annulant en xx, : 


E. AXE (A,x,xo ) T 
f(x)» [а TPS ETES: | йт*Ф,(Ал„җ)/@) pour x>x,>0 
Y (А, х,х,) = J, (Ax)Y, (Ax,)_Y, (Ах) 


v v v v 


(Ах) => Wi u) = 0 
(Ax) fonction de Bessel de première espèce 


v 


Y (Ax) fonction de Bessel de deuxième espèce 


v 


F(A,x,)= [а t (A.t,x,)/ (0) 


Xo 


e 
n AW. (A, x, x) 
f(x)» [ал X ;F(A.x, 
lët Gr axo «a axi Til 
Encas de discontinuité но) ло) faa ша ( А o) 
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Transformation sur l'intervalle [х„+°°) appelée transformation de Weber-Neumann : condition 
homogène de Neumann sur le rayon intérieur 


Un développement intégral existe également sur l'intervalle [хо,+°°) utilisant les fonctions de Bessel 
de premiére et deuxiéme espéce, avec les mémes conditions imposées aux fonctions f(x) (la 


+œ 
variation est bornée sur tout intervalle [xO,x1], et telle que l'intégrale est finie : Í dx| Рох со} 
0 


pour des fonctions f(x) dont la dérivée première s'annule en x=x, : 


А\Р,(А,х,ху) n 
x dt tY (A,t,x )f œ) pour x>x,>0 
1^ G Gs + Œ Gs 1 | i 


ra x )= J (Ax)Y,'(Axg)-Y, (аә), (Ах) avec Ј '(x)= EI 


Bs ду 
J, (ax) Jonction de Bessel de première espèce 


Y (Ax) fonction de Bessel de deuxième espèce 


F(2,x,)= Í dt tV (A,t,x,)f (t) 


Xo 


AY (2,x,x,) 
dA v 
з EAN ETATS) 


F(A,x;) 


En cas de discontinuité SIE +0)+ f(x-0)) = IL 7 DEM У; 


F(2,x,) 
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Transformation sur l'intervalle [хо,+ее) appelée transformation de Weber-Robin: condition 
homogène de Robin sur le rayon intérieur 


D'après l'article de R.K.M. Thambynayagam et T.M.Habashy, 2003 «A new Weber-Type 
Transform », un développement intégral existe également sur l'intervalle [хо,+еә) utilisant les 
fonctions de Bessel de première et deuxième espèce. Un autre article de 2009, plus récent de 
Z.Xinhong.T.Dengke « А generalized Weber transform and its inverse formula » en revendique 
également la paternité, en se calquant sur la démonstration de l'article original de E.C.Titchmarsh 
de 1924 « Weber's Integral Theorem », formalise correctement la transformation de Weber-Robin 
généralisation du problème de Dirichlet et de Neumann, notamment sur les points de discontinuité 
de la fonction représentait intégralement. Plus tôt l'ouvrage des édition Mir de 
V.Ditkine.A.Prudnikov, < Transformations intégrales et calcul opérationnel» évoque une 
transformation de ce type en page 71, avec un renvoie de référence qui s'avère hélas impropre 
après recherches des articles de J.L.Griffith. La première mention que je trouve de cette 
transformation est dans l'article de 1932, de S.Goldstein « Some Two-Dimensional Diffusion 
Problems with Circular Symmetry », en page 87, formule 227. 


Les conditions de convergence imposées aux fonctions f(x) sont les mêmes (la variation est bornée 


sur tout intervalle [xO,x1], et telle que l'intégrale est finie : je f Go x < оо), pour des fonctions 
0 


f(x) dont une combinaison linéaire de la fonction et de sa dérivée première s'annule en x=x, : 


4 AP, (A, x, x,) b 
LE Gare) Gei reit" e oto 


pour x»xy»20 


14 


Js (ax) Jonction de Bessel de première espèce 

Y (Ax) fonction de Bessel de deuxième espèce 
P, (45x, x0) = J, (Ax), (Axo )+ ЛАУ, Ux, Jj Y (Ах, (Ах, )+ РАЛ, än, H 
— W (А, ху„х,)+ hAF, '(А, х,„х,)= 0 


ôJ, (y) y'( = 250) 


d ia 296) 
avec J,'(x)= ду á ду Y '(x)= 


ду 


у=х у=х 


F(2,x,)= [а tP (3, x, ) f (0 


Xo 


Encas de discontinuité 


se) 7-0) faa : 
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En posant h-a/6, il vient : 


CT av (a) f 
т 1 Geteste Gs ЈА 0060 
af le ji Bf be) _0 E Ce fonction de Bessel de première espèce 
Axo )+ BY, (Axo) Y, (Axa 7„'(Ах„)+ BJ, (Ах )}) 
03-240) уш 0) 
" Ge w '(x)= SL 
ой, (Axo)+ BY, Qs) eA, Qs Ma Ү,'(Ах,)+ В Y, (Axo )}- Y, (Ax Ma J, Dos) BJ, (Ах) 
+ BU, äs, Mäe Y, '(Ax,)+ B Y, (Ах, )у—Ү, (Ах, )\й J, (Ах) B J, (Ах) 
= aA B( Y, (x, M, (ахь), (2x0 Y, (Ах) + Baa (Y, (ахо), (хь) 2. Us Y (2x) 


fonction de Bessel de deuxième espèce 
V, (Ax) = (7, (a Y, '( 
a 
, '(x) ES dal 
y 


+œ 


= Í dt tS, (At) f (O) 


"d Ap. (Ax) 
f(x)» |а I F(A) 
| (шол tele BJ (r, + Get Ux) FX Gs 
Et l'on retrouve les deux transformations intégrales Weber/Dirichlet (a0, 621) et Weber/Neumann 
(a=1, 6-0). 
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Intégrales définies obtenues par la transformée inverse de Weber-Dirichlet 


Р(х) [dt e ss) 
P Ax) = Ax) (Ах) Y, (Ах), (Ax, ) xn i 
f(x)» [da 


A, (A, x, x,) ( ) 
(x) + (Y. x) 


Prenons la fonction puissance: f(x)=x" La condition imposée pour la convergence de la 


transformation intégrale est notamment la finitude de l'intégrale : 
3 


+œ +00 1 

DES 200a 3 

[а= [dit 2 «o ите ? =0= u<- 
t—>+00 2 

2% 2% 

Nous allons voir que cette condition est trop restrictive et que la transformation a lieu également 

pour puissance inférieure à -1/2. Soit donc la transformée de la fonction puissance : 


1 +00 +00 
DEET elei [arer s 0-2 s) ја об ol 
AXo AXo 
Et utilisons les intégrales indéfinies des fonctions de Bessel en relation soit avec les fonctions de 
Lommel de première espèce, soit de deuxième espèce : 


(dz I, o= dev -I Osr -Ia Os] 
f dz z'Y (z)= EI +v- 1) Y, ()5 1,1 (2)-Y,. (2)5,, (2) 
[dz 227,2) = lä v MEM -I,a Sye] 


jezx.o-Aav-06os,..0-x.05,.0]. 


voir :A. Erdelyi. H. Bateman « HIGHER TRANSCENDENTAL FUNCTIONS VOL Il », page 90, 
formule, G.N.Watson < A treatise on the theory of Bessel functions >, Cambridge University Press 
1944 en page 350, du chapitre 10, section 10.72 < Functions expressible in terms of Lommel's 
functions >x, A.P. Prudnikov. Yu.A.Brychkov.O.I. Marichev-Integrals and Series - Volume 2 - Special 
Functions, section 1.8.1, page 37 formule 3. Les fonctions de Lommel de premiére espéce ont un 
comportement asymptotique presque similaire à celui des fonctions de deuxiéme espéce, à ceci prés 
qu'elles oscille fortement autour de la position asymptotique qui est une fonction puissance du 
paramétre u-1. Dans ces conditions il est nettement plus pratique de prendre les fonctions de 
deuxiéme espéce dont le comportement asymptotique et plus « régulier ». Ici prenons donc les 
fonctions de Lommel de deuxiéme espéce, qui sont en relation avec les polynómes de Neumann, 
pour des valeurs particuliéres des paramétres : 


s; (z) Lommel première espèce 


S, (z) Lommel deuxième espèce 


1 
O, (x) = p 1 
1 4 x Och: х2 
S 2» (x) -XO,, (x) avec 30, (x) LES В )= 2т+1 (x) avec 
X x 2m+1 O (x) ü E dc 24 
1 16 192 QAAE T. + 
O, (x) SRE ias” " : 
X ` X ` X 
5,(х)=0 S (p) =E 
50) 5, (x) pot 9 2 x? [O (x) polynômes de Neumann 
X ] = 
он 4т 8 96 S,(x) polynómes de Schläfli 
S, (x) = m + SC 
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Pour ces valeurs du moins, on observe le comportement asymptotique suivant : 
vm meN LimS,,(x)21 , LimS,,,1(x)=0 , LimS ,,,(x)-0,. ce qui est en accord avec 
x->œ0 i хоо d x->œ0 i 


le comportement asymptotique connu des fonctions de Lommel de deuxiéme espéce (en fonction de 
puissance) à savoir : 


uc) (=) -v ens EEN “е 


2 4 6 
X X X 


Dans ces conditions de paramètres coincidant avec les polynómes de Neumann, l'intégrale 
indéfinie s'annule toujours en +œ pour u«-1 (et plus précisément u<-1/2), comme suit : 


Ze ud -Ax J, (Ах), (А) - Qu v) 7, (ax, )S SC 
Da SEN бъ) | 


n u+ — Ax У, (Ах), (А) (u+v)Y,( uh SUE 
d Sal? D L | 
SITE 

=> (g. Ab) о) av), ich, Ac 

"m G dë Us. (4x )S antv (4x,)- (u +v) J, (Ax, Bayi (Ax, ) F 


ES ur 
n] nf Cu) an) Y Dh. as) | 
а a xU y Se (Ax) 2, (E. (Ax 


+ Ax (и +у)5, Qo y, (x, H. Qo) - J, (Axo ) Y, (Ax, )}+ 

— А5, x JY, Qo, H. (Ах), (Ах Y. (Ак) 

Comme S, ax) (Qo) S, a, (Ах AE (ax) 

= А(х,)= Ax (u-v)S, äs NX, (ax), (Ах, )- J, (Axo ) Y, (Ax, )}+ 
= А5, „(Ах NY, Qo, H. (Ах), (Ах, Y (А) = 


(AW. (aa) УЛ Gr, 


v v 


F (Ax, yel- Ir (Ax, Ja (Ах,)- J 
1 el, 1 

Allez? Y(x)ex 2 2 Дх e x" 

Xos Lay (4x) 


1 
Sí и<-1= Lim A(x.)=0= Р) Der: 
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En calculant l'expression sur Mathematica par les substitutions adéquates, оп parvient d 
l'expression asymptotique suivante de A(x1) qui confirme la validité d'utilisation des intégrales 


définie : 
DEE y PA (Ax, X, (Ax, )- n (Ax, у, (Ax, UE 


А(х,) = É д7 дү, 
Fer tee Am ESP 
x=Àx, x-Ax 
1 1 1 1 
J, (x) œ x 2 Y (x) œ x 2 ^ (x) œ x 2 one) ox? 
x x-Ax X х=Ах, 


= Alx,) œ gh = Lim А(х,)=0 pour u< - 


Il vient donc : 


| 1 +00 > 
Si H <= 2 = [at і" I J, D = Ах, Ua (Ax, 18,235 (Ax, )- (u * v) J, (Ax, ы (Ax, } 


2% 


| 1 +00 d 
Si p<-7> [ate Y, (e) = АД, (ах, ) 5а (Аха) (и +v) Y, (ах), (Ах) 
Аху 
DECH 
tal бене Pt Gs yr. (as) 7, Qs, (sJ 
Prenons maintenant les formules qui relient les fonctions de Bessel à leur dérivée, ainsi que la 
valeur du Wronskien des fonctions de Bessel de première et deuxième espèce, il vient : 


эы = da. (4) 72 (A) EG) г ar. е (J 


= Y, (4x, J, ASS (Ax, ы (Ax,)= Y, (Ax, 


—J, (x) 2—0 


14 


3J, (Axo) 
) Ох 


"m J (a VH) y (у) AX) _ 24 _ 2 
puis add EE. TAX) Хул 


Se | ( )#Ь@%)_, (a) 28 puy (Àx) 


=> Е(Д)= 


Y (Дх 
SE VV & : дх aA"? 


Formule de la transformée de Weber-Dirichlet de la fonction puissance : 


28 „a. EEN ТОМ. 


url 


E dt t^ (Y (а) (ax), (20), (Ax,)) = =a 


Xo 


Il vient par exemple а. v=1: 


Sus) o FQ) -fat J-J, (A, Ax) - —— 


)= [art (4, (2x "AG ty, (Axo) $ 2 
Е л XV x, 
La transformée inverse de Weber-Dirichlet donne dom l'intégrale définie suivante : 
vu < à f ag es P QS Gn) 4 GU) mt si > x 
2 74 A" (J. (Ax, + (Y, (Ax, )) 2|0 si xx; 


v 


Et pour u=-2, v=2 : S, (Ax,)= = 
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Cela donne par exemple dans le cas particulier oü u=-2 et u=>1: 


x> x = [ал Su 0x) Y, (Ах), (Ax,)- J, (Ax)Y, (Ах) _ Zu 

je essa 72 
Ө n Y, (Ax)J, (Ax,)-J,(Ax)Y, (Axo) | « 
ver T ex "Se 
SE Deel адаы) sl 

IC US QR 2 

eT Y (ax), (x), (a) (Axs) л 
u= > а (A) TG) 2s 
A _ 1 ` dÀ Y (Ax) (Ax,)-J (ax) (Axo) _ T xo 
RE |а иб) Б. 


Cas v =1,v = 2 < Prudnikov.vol.2.formule.2.3.17.1.page.275 
j Z Y, (Ax)J, (4x, )- J, (Ax, (Axo) ` 263] 
o 4 ` I, (Ax) +(X, Ms 2 


v 


Vv 


x 


v 


Il se trouve que la formule obtenue marche également pour le cas u=0 (même si ce ne sont pas les 
conditions formelles d'application de l'intégrale indéfinie). 
x»x 4-0 v-0o S, (âx,)=1= faa YA). Gs) (х) Gi) =” 

» À (J¿(Ax))+(Y(Ax.)) 2 


Cas v 20 € Prudnikov.vol.2.formule.2.3.17.1.page.275 
аА Y, (x, (x es (2), (Ах) T ( xo y 

Vv Í 2 > = | | 
00 QJ) + (x) 2\ x 

Au passage l'intégrale de Prudnikov représente une transformée inverse de Weber-Dirichlet et nous 

obtenons l'intégrale définie suivante en inversant la transformation : 


V.X> X T À Y, (Ах), (Ах„)—./,(Ах)Ү, (Ax) _ z Xo 1 
ve Jos SD e AIS 
=> fa x( 22 ] at 0) - (or Qs) fa "lacht (a) (ar (Ax) —— 


En réalité cette intégrale est facile à calculer à l'aide d'intégrales indéfinies trés connues, il vient : 
Гас) a) fae T eTa) 
X—+o 


1 
xJ, xe x? => Pour v > ; Lim x^" J, ,(x)=0 


l-v LA 
IL x" J (Ax)= An) | dx xY, (x)= TX n) 


1-v 1-v 
_ Xo (7,0, 98) (ax) El 2х, > 2 
x Ох 
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Comme оп a démontré facilement la valeur de l'intégrale définie précédente, en appliquant la 
transformée inverse de Weber-Dirichlet, on démontre la formule de Prudnikov : 


ER 72 Y (Ах), (ахо) J, (Ax)Y, (Ax) _ SÉ 


2 vi QS) 2 
Pour toutes les valeurs de paramètre ой sont établies les correspondances avec les polynómes de 
Neumann , il vient pour u=>2 : 


Si; (Ax,) = Jak, (ax) Р (Ax, 


x 


) == AX On (Axo) 


2т +1 
S T S, (x) vs, (Axo )- Jan (Ax)Y, (Axo) _ T 
EE =) Ра O 2 
Z T dA 12 Y, Ax, (Ax )=J,(Ax)Y, (Axo) _ m Xo ; 
x Gex) + (v, (ax) SE 


Фал Y (ax, бах) J Ax) Ax) т 
salaam о a 
j Y, (Ax)J,(Ax,)_J,(Ax)Y,(Ax, 


A Häff 24" 


м“ 
lI 
a 
[=] 
N 
P SNR 
= Lë 
SC 
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Et pour les deux cas incluant u=-1 u=0 : 


Ax, O Ах, 
S, 2m (Ax,)= Ах,О,„(Ах„) So ma (Ax,)= 2500) 
Y, Ax) am. (Ax,)- Jut (axe) zm +1) 


(J PR E Ae, +(Y,,, (Ax Axo) 2хх, 
je 


x>x, H=-l v= m cs Í dâ 20, (x) 


rasa жЕ DA) USQUE a) m 
| 17 í >) (J: Ax) + (7, (2x ))° 2 x 
Comme | Y (AxJ (Ах)— J (Ах) (Ах, EE 
A (J(Ax,)) +E (ax im 
рад) а) _ z x ( 3i 


0 A (Z (Ax, ) us (Y, (Ax, ) 16 x 


_ Е r Y, (Ax)J „m (Ах) Ji (Лә) (Аа): T 
x>x, u-0 у= JO on) US С 2 
Uic з= 4 lest Ee = 
À < (7, (Ax, ) + (Y, (Ax, ) 2 
ra Y, (Ax)J (Axg)- J, (Ax Y, (Ax) _ z Xo ў 
2. o A. 


А — Uy +(У,(Ах,)) 


Ge 1> 1: 16 , 192 Y, (Ax, (Aa =Z (AAE) _ x 
WE AA 


ra Y, (ax JJ, (Axo )- J, (Ax), (4x, ) Л 2 X 
| 8 Xo [ x? | 
( 


Or jÉ e S EET 
o À (x) +(Y,(Ax,)Y 2\ x 
| 12 меме мыны. л al s 
A el Wëss RN 
Фал (ax), (x, )-Z (x) (А) 7 fi „ү 
game > 20 [ (e) 


о «ax 38 


) 
em ET 
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Intégrales définies obtenues par la transformée inverse de Weber-Neumann 
F(2,x,)= [аг (a,t,.x.)f() 


P, (2, xs) = J, (Ax) (Ax) =Y, (Ax), 0) e WW 
х)= [da 


аздык) ra) 
(J, Qo) +(X, (a) 
Prenons également la fonction puissance : f(x)=x"- La condition imposée pour la convergence de la 
M oemanon intégrale est notamment la finitude de l'intégrale : 
+оо cd 43 
3 
[are Jt = ac^ 2 <o> Limt' 2 =0=> u <-> 
2% 2% Lui: 2 


Nous allons voir que cette condition est trop restrictive et que la transformation a lieu également 
pour puissance inférieure à -1/2. 


1 +00 
Gex = Pal ze јао -s a fac х0} 
Axo x 
En utilisant les intégrales indéfinies des fonctions de Bessel en relation avec les fonctions de 
Lommel de deuxième espèce et dans les mêmes conditions d'application de la formule de l'intégrale 
indéfinie (soit и <-1/2) , il vient : 


jacent a) ax, US, s) Qv) I (x) Ae 


2% 


ee | ad s... (A) BE {к rer às) 
Y (Ax ACTE (Ах QX (о) _ 22 _ 2 
v 0 дх v 0 дх TAX, s 


DIE del DO 40-1763) ja) 


Axo 


Fx)» Â =. 
ЕЕ ЕСИ 
= F(a, e )= 1 E Abe (Ax SIS EL )- Y, (Ax M, (Ахо))+ | 
ж +Ax(u+v)sS S ud ax Ax, MY, (Ax,)J, ax 9)- J, (Ах, D. (Ax, H 


= x) CO 88, a) V Sua Qo) y (A) (ax (A) 


quel S 
= Р(А,х,)= (и +у) Ах, A P 5, (Ах ty (ax Gas )- б) Qs, 
F(2,x,)= — V S, 4, (Ао) (u +v) S, (Ах) 


xc 2⁄2 
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Toujours dans G.N.Watson < A treatise on the theory of Bessel functions >, Cambridge University 
Press 1944 en page 350, du chapitre 10, section 10.72 < Functions expressible in terms of Lommel's 
functions », il est donné les formules des dérivées des fonctions de Lommel de deuxième espèce : 


OS, (x 
5.6) | s bu (x)+ (и +V- UR, 
= Ce E =-* S igo (x)+ (u * у), (х) 
x X 

08, Ax І 

+ zt | = ИУ) 
X-X 0 
05 (Ax 0S, (Ax E 
=> v Spay (x) Ax (u +v)S, „-i(Ax,)= x, S 2) ml Е r 
| ) x |. X) |= D bu 
D'où le résultat assez court suivant : F(A)— 2 Pin O) 
ТА" ду EEN 


Formule de la transformée de Weber-Neumann de la fonction puissance : 


ође Ja” U, (ah, as) (a ne 2 Buse 0) 


La transformée inverse de Weber-Neumann donne donc : 


Rum EN w 737 
A" Oy у=Лх (J, '(Axo))" + (Y, Län, 0 si x«x 


Ou encore plus exactement avant de passer à une limite donnée : 
Y (Ax eJ, (y) -J (Ax) дү, (y) 
| Ор Ae, lise Е H si X> Xo 
2 2 А 
эй êJ, (y) К er (y) 0 si х<х, 
ду у=®ху ду у=Аху 
Si maintenant on pose x=x0, puis x0=1, il vient : 


-J (Ax, aY, (у) M Y 
y=% ду у=Аху TX, 


Yy=ÀxXq 


2 + 
xis 


— 2 [da 1 ëch 
T 


qrt! ду 


Comme Y, (Ax, E 
ду 


k | 
E 4 j E 2 3 у=Аху —1 
7 êJ, (v) ‚| 2x0) 
Oy Yy=Àxq ду у=Аху 
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Soit donc l'intégrale indéfinie logiquement valable pour les valeurs u<-1/2 : 


1 Sua, À) m’ 


А QG) Дуу 4 


Pour des valeurs particulières des paramètres, les fonctions de Lommel sont associées aux 
polynômes de Neumann : 


иэ faa 


1 
1.4 x 06)= == 
S, (х) = x0, (x) avec О,(х)=—+—; MN )- belt) avec 
X x 2т+1 О пт ы 
1,16 192 PAT 527 а 
O,(x)=— 73 T D * UR 
A ME 
$()-0 $()-7 
SG) 85 (x) OUT 2\X) > O,(x) polynômes de Neumann 
x)= avec 
ge 4m 8 96 S (x) polynómes de Schläfli 
S,(x)=— Res 
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Il vient les diverses intégrales définies dont on remarque qu'elles s'appliquent également au cas 
u=0 : 


Р di 1 


1 ‚ 2 _ Фал l x 
po vota s) Ss 0-5 [GET s 
u=-2 vom Su) => Sa) 
Sab. 24 1 NN JPE 1 л? 
dk TE s ^1 а ay 1 
| 1 _ T 
74.2 SOY QUO] 384 
1 | 1 фал 1 л? 
шы уй бл о 
1 | 8 Фал l л? 
и=0 у=2>5,,(1)=2 EE FE E: GI DG 32 
и =-—1 v=3> SS) Le 5 (2) = 2 
**da( 24 1 к К к 1 E 
E Can zan pn BUGA) 12 


717 AMECA 
EE E Pos e J +22) 

“dA(. 24 1 xi. 0. o 1 л? 
= fen Zem RO EE say Gray] зи 
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Pour toutes les valeurs de paramètre ой sont établies les correspondances avec les polynómes de 
Neumann, il vient : 

Р : | OQ, al) A0, E 
S ол (4)= o,,(A)+ 4O,, (4) KÉN (A)- > 1 ) = ( ) 


2m+1 
1 ENT ER л? 


Сотте faa 
0 


и=-2 у=2т= [a4 | 
4 


) 
=z roy «Gy, 


u=-1 у=2т+1= (dA 
0 


* 1 O,, (A) 
=0> [da = 
u-0 u A 


2m 


Dans ce qui précède on a utilisé la formule suivante : 


(dil 1 AES 


S PANVE aN, 4 


Étudions maintenant la transformée de Weber-Neumann de la fonction puissance négative x": 
[des All aa) [dee T e (а) 


1 
xJ, \(х)сх? 2 Pour v > = Lim x'" J, _,(x)=0 


X—+o 


1-v 1-v 
| dx x" J (Ax) = X Tal) fax xY (x)= x i 
J, (Ax) = талх) у y (ax) Y, EE 18Y,(Ax), v y (4) 
À € Ax À Or Ax 
1-v 
=> ar (lt 0.) о r D= C, QI) - 1G) 
ye v Xx v 2 2v 
=” - (J (Ах, YX, (Ax, )- Y, (Ax, JJ, (Ax, )) = >= = 
À Se v ( o) „( ï) v ( °) A 3) ER Ах, T ху л Ах l+v 


Sea Gr, Qu Qs) - 3, Gr (Б йз” 


v v 


ы x л Ax 
Il vient alors par la transformée inverse, l'intégrale définie suivante : 
p (J, (Ax)Y, (ахь) -Y (Ax)J, (ахь) x 


а Об) оО) mex 
2 qa 1 T x 
C J (Ax, Y, (Qo) - Y, (4x, J, (x) = 
omme „( 3 v ( 3 it 0 v ( o) Л Xo Б | А, (J, "Us, )) + (7, '(Ax,)Y 4v 
dA 1 EU YÉ. : А 


aO) | O) 
Oy Voas Oy 


On a donc démontré la formule utilisée pour établir les diverses intégrales indéfinies. 


à E ES 


| (A V5 (ДУ QUO 4 
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Transformation généralisée de Weber-Orr vers l'arrière 


D'après l'article de C.Nasim de 1989, « Associated weber integral transforms of arbitrary orders », 
une transformation intégrale généralisée de Weber-Orr est possible, et construite selon la formule 
suivante : 


F(A,x,)- Га xY. (4,x, x, )f (x) 


V „(%.,х,х„)= J, (Ax), (Ах) Y, (ax), (Axo) M i y (À ) Е 
FO Ја еу ^ 79 


Plus précisément l'article étudie Іа transformation intégrale ауес le noyau suivant construit à 
partir du paramètre u-—v-a : 


к(а) [as х%,.„„(д,х,,)г 0) 


Y, (охх) =, (Ах) (Аха) Y,a (А), (А) Eo We ue 
Jode [a cr as (кай C) 


L'inversion y est prouvée lorsque les paramètres sont définies tel que u=v-a, ауес: 


бое Чы. si aeN—0«a« nm 
2 4 2 4 


ыш E i 
2 4 


+ |=20<a<2=a=0,1 


3 
4 


[v 
2 
Iv 


"n =3=>0<a<3= c = 0.1,2 
2 4 


Transformation généralisée de Weber-Orr (cas u=v-1) 


Y 


v-Lv 


(4; )= DU (4x, SY. 


Za: Joa XY (Ах) Y. (4x), (ax) 
(kl: Jan тоол сш. 
Autre formulation 


ы '(Ax)=J 


À ôx d 

En notant J, (ax) 9,0 
ду 

(Ах, хо) = J, Qr, (Axo)=Y,'(Ax)J, (Ах) 


=> Y 


v-Lv 


Notons que See EE (Ax,)-Y, Cox, M, (Ax, )=- Z Wronskien 
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SÉ f, qtd] а) 
STEEN 
De plus Al EA At (а) v(a = 2 y, (д)- a) 

X dr б) sn e. (A) (A) (A) 
SUE 


Rappelons les intégrales indéfinies suivantes : 


faz zu (2) = EVI +V- 1)J, (2)5, ,,4)- J, (2)S, , (2| 


S, (z) Lommel deuxième espèce 
[dz z*Y (2) = (А +v Ate, tt У, eg, ell 


Calculons la transformée intégrale de la fonction puissance x": 


[at po di (t) = Ах, (J, (Ax, dE (Ax,)- (u +V – 1) Ja (Ax, )S uv- (Ax, )) 


Аху 


[ar D (t) = Ах, | (Ax, Jis (Ax, j= (и du 1) Y, (Ax, jus (Ax, )) 


2% 


Fa ут үка", i aeo] 


F(A )- Ах, Y xU. (Ax, Жа (5-6 +v -1) Jua (Ax, )$ 22 (Ax, j- 
ш Е A - J, (Ax s (Ax, )S ava (4x )- (a +v -1) Y, AR )$,„_2(Ахо } 
EN FU )= AU y, (AJ, (45) J, (Ax, ўҮ,_,(Ах)}+ | 

dM - (u +V -1) KEE ЙҮ, (Ax,)J, (Ах). (Ах, Y, Dal 


v 


2 Í 2v-1) (uv -1)85,, (Ах) 

=> F(A,x,)= — | S ava x9) u, 

V Pç L PRENS A 
EE СҮТ, тиру fax, (и +v —1)$, „_›(Ах„)—2(у AS, (Ах } 

0 

> 2b (u+v AS (Дх, )=29 - DS. Dy 

=> r(x)= [axe (d, (a Qn) 1 (a, (x) E Ut VS as n 2G D (te 
0 


Xo 
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D'après l'ouvrage de Watson sur les fonctions de Bessel et relatives aux fonctions de Bessel, nous 


ER hX 
ER = х(и +V -1s, a EE Bus (x) x E = х(и zk 1s, Jet: (v m US. al 


05 х 
> x(u VE ns (x) ES 2(v Zei TE (x) =X vd (v ER (х) 


Сотте 50 =V S s (x) = х(и КЕЙТ 1)$ gc) 
x 
+ an (х) (v 1) S atv (x) = x(u —V + 1)S,, (x) 
x 


> x(u +v— Does (x)- 2(v = DS ал09) = х(и Sc 1)5,, (x) 


On arrive encore une fois à un résultat fort simple : 
2(u -v +1 
F(A,x,)- e Ds, (ax) 


F(1,x))- [dex p cb (x) Y. Delt s ТЛ» (a) 


' л Д 
— x> x, [aa Sy (2%) JG), (0) - (Ах), (дж) — E x" 
0 A (J, (ax) + (Y, (x) 2 u+l-v 


Prenons maintenant la valeur х=хео, il vient : 


2 + dA S. (Ax,) n? x "t 
J. (Ax )Y (Ax )- Y. xd (Ax) = _- <> ы = S — 
m ^" qi o) ge UL o) x TÀ => [ Д? (J, (Ax) + (Y. (2х,)) 4 u+1-v 


y 


Au vue de l'intégrale la valeur de x0 est non essentielle, et l'on peut prendre la valeur x0=1 : 


Í 4А $209 a 1 Í dÀ S, (A) __ z? М 
2⁄2 QG) + CU) 4 u+l-v А (JAY «Quay = 4 u+2-v 
__1— [45 S, (2) ES о À S. (À) _ T 

da | ОО) w ^ ‚ | X (GA) +(X, 4(-1) 


Si l'on se donne comme paramètre ceux pour lesquels les fonctions de Lommel de deuxième espèce 
coincident avec des polynómes de Neumann ou de Shläfli : 


1 
Oo (x) == 1 
1.4 x а 
S 2n (x) = xO,, (x) ауес О,(х)=—+—; sl )= Safe ауес 
X x 2m+1 3 24 
1 16 19 Ox)- Zen 
O,(x)=—+— 5. 
bd x x 
4 
S (x)=0 S,(x)=— 


x O,(x) polynômes de Neumann 
S, (x) polynómes de Schläfli 
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Alors on peut écrire : 
x = dÀ l om 
A A (76) Qua s 


SONUS NET NV. 1 л? 
=-1 v=4 S,,(4)=—|1+— 1 Ü 
u v a4(4) H EEN a lageren 16 


ial s,ex(« Jo FS)" s 


и=-1 v=2 S ,,(4)- 


dÀ 1 р | 
А (s, (a) + (yQ) 8-1) 


+оо 
Plus généralement оп constate que : | 
0 


Pour démontrer le résultat de cette intégrale définie, calculons la transformée de Weber-Orr d'une 
fonction puissance x^" : 
[desc om Ea ЩЩ cs e) so" a) 
3 


x*" EE x? = Pour v > i Lim x^" J, _,(x)=0 


X—+o 


2-v x us (ax) 2-v x” Y, ., (ax) 
[à xJ, (Ax)- = qucd | ax x^"Y, (x)= B 
= fax xt (Ui (ax), (2х, )- I (Ах), (Ax, )) = => [do (2х, Y, (4x,)- Y, 2 (2х, V, (Ax, )) 
Comme Y, (ix. (о) бз) =) 
> jac Ga (ans) y (Arr, (âx) D 


Comme on a démontré facilement la valeur de l'intégrale définie précédente, en appliquant la 
transformée inverse de Weber-Dirichlet, on démontre la formule suivante : 


"t dA J, (А), (Ax, ) У, (ax), (Axo) л Ж 74 , 


dä (¿xY GS) => 4 (v-D) 
Cela donne pour x=x0 et x0—1 : 
3 Фал 1 л? ху dA 1 л? 


"^3 [ж оноу взу) Too %-) 


0 


3 
Vv»— х>х, 
2 
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Transformation généralisée de Weber-Orr vers l'arrière (cas u=v-2) 


F(A.x,)- )= (dx x( J (А), (Аха) (x, (x) f (x) 
e Jual), (x) Y, QJ. x.) j 


| (¿Q P Ets О 


v 


Оп a besoin de quelques formules, comme suit : 


LD (n) =Z s Ga) BLE = у, (Q3) 0) 
A (a) 20220), Mad (a) Y, Most 292, аз) y as) 


> Y, (Ax, л (Axo )- J (Ax, JE (Ax,) = 


v 


26-26 (ax y, бахо) (axe) sU) C, C a) о As) 


x) [are a.) (ax) arr" X. o 


Е Ах, о = 
-J (x У, äs. S) (uv - 2) v Y, „(Ах Ss (x } 
Xo cu [xu ;(4x,)- J, (Ax 118 (Ах DIE 


—(и+у—2)$,„,_ „(Ах „ЛҮ, (Ax, Ve „(Ax 0) J (Ax DAE „(Ax sch 


lh aer xo ca У 2b sm) 22 -1Йи +у - 25, (an) 
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D'après l'ouvrage de Watson sur les fonctions de Bessel et relatives aux fonctions de Bessel, nous 
avons : 


@) 2v $,, (x)=x(u+v-=1)S„, аб) XQ -v S, ca) 
(2) 5,.2„(х) x" (uv +u -v € US, (x) 


6) «вуз, (= uv RA 


=y S uy (x) = x(u == 1)S pava (x) 


бф Е -10 2) 5, (х) 2x(v - 1X v – 2), (х) 2х(и - Yu -v +2)5,, (ә) | 
=> 2x(v -1)u +v —2)$„„,(х)—- 4(v -2)(v —1) Su. (к) = 2д(у —1Xu-v + 2)$, v(x) 
е. x(u-v-18,, ek Ale UR, a [x)=x(u-v+1)S. (x) 
бешш решен) Geck a DR e] acea 
Sale xet —(и De v+ 2)S,, e (x) 
OS latices, oe 7 rs 
u+v 
Calculons 2х(у 10и +у – 2), ,(x)+ (x? - 4(v =1)(v - 2))S auv 2 (x) = 
lat. ele и -v 295, 2] 


pu. 
= ds. (x)+ 2(v 1) базу 


(uv) 
_ xbv = LSC? + x(u +у)5 ауз (x)- 2(v = Us, (x) 7 x° v “y (u-v+ 
Ge (и + v) e (u + v) It 1)х (и SR (x)! 


Compte tenu de ces relations, il vient : 
(x? - 4(v —1)/(у – 2))5, (х) + 2x -1Xu +v -2)5,, ,(x)= 
= cke, (х) + 2(v = (и SE 2)S d (x)! 


= G =) lv — Dx" + (и -v+ 2) nr (x)! 

Рб) = (о, sche 207 и v 25, lA) 

Р) aq Rv (и 2)8, n) 

аа) aaa, uapa Gaya с n 28, Qn) 
A лл" (и + v) 


En inversant la transformation généralisée de Weber-Orr, il vient ; 


Gro (y) GG) + (v, (Ax, 
Jes (Ax)Y, (4x,)- Y, „(Ax 


v 


) m +v )x" 
Q (фый 2 


v v 


"VES í dÀ 2(v-1)2" x," + (u =v *2)5,4, (Ax) A äs. (Axo)- Y, (x), (Axo) 


+œ dA 
> [acte е (a 25, 0) 
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Transformée de la fonction puissance x°” : 


Calculons la transformée de Weber-Orr généralisée d'une fonction puissance x°” : 
jaa, ux)-— 2a) [dex ux) (5) 


5 
x" J. a(x) x? => Pour v > ` Lim xJ, (x)=0 


3-v 3-v 
keam а [ах "Y, ,()-- xY, 4x) 


> jac обаа) (ауу бао) (о база) Qs) - Y. (A, (s) 


Comme Y (x s) (E s) — [ SE 


14 


> јазот кь) ëch а) 229 [ av -1v - 2 


Xo 


= je Ж” (J, (Ax)Y, (Ax, )- Teo (Ax), (Ax, ) ER 


Xo 


TX Xo 
Par la valeur de l'intégrale définie précédente, en appliquant la transformée inverse de Weber- 
Dirichlet, on démontre la formule suivante : 


5 (44 (12.2 J, (Ax) (Ax) Y, (х) I). 7 „ 2( Xo Y `. 
nä en a oa SA 


Posons quelques notations des intégrales indéfinies suivantes : 


Mes” J, (ax )Y, (Ax) =Y, (Ax), Qs). z (= | 


SW жоол EE 

TO) (x 3 jp J, (Ах), (A )- X. (x 9, (Axa) 

ees а а) ос | 
2 рад J, (Ax)Y, (x, )- Y, (Ах), (Ax) _ л Xo fe 

pes) READS SUE cap us) 


v 


dÉ X je J, (ax), (4x,)- Y, (x) 
UU SV Ga) +( (Ax) 


I$ 


Vv 
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On arrive aux résultats suivants : 


v-2 
Alors x, I0 (x,x,)- 4(v E 1v Е 2)19 (x,x,)= SCH 9 


De plus A (x) 26 (ue) 4 (2) A 

уб nr Par of ce 
steil ax.) 

GEET EE 


Ce qui ипе fois réinjecté dans la formule de la transformée inverse d'une fonction puissance, 
donne : 


"t dA или 7,_›(Ах)Ү,(Ах„)—Ү,_,(Ах}, (axo) _ T и 
х>® [жи (ov 1) 2 + (u-v 2, Lä 5 RE (ижи)х 
ИУ [аА J, (Ах), (Ах) Y, (Ах), (Ах) m (u +v) x" 

0 [=ч S is (0%) (7, 2x, )) + (Y, (ax, Ir 2 (u-v 42) 


Passons maintenant à la limite x=x0 pour la fonction puissance x°: 


(Z J, x), (Аха) Y, (Ах), (ax) _ T = Е N 
dä (J, (ах) + (Y, (Ax, )Y &v-1(v-2) "Vx 
4(v -1) 


Comme Y, (Ax, )л,_›(Ах„)—.7, (Ax, Y, _,(Ax,)= rx 
= e 1 А in 
o À Q,Qx) +(X, (Ax) 320 -1)'(v -2) 


Soit en s'affranchissant de la valeur x0, il vient : 
5 í 1 л? 
À HOT.  3(-)(v-2) 


Уу > 


Il se trouve que le passage а la limite х=х0 ne produit pas de contradiction dans ce cas. Mais il 
n'en est pas de méme pour une fonction puissance quelconque. Comme on va le voir par la suite. 
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Passage à la limite en x=x0 pour la fonction puissance x! : 


Comme Y (xs) 2 (У (a 4027) 
de d^ J, ix JY, (Аха) Y, (Ах), (Axo) % _ x (uv) 
0 | Said) (J. (Ах, )) + (У, Go, 2 0 (u-v+2) 
= | a S uuy (Ax) LU y (u +v) => contradiction 
SA Häss, `` 8 7 (v-1u-ve2) 
* di $,4, 4) m? l 


O stablit t | а 
Кы E Cre QE 7 2 


Donc par ce simple passage à la limite x=x0 dans ce cas, on ne parvient pas à établir la formule 


AQU) +(X, (4) 4 (u-v+2) 
Il vient les formules suivantes que l'on a déjà établit avec la transformation généralisée de Weber- 
Orr du cas йл 1: 


d Su) | ж 
2 (Z, (2) 


+00 2 
correcte. Si l'on part de la bonne formule à savoir : | 42 Suy (4) ij 1 


legi | 
0 


274 PANE s 


u=-2 v=4 540-540 - (1 ]2 | ett 
0 


ЖП GP ay is 
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Transformation généralisée de Weber-Orr vers l'arrière (cas u=v-l) 
Commençons par formaliser quelques résultats sur les intégrales définies que l'on vient de calculer 


Règles de récurrence sur certaines intégrales définies 


Soit l'intégrale définie suivante : 1x, Xo ) = qs Z= (ах), (Ax, ) IS (Ax)J, (x) . 
И Tu (J, (x) + (Y, (x) 


v 


Appliquons lui les règles de récurrence des fonctions de Bessel, il vient : 


200-1 200-1 
дан) D, а) a) т, a= Er, (a2). y (2) 
u-2=v-l=>u=v-1+2=u-1=v-1+1 
2lv -1+1 2{v-1+1 
eh HE „()-у, (09) x (9) Ey | y, ба) 
2lv -1+1 2lv -1+1 
ea, OE, tel DT a) ana, 
аг VIC EN 
0 v 0 v 0 
2v-1+1)*F dA. J, 1. (Ax)Y, (Ax,)-Y, 1 (Ax)J, (x) v +1-1 
x | a (дыў (Ya) Ia = HD m) al 
Se el A 
dA J, olx), (#х,)—Ү, 10 (Ax)J, (x) If? (x,x,) = 2(v-1-1) -1-D rea. x,) - 412. (x,x,) 


VA (J, (Ах, )? + (Y. (Axo))' | 


v 


Appliquons ces règles de récurrences aux résultats obtenus auparavant, à savoir : 

EEN 
» 4 — (0x) +(X, (Ax,)) 2\ x 

(2) 19) (x, х, ) = Í dA J, (Ax)Y, (,)- Y, 2 (xy, (x) 
o 4 (J, (x) + (Y. (Ax) 


v 


(3) 10) (x, x) = (Es Jua (ax, (Ax, )- Y, (Ах), (Axi) Z Е А 


2 OA) +Œ), ` A0 
и) (к) [42 3 0) - AM, An) л (Y 


A (+ „ж 8-и 


v 


а) ьи mY 
x 


je pis on) 2220174: 70022) 
X 


v-2 v-1 

"o ___7 EN kara ¿ge H m 
ceu) tae s 

‘pda J, (ax, (Ax, )- Y, (xy, (Ax, ) 
SS | 2 2 2 

de (J, (o) + (x) 
De ces résultats, nous pouvons par exemple faire l'hypothèse : 
0<leN<v Alte, )- | Е J, (Ах), E Ti (x, (Axa) =0- 

| 0 À (J, (Ах)? +, (x) 


v 
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Dérivons l'intégrale précédente : 


аш) AJ ae), (ан) 913093 ie) - ы) 


At x )- Í dA J,, (Ax), (2х, )- T, (Ах), (Ax,) = 
Тос жоош 
E (darr (R) =" Ge ard 1-1 939-09 о) 
Es PR 


EE 
Le eg NR D 


=I! HNx,x,)=0 


1+1,у 


Donc оп arrive аи résultat suivant : Ix, xy) = Í di Es (ax}r, (4x,)- Y, | (xy, (2, ) 
` de (J, La + (©, Qs) 


v 


Si on applique la récurrence à l'intégrale identiquement nulle on arrive au même résultat: 
(1) DPG у= 21D MESS )— T (x, Xo ) > = ) = I, (x, x) 


u 2(v —1—/ 
(2) D et.) = MM ,(x, x,)- P (234) 
Comme Ix, x,)o I (ж %)= 0 (2)= n (жу) =0— I" "(x,x.)= 0 
Ce résultat reste valable tant que l'intégrande garde une puissance en 1/4. Mais la règle de 
récurrence établit ne permet pas pour l'instant de calculer toutes les intégrales définies de la forme 
proposée. 


Calculons la transformée чо de Weber-Orr pour u=v-l pour une fonction puissance x" : 
[ax об )= =x “J, -(1+1) (x) [ах x" Y. ,(x)==x Ж + FU. ) 


1 
x J. paye Je x "T' => Pour v > +5 Тн. EN )=0 


x—> +0 
vil =v 
J, -(1+1 (ax 
À 


= je p (J, (Ax)Y, (Ax,)- ) (Ax)J, (Ax, ) Ж 


Le produit croisé obtenu fait intervenir les polynómes de Lommel selon la formule suivante : 


2 1 
uquy 4% X, (Ax,)- Y nay V, (Axo) m VP. Reg pu 
0 0 


xy cux) 
A 


D u EE D. (Ax, ) = Y, (ay, M (Axo )) 


| dx x" J axe 


Гаа) 


1+1-у 


Terme dominant J, ухо Y, (Ax,)- Y, qx, ) 


L ++ — 
Terme dominant NP | Se (es e 2) у) > |. 2^ (y (be =-2)---(v =!) 
0 0 T; 7 Al 


La transformée de Weber-Orr généralisée de la EE puissance х!” est donc la suivante : 


+оо ER 2. -=v 1 
NEEN 
Xo 


Xo 
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On peut donc déduire de la transformée inverse l'intégrale définie suivante valable pour х> х0: 
Í di R | 1 | J, (Ах), EAR YAR J, (4%) = IE хх” 
o À "Län (J, Mäe +(X. Qo) 2 
Si par hypothèse tous les termes non dominant du polynôme de Lommel conduisent à des intégrales 
identiquement nulles. П ne reste que le terme dominant qui donne l'intégrale suivante : 


T di І J, (Ах), (Ах,)— Y, (Ах), (Ах) _ Л уу 

ТИ Б а 27 

> j dA J, (Ах), (x, )- У, (Ах), (Axo) s л xa 
UP о) susc ER 


Ce qui est le résultat recherché. Reste à démontrer que les termes non dominants de l'intégrande 
conduisent à une valeur nulle. Il se trouve que les polynómes de Lommel sont soit des polynómes 
paires ou impaires, les termes successifs sont donc séparés à chaque fois d'une puissance double. 
Or nous avons démontré par la régle de récurrence que les deux intégrales suivantes étaient 
identiquement nulle : 


Hits | 2 (J, (Ах), (o )- Y, (4x)J, (Ax,)) 0 


cu me = 
= ж s U Go) + (b (tx) 
Les deux intégrandes sont donc séparées par une puissance double de 1⁄2. Par ailleurs le terme 
immédiatement inférieur au terme dominant dans l'intégrande de départ est identiquement nul : 


Terme antérieur au terme dominant dans Rz) oc = —»oc | 7 уч (у, Hat: Y (x), (x) =0 
À d (J, Qa) + (Y, (4x) 


Ах, Š 


Оп a donc bien le résultat que toutes les autres intégrales sont nulles hormis celle du terme 
dominant. 


On arrive aux résultats suivants : 


Vv,l tq O«IeN«v 
I(x, x)= Í d J, (ax), e (Ах), (Axa) 0 


1 
9247 Ga) +E, Ax) 
ID(x, xo) | di (I.Z (Ax)Y, Hi Кш} (x, (Ax, )) =0 
| 0 A" (50) +( Ee) 
1 Dx, x.) = Í га (J, (а), Val: = (xy, (x) = 0 <= en dérivant autant de fois qu'il le faut 
QA TP (J, (Ax, )Y + (Y. (Ax) 


avr К Кт : eng 
Ate EEN T aenea" Gr 
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Si à partir de la dernière formule on passe à la limite x=x0, il vient : 
ID pls | Jua (Ax, Y, (x, )- Y. (Ax, Y, (2х) = л xy 
| dy (J, (x) + (Y, (Ax) 2“Чу -1kv -2)---(v -1) 


2 1 
De plus Jia (Ax, Y, (Ax, )- Y, (Ax, J, (Ax, ) md e Riva =] 
0 0 


1+1 


À 1 1 л? 
Si Zheng ex) 2" (v-1Xv-2)--(v-1) ° 


1 л? 
S [R 1-1,v+1 (2 Јаз = 2 (v -1v -2)- =D) 


Tirons partie des relations de récurrence sur les polynómes de Lommel suivantes : 


6) a, (1]-o «,(1]- 


1 1 2(v -1 1 
(1) Riva (4) = Sal J = (ë n ( | 
x x x x 
1 1 2(у +1 1 
(2) К, DE sl | zx (g s. | 
x x x x 
1 1 2(v -1-1 1 
(1) = Riva H + fs | = _ dÉ | 
x x x x 


Réalisons une combinaison linéaire de deux intégrales successives, permettant de reproduire la 
relation (2), il vient également un résultat notable: 


+? dâ 1 1 л 
Í SZ R o 


SAT 1) (¿QJ (Q) ` 2"2(v-1Xv-2)--(v-0) 


[s A. 1 c 


MUT 1) LI T Q  27(-)v-2)-(& -0i -1-0 
= [ж\з SE di _ 7 8A | WG) SC Gr" 


1 


EE 


v 


2 


1 c dÀ 1 1 
= uut [r "(a осу у "° 


Cette dernière relation est très pratique pour trouver des nouvelles intégrales définies : 


dA 1 1 
p dp Ua ` 


TE E 1 EE 
Nous savons que : 2 $4 (0 Y «(GG &(у—1) 
ss dA 1 " л? 
2 $4 AVEA 326-162) 
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Prenons [—2, les deux intégrales précédentes vérifient bien la relation établie: 
1 4(v AER 
AIDES 
i À À 


de A) (s, 0) +6.) o À x 
=. — 1 = 1 [es 1 
02 HAT +(X, (a) 4v-1v-232 HIT +(Y,(A)Y 
Prenons [—3, il vient une nouvelle intégrale : 


TEE 4(v -2) | &v -1)(v -2Xv -3) 
"UA À 2 


‘Фад 1 1 a. rdA(2(v-1(v-3) 1 1 Е 
(allegueren crc 
EN [2 1 _ 1 (Z 1 

, A QU) «(GG 2у-1у—3)% 2 LD Y 

РАА 1 Е л? 
7 | X QG) «(Gy = 64(v-1)'(v-2)v -3) 
Prenons [—4, il vient encore une nouvelle intégrale : 
Pedum ВЕ 12(v -2Xv -3) 16(v 1Xv - 2Xv —3)(у -4) 

4,v-4 à A EN 

үй 1 1 о [44 1 = z? (5v -11) 
la ^1 Joy «xy | X (J (a) «(y 512(у—1)(у—2) (у -3)v -4) 


Transformée Généralisée de Weber-Orr vers l'arrière de la fonction puissance x" : 


Prenons également la fonction puissance : f(x)—-x" . La condition imposée pour la convergence de 
la transformation intégrale est notamment la finitude de l'intégrale : 
+оо +оо m v5 3 
[ai^ к= [ati ? <o> Limit? 202 u <-> 
Ax; d. t—>+00 2 


Le comportement asymptotique des fonctions de Lommel de deuxième espèce à + est connu: 


sosi (ui =v? [ш 1) ch з) "о 


6 
x 


1 1 1 1 
=> > J, x > OY, x => 
œ x 2 Y (x)œ x 2 а ji d his 
v 
Ох Ох 
Rappelons les intégrales indéfinies suivantes faisant intervenir les fonctions de Lommel de 
deuxième espèce (mais également de première espèce, toutefois ces dernières sont moins commodes 


à utiliser): 


Celui des fonctions de Bessel également : y (x) 


[а z"J, (2) din en AM (20), del eg, C)] 


S. (z) Lommel deuxième espèce 
[а z"Y, (2) = zhu +v KAS 7 Y.,025,, 2] 
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Il advient donc les expressions suivantes des intégrales définies liées à la transformation de Weber- 
Orr généralisée, uniquement valable pour u«-1/2 : 


[dz ar, dm v 19008,0074 408,4, ell 
Га z^ Y, = «и +» Dy, 8, a0 el, eil 


[at oke J, D = Ах, EZ (Ax, NS aa (Ax, )- (u py 1) J. (Ax, IS E (Ax, )) 


Аху 


[ace Y, (@)= ar (2), а (Ао) (иу D Y. (ахо), Qo) 


Fa) far A 40-4 Gs far re" ul 


Xo 


= Flac) he ri jan -s0 jani) 
2x0 


Ахо 


J (Ax Ax, YX, v-I- x oS, asas) (uv - DY, (Ax,)S,,, EET 


= КАЕ = rx) a (Ax )S 4, 00.) (u EN 
= Fla x.) = Ke о E (4x,]- 


Е |= (uv -1) 5 ia JU, (Axo Y, (Ax,)- y. (Ax )J, (Ax, y 


v 


2 1 
J, (Ax, Y, (Ax, ) 2 E (Ax, J, (Axo) EGAS Riava =] 
лАх, Ах, 
Сотте 


=, (Ax, D. (Ax, )- Edi (Ax, )J (Ax, ) SSC 28 Rio = 
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L'expression à simplifier est la suivante : (u+v 218... (2) i (x)=R, (ijs ax) 
—Lv-l- x u,v-l- ,V- и+1,у— 


Utilisons les récurrences sur les polynómes de Lommel : 


€) a, (1)-o a.(1)- 
1 


x NEE 
1 1 vcl 1 
(2) hal ES x ` ) 
x x x 
2(v - 1— 1) 1 
E J+ Ran- 1— E ps | | 
x 
1 2(v 1 
()=> GE ka. E j- К, BE | 
x x 
2(v - I). 1 
К, 3,v- "E E 1—1,у BE Е К, (2) 
x 


Ainsi que les récurrences sur les do de Lommel : 
(1) 2v Dus v (x )- x(u +V- es (x)- x(u x di 1)S лы (x) 
(1) — Xv - sx (x) = х(и HAE Ета (х) =, x(u =V Г} ы (х) 


(x) " Xv — (x) ( 


+H\U-V+ DIS esa (x) 
(1) = 2(v =} uz US oc (x) = х(и rs І + 2) ema (x)- x(u ME ЕА (х) 
S (x) a x{(u +V -1 T 2)S 056 (x)- (u SVT DS aos (х)} 
u+2,v-1+1 
2(v - 1-1) 
(2) = Dee (х)= SS —(и my l + 2Xu- V+ DS, ua (x) 


e cut) 
= j e Z u+2,v-1+1 
> (u VE ) ВЕ о: (и +v-/+2) 


— (и+у—[/)$ 


u, v-l-1 
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Il vient : 


l 1 
Roos x ш+у- DS osx) Yo LES 
x 
1 \[2(v -))S a, x 
= К. 1,v +11 2 ( ) ит. -( + (а v +1)5 Da ы] l- Ra 1 d 


X 


2 l 1 1 
(u- v +1)S uv- p. = ) n. 1v+1- (2 )- Кр : Js... 
x x x 


> tes NE SEN (х) 1 
= Raiva + Rip x: Sas to) 


(u+v- 1+2) 


\ H +V — l xx 2)S EE (х) 


x^ (x)! EN 
Raiva $ оу А+ Куә n 


(u+v-1+2) 


1 e 
w 1+1)х^ ERR. Lut E J u= v +1)S uv Jal) 
2(v -1+1 1 1 
[e Rares) R 1,v+1- É Je kk: 1+ 2)$ 525 (x) 


X 


2(v - Ix (u *v -14- 2) 


1 = 1 
ЖЫ; -l+ 1)x” + XR aval E Ja- V+ E (х) 


1 
+ ХК, 3,v- ai |7 KS de І + see (x) 


2v-1+1Xu+v-1+2) 


Le calcul a l'air juste, mais en réalité il n'est pas tellement plus simple. Le seul mérite c'est de 
diminuer l'ordre du polynôme de Lommel : 


1 1 
(и kk Dass H Sav (x)- Jk (х) = 
x x 
+ 1 1 
2(v = I+ M "Rasa r х(и -V+ дв... E) Se (x) : 


1 
+x(u+v-l+ Жыны e (x) 


2(v -1+1)(u+v-1/+2) 


Vérification si l=2 : 


1 1 1 2(v —1 1 4(v —1{v – 2 
[22 R. 349-143 B -0 D = A. (t T ( ) Rs | SS ( X 1 


us, ach AUD 5 a El 


2 2(v –1)х“ + (и -V + 2)S ne (x) 
(uv) 


= 2x(v -1u +v -2)s,, ele be - (v -1)Xv —2))5„.„.„(х)=х 
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En inversant la transformée, il vient l'intégrale définie suivante : 


Y 


v+l,v (2, x, x0) = J, (Ax)Y, (xs) Yu (Ах), An, 


Vxx>x, f(x)=x"= [44A Е(А,х,) 


o (J, (Ax, IT + (Y, x) 


1 1 
2** dÂ (о +у -!) dE jte, [Es a fe (A,x,xo) 
= Vx x>x, xs Í 0 С 
л 


E DEG Can 
1 1 
" E a +v -!) EE ү (2,x,x,) m 
Vx x> | ) L 
> vers [e GER d 


Transformation généralisée de Weber-Neumann vers l'arrière 


Ce qui suit est tout à fait expérimental, et sans aucune preuve formelle, mais imaginons que 
l'inversion de la transformation construite comme telle, soit possible : 

F(2,x,)= [а x DE as fo) 
Y (Aux xs) Ах), (Ах) Y, (x), (Ах) x 


/(х)= | ала = s. 


Lef + (F, (o) 
Plus précisément regardons les transformations avec le paramètre u=v-l : 


En 


PU. x,) = | dx LE (A,x,xo )f (x) 
j (4, x, Xo ) — J, 4 (Ax)Y, (ax) m (Ах), k (Ax, ) Ze Ü 
f(x)= [алл +, (хх) F(,x,) 
= E 
о (х) +E äs lf 
Calculons Іа transformée généralisée de Weber-Orr pour u=v-l pour une fonction puissance x" : 
fax E E (х) = а) fax X. ДЕ (х) = a Y sax) 


1 
LS VORE Da oc x. => Pour v >1+ i Lim x" J, у(х) =0 
ae e (Аа) 


) [ax жу (x)= — 


+œ 1+1-у 


=> [ас x" Q, (Ах), (ахо) Y, (ax); (a =, x) os) X ay QV, (s) 


Xo 


[ax xu (Ax)= 


Jeux 
À 
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Le produit croisé obtenu fait intervenir les polynómes de Lommel qui sont en relation de récurrence 
comme suit : 


x 
1 1 2(v -1 1 
(1) Riva (3) + md | = i dl | 
x x X X 
1 


DUC x x x 
1 1 1 2 1 
=> Riava H = Riava H + f odd H Fem Riva H 
x x x x x 
Le E L 


Le produit croisée se présente maintenant comme suit : 
аа) zf ah, Ais LL = 3 v, (22) y) 
Ох Ax Or Ах 


= J, '(Ax,)= (о) 2, (x) Y, '(Ax,)= (x) Y. (xs) 


| iste Ex Ga) 


Y, at (ln. (x) . ‚(Ах ) 


(oat, ym (Ax, )- Y, (4) us W (Ax, ) 
> Јо )Y,'(Ax,)— У, q Ax, (Аха) = v 


= Uns) а). 


=? E AN у, (a, )- Y, (y, M, (Ах, ) = 


0 
2 1 
WEEN 
0 0 


2 1 
Comme 4J,,, (Ax, RK: (Ax, ) =Y; 44 (Ax, J,a (Ax, ) GE Ryu [=] 


2 1 
Jta JY, (2х, )- Y, ua) (2, )J, (Axo) vix 8. (a 
0 0 


2 1 1 
д) aas M Qe Gnd Fd) 
л. б б Ах 
е ; 2 1 v 1 
AJ (Ax, Y, (Ax, )- Y, (Ax, JA (Ax, ) E e HE DE WEEN Ra Z) 
0 Xo 0 
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Cela donne l'intégrale définie suivante : 
+00 oe, ; 2 lv 1 1 
= |a (оа) Axa )- а) (а) CH | 2 R =] dE 


Оп peut donc déduire de la transformée inverse, si tant est que l'inversion soit prouvée, l'intégrale 
définie suivante valable pour х> х0: 


вание 


А Ах, Ax (7, (Arx Qc (ax, y 27 


0 


Xo 


Conjecture générale qui reste à démontrer 


"t di J, (Ax), ax, )- Y, (А), '(Ax,) = 
PEN [uu «GUY о 


v 
De méme en dérivant par rapport à x, il vient : 


T di J,a (АӘ)Ү,'(Ах„)—Ү,_,'(Ах),'(Ах„) EN 
ем [ari Gaia) OO 


Cela а ипе conséquence importante, comme les termes dominants des deux polynómes de Lommel 


E | 1 лез е, 


1у—1 1.4 
Ах, A Xo 


Ra ) _ 27 -1v -2)--(v - 141) 


Id 
À Xo 


sont respectivement : 


Tee 


Alors : (44 l J, (Ax), (Ax, )— Y, (Ах), (Ао) _ 
kE dë | (J, Qo) + (Y, '(Axo)) : 


0 


L'intégrale indéfinie prend alors la forme : 
Td ï 2 (йж) Un) = Cy A) E v-l Iv 
е (О) + „(ыу AT" 
0 0 v 0 0 
On constate pour les mêmes raisons que seul le terme dominant du polynôme de Lommel reste, il 
vient alors une intégrale définie de forme assez simple : 


[ыле 2-62, В | (Ах), Qo) - Y, (Ах), (Axo) _ T opel 
| ж al baffen 


v 


0 


Í d^ J, (Ax nA JAY v-l (x), 0x) _ л ун Г-У 


> ш - x 
IU GR) 2v 0-2) LD" 
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Maintenant le passage formel à x=x0 peut être réalisé sans dommage, et l'on obtient : 


T 


+ ` 2" v(v — D -2)--(v =) 
v 1 
R — |-——R SC 
+œ dA | = Ах, sal ) 
=> Í 
+ 


Ах, - л? SE 
0 À (J, Qo )Y Q7, (2) 2 (e -1kv -2) 6-1)? 
1 v 1 
fi dA t TEN) л? 
4" CAN CS) «GG 


` 2? y (v -1Xv -2)-—-(v -1) 


1 


A (у, (A) «Gy 8-1) 
2v(v -1) 
| ад | 2 -1) _ л? 
A LA «(OG = 16v(v-1)(v-2) 
+? ад 1 л? л? 
Pour 1=2 : orle DJ A (4, (Q)Y +(У, ANY 7 16v(v 10у -2) 8и (у 1) 
ал 1 "mE I sl 
эмел нуы E22 H 
_ [ d^ 1 л°(Зу —4) 


SX UI «QUO ` Sävche hu 3 


dérivée première. 


Par contre pour des valeurs paires ou impaires trop importantes, il y a un problème de convergence 
de l'intégrale numérique probablement lié à l'implémentation des fonctions de Bessel et à leur 


Résolution des problëmes aux limites de Laplace en cylindriques 3D par les transformations de Hankel et les 
transformations de Kontorovich-Lebedev- p 212 


Гаа) = ахо (0, (Ако) а (о) (ижи 0) 7, än, än, 
Ах 
Í dt p Y, (0) = AXo (as (Ax, E (Ax, ) is (u "V< 1) Y, (Ax, ДУ ўе (Ax, )) 
Ах 

Js Ju (xo) =z län! Y, '(Ах,) = e (Ах) 


Xo 


SESCH A Apel 2 


= Fx) zs GU [a J, ts face v, o 
Ae; Ts аа s) Qt D (S, n) 
EA J, (Axo NU, (2х), (Ax) (и +v =D) XQ), (Ax Si 
" T Ax, k Ee Y, (Ax, Je... (Ax, ER (Ax,)- (u +v -I)J, 4 (Ax, jeux J- 
x — | (4x,)- red (Ax, Jic. (Ax, Sas (Ax, )- (u EVE 1) Y, (Ax, EE } 
K E < v-I- x SÉ, (x )}- | 
F(a j= Xo - (u +v -1)S py- Ax KI, (Ах Y, ax »)- Yo (Axo J,- (Ах 5) 
"le [Suas (о aa JE а) а) (A 
а sU E Eech Y, As, H „(Ах " 
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Ce qui donne le résultat final : 


A (2х, Y, (o) =Y, _; (Ax )J, (4x0) =- = + Кулы E 


2 1 
> Л, (Ax, x, (4x,)- Y ua (Ax, V, (Ах,) LI R, d | 
Сотте 


1 1 
BS gis (Ax, HEN Cé (и du 1) Syst (Ax, ЖИР 
2 x Ах, Ах, 
= F(A,x,)= 0 


lli u+l 
лАху À 


v 


1 1 
y" Ах, ү (Ax, Ж = (u dy 1) S угу (Ax, e Ei 


1 1 
A (u PKS 1) Suv- (Ax, ) = = SS Riava GJ Em 
=> F(A,x,)=—— 


nA" 1 v 1 
B S ust (Ax, le E = Ax, PS 


Une fois cette fonction inversée, on obtient l'intégrale indéfinie suivante valable pour x> x0 : 


1 
К, et ) 
` Ах 
(u+v-1)S „1 (Ax, | 
+00 H ' ' 
Í di E Ж | 1 | J, (Ах), (axo) -Y (Ах), (4x) л 
Ди"! —1,v+1- я 


0 
1 2 1 
a Dv (Ax, EN Ах, sl ) 
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Une formule de transformation de la fonction Puissance utilisant un noyau mixte vers l'arrière 


de type Weber-Neumann 


Soit le noyau mixte suivant : (4, x, x,) -J, (ах), (Ax,)- Y, "Gel. (Ax,)= 


1 OP EE (2, x, Xo ) 


À Ох 
Et la transformation de Weber-Dirichlet vers l'arrière : 
Vuoi), p Du m (a), (ax) 
W, Lei РО.) = [ах xW, is Qs) PG 
о) Fa. TT 


CEE) 


Par la transformation de Weber-Dirichlet vers l'arrière la fonction puissance devient : 


2 1 1 
F(2,x,) = His b" j- aA? L du 1) Rivaz Les) Ae ee (Ax, ) x 


Ах, ÀX, 


Prenons un noyau mixte et regardons comment toute fonction est transformée : 


V, (sx) eJ Qo Qs) - Y, A, (Ax) - > pop ааа) 
W, iy {f(x)} = F(A,x,)= Ía х\Ф, |, (А,х,х, )/(x) 
-zje OP, pen) ГА! )- ro. ‚(А„х,х,)] EL o, (A.x,x MfG) rte, 


T VE (xx) 
K Lr (oe lie? T| as O 


(Р(х) = xo f (Xo)P, (A, xo, Xo) 1 W bon )+ / Ae 


v-l,v А 3 v-Lv x 


. 2, 1 
Comme : w. (4, x, x) (Axe) s) Y. (Axe M, (Axa) = ed 


Avec la — puissance, cela donne : 


Е 1 1 
W, x '| a qurl Kal. H+V— [= 1)R 1— DE (4x,)- slk. (Ax, ) 
0 0 


2 AIX 1 2(u +1 1 1 
Hu be j= = Sg = " Ho +v-l- 1) E (Ax, )= HE 
0 


ÀX, 
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Transformation généralisée de Weber-Orr vers l'avant ауес ип ordre entier selon A.V.Gorshkov, 
2019, « Associated Weber-Orr transform » 


Un article tout récent donne une formule d'inversion de la transformée de Weber-Orr vers l'avant, 
mais uniquement dans le cas ou la transformation intégrale est construite avec un noyau dont les 
ordres des fonctions de Bessel sont entier. H s'agit de l'article de 2019 de A.V.Gorshkov 
« Associated Weber-Orr transform » sur un décalage 1 vers l'avant. Je donne ici le résultat de 
l'article : 


v-k-l keN> Paian) = 7, (Ах), (Ах) - Y, (x, (Ах) 


F(A,x,)- je x5 (A 2 Ce) 


xo 


E 2(k-1) +œ +œ 
eebe => ы 
no x e 0 С (А) * (Y, (Axo)) 
Гадка И) E _ 
sxo) 2 = К=0 ou К=1 
0 (J, (Ax, )) +(Y, (Ax, )) 
En décalant d'une unité k->k+ 1, il vient : 
Pona (4,x,x9)= Jp (Ах), (Ах) У, (Ах), (Ах) 
F(2,x,)= нш 
— M" Xa (A,x,xo) 
ds dAAF (А, 2 k21 
ï [a L) + ar N 


fe)- Poux (4, x,x,) 


(J, (Ax, Ur (F, (25,3 


jaa A FG x) k=-1 ou k=0 
0 


La formule est simple et a le mérite de marcher parfaitement avec la fonction puissance CT" 


V, au(x) = J, (Ах), 0o, ) - Y, (Аж), (Ах) 
Sg? +œ _ 2 -1-k 1 
Pour k21 (х) = х" => F(A,x,)- ` ag NEN REO EP =0 
EE 
f(x) x =x Jx € L (x, oo) «= fa x =——° = 
1 1 
Sarp, VA È 
2 2 
_ 2kx E (x) _ 2k x мазе ds 2k x,” bu 2k x,” 1 
Ai [art] ei [ ZA — M | ve = S зк = ыт €4 d 
0 


Pour k=0 ИЕ (4, х,х,)= (x) 0, )- (x) (Ах) [as (x)=-J,(x) 


1 So А 
f(x)Vx -——€l, (x, +00) l'inversion n'est pas définie pour cette fonction 


Ух 
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Voyons maintenant si lon peut générer des intégrales définies à partir de cette transformation 
intégrale. Prenons une fonction puissance quelconque x" : 


Í dt t" J, (e) = ÀX, (J, (Ax, Biar (Ax, )- (u +k+ 1) Jin (Ax, Je, (Ax, )) 


Лх 


Í dt pe Ү (c) = Ах, (v, (Ax, JS ga (Ax, )- (u + К + 1) Ya (Ax, Js (Ax, )) 


2% 


пал) rn ae к) f de en n| 


Xo Xo 


1 "t + I: 1+ 
> еб) sar t" “Лы (0)-J,(âx,) | dt m D 


Лх Axo 


Ах, Y, (Ax, Yu. (Ax, JS es (Ax, )- (u + k + 1) Jua (Ax, Je, (Ax, )}- 
=> F(A, x,) = Pr 

À e. J, (Ax, yl, (Ax, IS еы (Ax,)- (u +k+ 1) Ya (Ax, 5,4 (Ax, ) 

xou +k +1)S „(Ах 
> pla, x ) З а) аур s) 1 Alt 

2 1 2 
Comme delt Y, (4x,)- Yai (Ax, V, (xs) = ÀX, P J = лАх, 
2 1 

— F(A,x,)- (u+k+ Je, (Ax) 


r qut? 
Dans ces conditions l'inversion de la transformation donne l'intégrale définie suivante : 


bom (2, X, Xo ) = J a (Ax)Y, (Ax,)- Y, (Ах), (Ax,) 
2(и + k + Us. (Axo) 


F(2,x,)= SE 
„u _ 2k Ce (^ P 2(u - k +1) | As Dua (Axo) м, (А, хха) 
M 1+1 2 2 
x A л d 2^ (7 A (Ax, )) + (У, (2х, )) 
Е 58 2k x di 2(u +k+ 1) Та S (axo) Poux (xx) en 
0 


AT GA) + (V, (Ах) 


PA ç (Ax ) y (2, x, x ) л 2k x u+k+1 
dà uk 0 k+l,k 0 = Hu == _ => >. k>1 
> |a m ү ax.) +G) EEN 


Passons à la limite x=x0, on retrouve une intégrale définie déjà calculée en usant de 
transformation de Weber-Orr vers l'arrière. 


EN 
TÀX, 


Poux (A, Xo» Xo) - Ja (Ax, IS (4x,)- Y (Ax, JA (Ах,) 


) л u-k+1+2k 


o A7 ху (J, (Ax) + (v, GA 2(и+к+1) (и-к+1) 
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Transformation généralisée de Weber-Orr vers l'avant selon C.Nasim 1989, < Associated Weber 
integral transforms of arbitrary orders » 


Élément de calcul pour preuve de la justesse du calcul réalisé par C.Nasim en 1989 dans son 


article : « Associated Weber integral transforms of arbitrary orders » 


En section 4 de cet article l'auteur entend démontrer la formule d'inversion pour l'opérateur 
fonctionnel W et son noyau de Weber-Orr vers l'avant, qui selon ses notations s'écrit : 


R, (€, p,a)= (ёр), (6a)-Y„(E0)J (6a) — Won, if) [de e Ra (E. p.a)f(p) 


bi 


Pour cela notamment, il introduit les opérateurs fonctionnel d'Erdélyi-Kober et leur version 
modifiée : 


-2(a+n) 
Erdélyi— Kober 1„\/;х}= FE ја e (x =Y ri) 
~ 2(a+n) a — 
Erdélyi - Kober modifié Lal ;x;al= 2x- Ta) ја gx zi ' f(e) а> 0 (voir formule 3.5) 
ag d г! 
T.M pid = [arl di (p 0) 


C.Nasim introduit alors le noyau modifié par l'opérateur d'Erdélyi-Kober modifié : 
" p `> Í Ea *p -(v+k) £ JE i 

R, (6. p.a)=| | 1, AS. ;a dt(p? -1? ] t" R, Lët, 
„lé. p.a)=(2) T, le (Est, a) psa) PP Sal Esa) 


Le noyau de la transformation de Weber-Orr posséde une propriété intéressante : 


k 
E 2) lx R, (E,x,a)f= EX ERU a (E,x,a) voir Lemna 3.5 


Pour cela différencions le noyau R, il vient : 


а (а, osa) =E, ENa- EM, (а) = (P689 


+E (х) (6a) у, (Ex) Ea) = "(v Ex), (62) 7, (Y Ea) ER, i (Esa) 


uY, (£x) J (£a IS HJ, — (ga + 


Е Ё £). (&,х,а)}= SR. (Ewa) ИК, LE. х,а) 


ах 
SE ze “R „(&,х,а)}=&х R, ,(E,x,a) 


= (х&-+и}{н,„(х,а}}=&х“ Ra enn) m Z (R. Go] E Ras Goa) 
X 
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En appliquant l'opérateur différentiel vers l'arrière pour u=v on arrive à la formule du Lemna 3.5: 
pus Ë bu ue vi 

u-v— (2) R, (E,x,a)}= Ex MC iy (E, х,а) => 22) {х к,„(&.х,а))= СРИ R; wu) 
х ах х ах 


On peut également construire un opérateur vers l'arrière en utilisant les formules de différentiation 
des fonctions de Bessel J et Y, à savoir : 


LED d[1,69- 7,69 o Le Gon) ик, Goa) Sep (ха) 
d -ua[ „4 -u 
RE M, G oss (i-na оа), ens) 


1d E B 
(iz "R. (E x,a)] - Ex G UR a (£,x,a) 


Cette formule de différentiation va être utilisée lors d'une intégration par partie sur le résulta 
essentiel de la partie 4, à savoir le calcul expression 4.5 du noyau modifié par l'opérateur 
d'Erdélyi-Kober modifié : 


xa y k _-(v+k) " 
STEE EE v R. x,a); p:a а}= mL - x?) x""R (E, xa) 


z“ 


Pour cela concentrons nous sur l'intégration partie de l'intégrale, en première étape : 


1 d ç dus | 
Сотте É zen. (5а) x" R. osa) э Ra ,(E.x,a)i=5x*"R„,(5,x,a) 


и=у R. (Ex) Ex Rosa) 
X 


[etre nee) Lact ef" ea. (a) 
Alle enema jac eoe) 


k-1 v+ 
La e eaa) UV en. eler ef 
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Réitérons pour une deuxième étape cette intégration partie cette fois sur le second membre de 
l'expression droite : 


ш, DAP Ras osa) o E Ras (nsa) 
X 


Al YE 
p lei e : qd 
=$ fax (p? -x f vin (Ex, a)= 


E » jet dx E: yo S AR (Ce) 
k CEGA Nat ln 
‚= к, (баа) E Ap een ach 
I B = = ахх", (E, хао? х?) 


Оп voit parfaitement se dessiner Іа récurrence jusqu'à l'itération k-1 : 


p 


Jax(o? -x Y 8, (Exa) = 


a 


Jl, с, Le 


Eas E? 
E k-2 2 2y-6).— vx k-1 
k -1)(p?- k 
E l DEN ш Ru @ (£.a,a)+ (k TE D m [aen E eR T lË, x, a) 


= so 72 v+2 
"ee E eech: 2(k Uer? SC SS E I es (E,a,a)-...- 


ENEE ҮТЕ ачка) 


р 
v4 1 v+ v v+ 
Or faxx aas esa) [к "Ray (E,x,a)[ => "К, (£,p,a)- a ER sky (£,a,a)} 


2 2Y- v+ 2 2Y-2 v+2 
LA 2(К—1 = 
a a" n Geet A Kai: R, (6,a,a)-...- 


š 
p iva 2*- 2 k - 1) = vc-k-l Des k —1) a 
= [а (о? -х?} E. 'R „(&,х,а)=+— l e 1 ы d —— € "Ryuu, (6,a,a)+ 
k-1 
+ LS (E, p.a) 


p " ach 9m (52 _ q2 Y^" qvem (fe 0 2" k=-1 v 
> [dx (p? =x?) "x 'R, (6,x,a)=-Y_) (o ar] "a ( Lo ëch e ep “Res (E, p.a) 
a mer à 
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D'où l'expression obtenue par C.Nasim qu'il injecte dans l'expression 4.4 qui réalise la formule 
d'inversion de la transformée de Weber-Orr vers l'avant avec un indice k entier : 


R, (&,p.a)=2 MED falo 2x) tr, (Ex a) 


E pu ES т=К Dm (p? < а? -т ау" (k – OP 


т=1 E" (k = m) 


=k m-k vim 
zu (vk) N G J п, ? _ aY" R, „ (E,a,a)+ te, 
m= d 


ZE „ы 
——————— EEN 


Nos notations : 
On revient à nos notations par la substitution : 


a=>x, ё—>А px R, (Gala 


uwv 


(4,x,x,) 


W. {f}= [do p Ra E p. a)f(o) 9 Fla, x)= Í ах, (Ах) (б) 


Xo 


Résultat : donc d'après l'article de C.Nasim de 1989, < Associated Weber integral transforms of 
arbitrary orders », avec un noyau de Weber-Orr généralisé de la forme u=v+l, l'inversion de la 
transformation de Weber-Orr généralisé vers l'avant (formule 4.6 après substitution de mes 
notations) s'écrit : 


et x, х)= E (Ax, )- Y, (Ax 4, (Ax, ) 


F(,x)= je XV, (Ах) (x) 


Xo 


Lye = Date TX 2 d 1+/-т ы KE Xon) e 
| | (I—m)! [ал RU de (Ax) + (Ax, T 


Nn Mos LA x, Xo) 
Meter 


т=1 


f(x)- 


L'expression peut encore se transformer un peu : 


v+m 


Erne: Ua) ME EEN 


лАх, Ах, 
хе m-l 2m x -1 ` V+m- (е Sak UP +00 1 (4 m ) 
= $y 0 0 — 4 m—1,v+1 2 770 2 + 
a 1-т)! ` N Ax J (J, (Ax, Y (Ах 
r=) zl ©” CEET 


+ [44A Е(А,х,) 
0 


(J, Qo) +(X, (х)? 
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Le cas particulier de l=] est le plus simple, il vient : 


wa oes) as (s) lan Qn) A ll Z 
o... охь) Ach Qs) 1 Delt (am) 
F(2,x,)= Jaca Goss) 
^ E ` F(A, x) m SP as Xp) x 
fee I cras «a ax 1^ Са e 
Continuons le calcul : 
_ Му (А, y. x,) ei ET x, Xo) Б 
r= а оуу О? ere re A 


Intervertion des M MN 
Pou (2, YX Xo) ni P out (2. xX, Xo) 


ас чы” anren eer 
Dosrt-omculwler je SP ИЕ Y. Xo) “| J (Ау), (Ax, )- Y, (Ay), (Ax, d 
оу на Т Оа а) 


Si l'on fait une trés grosse supposition, alors on obtient la formule de A.V.Gorshkov (voir article 


de 2019 < Associated Weber-Orr Transform > généralisée à toute valeur de v 70 : 


+00 2v +œ 
[а J (Ay, (Ax, )- Y. (Ay y (ax) = ZVX) y = v _ 2v Xo ja SCH w Ë G, w) 


ag E (7, (ахо) Qr (ах) 


Démonstration partielle du résultat : fa a (А), EE Y, (y M, Un, DE Tee 
0 (J, (Ax, )) + (Y, (4x0) 


J'ai cherché un peu partout dans la littérature oü je pourrais « prouver » cette assertion, et j'ai 
trouvé dans un article de 2009 de Z.Xinhong, T.Dengke « А generalized Weber transform and its 
inverse formula »ипе formule utilisant un parcours d'intégration sur le plan complexe qui me parait 
commode (voir formule 2.1 de l'article), soit : 


N, (А.х, x) e J, (AxB Y, (Ах„)+ AY, (ахь) УА) J, (Axo) + a, Län H 
D (2, x)= {B J (Axo)+ a4, (ax € 4B Y (Ax,)+ aA Y, (Axe) 

ү” Nos). Z x ох” 

3 Å D, (Ax) = 2x,B-av 
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En posant a=0 et B=1 : 
TA J, (Ax)Y, (ахо) Y (Ax (Аль) ms 


14 


„4А р(х) ` 27 


v 


En dérivant formellement cette expression par x, on retrouve l'intégrale que l'on cherche, mais c'est 
la ой potentiellement il y a une « faiblesse » dans la démonstration, car la dérivation à l'intérieure 
de l'intégrale est-elle licite ? quoi qu'il en soit on obtient : 


eg All An Gs) (an) (a) (a) 


n ER e о" 
v 12 J (Ax)Y, (4x, )- Y, ( d Ju (Ax)Y, (4x, )- Y, (Ax), (Axo) ET T V val 
TE s UG) доо? 


= yx, vx je б M | 0) Л у =v- 

2 V, (os) + 0)? 

> [aa T Rr (Ac) T. (Qr 0) 
(J, (Ах) +4 Ах) 


= улху X 


се qui est la très grosse supposition réalisée, mais qui me paraît cohérente. 


La formule de généralisée de A.V.Gorshkov marche parfaitement avec la fonction puissance x^" , 
selon les mêmes calculs réalisé avec v entier, mais pourvu que v21 : 


+оо -1-v 
Pour v21 f(x)=x"" > F(A,x,)- [а x” EE s — -0 
л 


Xo 


à +00 4; +) б io) zd 2k xy" +оо 1 
отте [а x EE PE fx =x ""/xeL (x, +œ) => f(x ei [af] SS c.q.f d 
Xo 75 


Calculons la transformée généralisée de Weber-Orr pour u=v+l pour une fonction puissance x” : 


1 
x" J. (x) œ x ? — Pour v > -1+7 Lim x" J, (x)= 0 


1-1-v 1-1-v 
fax х -I- "J. (Ax)= x 7а) [а xY, (x)= X Y, (Ax) 


1-1-v 


> [a x" (Js (A), ако) YA), (Ax) = "Mali Qs) ua M (Axo)) 
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Le produit croisé obtenu fait intervenir les polynómes de Lommel selon la formule suivante : 


2 1 
JG). (Ахь)—Ү, (Ax, EEN 
лАх, Ах, 
Сотте 


=> J, (Ах) (Axo)- EEN 


Ах 


3. Ax E зь), s, баа) 22-8151) 


Оп peut donc écrire l'intégrale définie suivante, décliné selon différente valeur du décalage avant 


l: 


Ө Bee 1 

UE EEN 

Pour 1=1 [ax x" (Z... (Ax)Y, (Ax,)- Y, (Ах), (Ax, )) =0< ES -0 
Xo 


Xo 


+оо SH 
Pow 1-2 [arx (J. (AJE, (x) X, (axe 22 =) 
0 


Xo 


+œ 3-v 
рош 123. faci Qs) Qn) as (ape 38 D ед (A) 2620 
T Xo 


Xo 


De maniére à tester la formule de C.Nasim ou de A.V.Gorshkov généralisée , réalisons une 
transformation vers l'avant d'une fonction puissance quelconque : 


You (2, x, x)= J, (Ах), (Ах) Y, (x, (Ах) 
DE wee 
_ 2V x,” Nn You „(2,х,х,) 
ra =] O dii : F(2,x) 
Ja Гал eg ma p ra 
Ond аргыый ы ИЕ » э) Ts, (ax) 
VIA 

2v X, “Ж Les 2(u +V + Sal Y 1 (2,x,x,) 

= Hu 0 d u-v d S, v+l,v 

dicm л | та Ius BO 


BE dy x," yurt 90 2(u tV +1) ад 5 (Axo) Pay xx) 
H =) l Xo Л A (7, (Ax, ) + (Ү, (Ax, ) 

= [aa Suy (Ax, ) É y (2,x,x,) _ л H 2y yd 

AT (J (Ax) (У, Ax) ` 2(u+v+1) x" u-v+1 


v 
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Si maintenant оп passe à la limite x->x0, il vient : 


2 1 2 
p e (2, x, Хо ) = Jd (Ax, x, (4x,)- Ya (Ax, V, (Ах, ) = EN Е 


TAX, Ax, ) = TÀX, 
pM w (Ax, ) 1 d л? 
ce GR dv d 
= Í ал Bus (A) _ л? 
аа ау «GP av 


J'ai testé numériquement cette formule, elle semble être juste ! 
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Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation 
de Laplace : problème de Dirichlet axial sur le demi-espace z>0 


Si nous revenons au développement intégral de la solution du problème aux limites suivant : 
AT(r,z)=0 re[0,œ0) T(+œ,z) fini 
zez] 70,2) „=Л() Te = AC) z e [0,40] T(r,2)..,= f(r) 
C(A)Sinh(Az)+ D(A)Sinh(A(z, —z)) |j. s | dÀ AJ (Ar)e SCH) 
Sinh(2z,) d 


II paraït alors assez naturel d'identifier les fonctions paramétriques d'après les transformées de 
Hankel respectives des fonctions limites, et d'aprës le théorëme de Hankel, pourvu que les 
fonctions limites respectent les conditions requises, on écrit : 


c(a) = je rJ (Ar) f.(r) D(A)- je rJ, (Ar) f (r) 


T(r,z) = je AJ (Ar) 


Imaginons que le problème aux limites soit sur le domaine z également non bornée, le demi espace 
des valeurs positives, il vient une solution de la forme 
AT(r,z)-0 re[0,œ) ze [0,+œ] T(+%,z),T(r,+x) fini T(r,z) 


z=0 e f(r) 
T(r,z)- [ ал AJyAr)C(A)e^ c(A)- [ dr rJ (Ar) f (r) 

0 0 
Si la fonction limite est de la forme : 


а. г< 


f(r)=0 к>, = C(2)=T, [dr rJ (Ar) < [dz z/G) =zJ,(z) 


Aro 


c(A)- 2 [а 22,02) = е Any) = T(r,z) = То [dA J (Ar)J An)” 
0 0 


Le profil de la solution sur l'axe z est alors le suivant : 


r =0=> T(0,2) = Ts, (dA J (0)J,(Ar,)e * < J (0) =1— £ == =n [dA J (Are * 
0 0 


0 


On peut calculer explicitement cette intégrale d'après l'intégrale définie donnée dans les tables de 
A.P.Prudnikov, Yu.A.Brychkov et O.I.Marichev- < Integrals and Series, Volume 2 - Special Functions > 
page 183, section 2.12.8, formule n°3 : 


[dre "J, (cx) = ЕЕ + Jp +c? | Re(p)» IIm(c) Re(v)» -1 
d Jp +c? 


Par identification des termes, il vient : 
+00 zl 
p=z v=l с=т, х=А= [ал ane” = 2 [z+ 2 +n 
0 Jr tn 


T(0,z) ` М n 


2 2 Е 2 2 
T, NEM DEG ) NZ t zd z^ +, ) 
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C'est exactement le profil trouvé dans le cas limite d'un cóne d'angle ouvert droit en coordonnées 
sphérique, avec un disque de rayon fini porté à une valeur constante, l'extérieur du disque étant 
placé à la valeur nulle, il suffit de remplacer par l'expression : 


Т(0,2) | x _ Vl+x? —A _ 1 


T, Ja x? Ai x? Joe re + x) 
z Т(0,2) 1 n 
x= ai T. = 2 2 = 2 2 2 2 
r J 
0 0 paliz (z +H +227 73 
Wé H Wé 


Par contre on пе peut retrouver la solution donnée dans le livre de J.D.Jackson, Classical 
Electrodynamics, page 92, formule n*3.178, pour un disque placée à l'origine sur le plan (x,y) 


TRU) M as 2r y | >r=0> NU. e Е | 
л 


T, г? +(r=nY +2 urn T, л Jz + у 


Ce qui est tout à fais normal саг le problëme аих limites du disque de J.D.Jackson est 
fondamentalement différent dans la nature des conditions аих limites. 


Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation 
de Laplace : problème de Robin axial sur le demi-espace z>0 


Cette fois sur le disque de rayon r, la condition aux limites est inhomogène de Robin : 
АТ(ғ,2)=0 re[0,») zel0,+0| T(4,z)T(r,c) fini 
g r,Z 

= = f(r) 


z=0 


C.L. + BT(r,z) 


T(r,z) = | dA АЈ (АР)С(А)е ^ = a Í dA PJ (Ar)C(2) + B Í dA АЈ (Ar)C(2) = f(r) 


P T(r,z)= Í dA АЈ (Ar)C(A)e * 
> [dA AJ,(Ar)(B +a2)C()= f(r)= 


c(A)- EH d dr rJ (Ar)f(r) 


"ES 


Si la fonction limite est de la forme : 


f(r)- T, Ped " T. "n i 
| ped MCN rJ (Ar) = |а zJ (z) =z) (z) 


Алу 


T, 


J (Ar)J (Ar,)e ^ 
c(A)- pm y [ az zJ,(2) = E ER ал”? T(r,z) = Ty, j dÀ EE 
Cette fois le profil sur l'axe z change de forme : 
-Åz +œ -Az 
Ja (0) J (Але anias О Т(0,2) = [da J (Ar,)e 


r=0=T(0,z) = Tyr, | dÀ 
| ад + B Т, A ад + В 
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Je n'ai pas trouvé d'intégrale définie dans le cas des deux paramètres non nuls, par contre lorsque 
le problème mixte est de Neumann, il existe une intégrale définie, donnée également dans les 
tables de A.P.Prudnikov, Yu.A.Brychkov et O.l.Marichev- < Integrals and Series, Volume 2 - Special 
Functions > page 183, section 2.12.8, formule n°4 


fa e "J (e= I = © p+ Jr +c 22 Re(p )> (с ) Re(v)» -1 
Ps 22 T E 
| б | Jere 


us La +c? + +c? 2 Е ce” CET. g " dog 
Ji elle frire] рр рер) dei rei 


Par identification des termes à la valeur a=0, il vient : 


eg e 


+œ —Àz 
p=z у=1 с=т, х= А = [dr tre = h 
0 À z+ + r 
Т(0,2) nf, Л(Ањ)е 1 ? Ge 1 
,z) h (An )e ra 
= 2 | аА - 245-12 T(0z) = ——— 
T, | 4 dX z; + z + r Sieg z+ l+ z2 


Application des transformées de Hankel/Weber à la résolution des problèmes aux limites de 
l'équation de Laplace : problème de Dirichlet axial sur une bande centrée de hauteur 2z, 


Si nous revenons au développement intégral de la solution du problème aux limites A ; 
AT(r,z)=0 rel[r,,œ) ze [= zo, Zo] Т(+оо,2) fini C.L.Homogéne  T(r,z) 


DE „= fir) MELIA 


(T 8 „ CO) Sinh(a(z + z ))+ D(A)Sinh(À (z, - z)) 
T(r,z)- | 4А du (Ar) Y n) - Y (Ar)J (An )) ОПОТ 


r= 


Dans l'application des conditions aux limites et par identification des termes des transformations 
de Weber : 


f(r) = | dA ^W, (Ar)F(A) 


DO AX (Ar) 


dt t V (At) f (t 
(J Ars) +(X Ar.) 2 Л 


ja tP (At) f (t) 


FG) OY еу 


Yo Ar) = Jr) n) -Y ArI (An) 
J (âr) fonction de Bessel de première espèce 


Y (Ar) fonction de Bessel de deuxième espèce 
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II vient : 
AT(r,z)=0 rel[r, œ) z e[l- zz] Т(+оо,2) fini 


T(r,z) = fi(r) T(r,z) = fr) Wo (âr)= J (Ar, (Ar, )- Y (Ar )J (ars) 


Z——Z9 2=20 


[ar tP (ae) f (e) 


= [ dÀ AW (Ar)D(A) 9 D(A)= (J ER )Y +(004%)) 


[at (PG) fn) 
(A) +A) | 
Comme premier cas particulier, supposons que les deux fonctions limites sont identiques : 
[ac e (as) f(e) 


DC + (Yo (Ax,))Y 
Sinh(2(z, + z)) + Sinh(A(z, — z)) 


E Í dÀ AW (Ar)C(A)> C(A)= 


w. (ar)= Jar) (ar )- Y (Ar)J (Aro) c(A)- 


T(r,z)- Í dA AC(2)F, (ar) 


Sinh(2Az,) 
4. Cosh(az) 
e | d AIET ES 12) 


Encore plus particulièrement supposons que les fonctions sont en plus constantes, il Е S 
AT(r,z)=0 re[n,) zel-z,,z,] T(+,z) fini C.L.Homogëne T(r,z) 


T(r,z) 


= 


2=—20 


c(a)= Y (an) аел) 2,0) [d (Y (ar) 


To To 


Агу 


J di t J (=> Lo | ju (ae) = => er. Ob: Lim t J (t) = Lim t Y (t)=0 


c(A)- =" (x (ar) (ат) A (Ar) = (ат) (ат) ат) (an) 


= T(r, z) = Ts [дА Län, EA) Yo (ar, (n Sr) Cosh(à z) 
; (4,0) +, Ur) Cosh(2 zo) 
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De plus il existe une relation liant les fonctions de première et deuxième espèce, à savoir : 
v (r) = Jar) (ar )- Y (Ar)J (Ars) > Pr) = 0 


Pol) = CAC) (аљ) – Y Gr) = Ала) (Ar) = Yn) (an) 
MEM 
Identité du Wronskien W. AJ, (x), Y (x)] = — < 

"o ye) TEQ)-xu6) 


=> W. AI, (x), Y, (х)}= Jio (x) =J, 6) (x) = A box )- J (x)Y, (x) 


OO ONDES (ик) ONCE 


=> ч'(%)= LA) Jn), nc 
2 


= J (An Y (Ar,)=J (Ar, )Yo(Ar,) = B 
0 


Се qui GE de trouver le résultat final : 

J Gn), (Ar)- Y (Arn, J, (Ar) Cosh(À z) 
AA) + (HA) ` Cosh(A zo). 

Pour le problème suivant, le résultat découle de suite : 

AT(r,z)=0 rel[r,,œ) ze РЕА Т(+оо,2) fini 

T(r, 2) = 220 ТО» 2) =T; 


T(r,z) =T, 2] dÀ 


T(r,z) - T, -T (r,z) 
AT (r,z) 20 r efr,,00) zel-2,,2| T(+œ,z) fini 
T(r, 2), — =T, T(nz),, =0 


Ee ЫЕ Y (Ar, )J (Ar) Cosh(A z) 
AU, Qn) +(Y (Aro) — Cosh(A zo) 


—T(rz)-T,1 2 je J (Ar,) Y (ar) Y (Aro (ar) Cosh(A z) 
ia AU.) +(X (Ara) — Cosh(Az,) 


= Dr. 3) =É Li 
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Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation 
de Laplace : problème de Dirichlet axial électrostatique dans l'espace plan de hauteur 2z, 
potentiel ай à une charge placée à l'origine des coordonnées 


Dans ce cas la charge provoque un potentiel de la forme : 


Е To _ D 
1,(r,z) Je 2 2 Jr? 2 
X +) +2 r^ +Z 


La charge induite sur les plaques du conducteur parfait annule le potentiel sur ces derniers en 
z=+z0 et z=-z0. Le potentiel se décompose en deux sous la forme : 


T(r,z) - T, (r.z) + T,(r,z) = = 
r +Z 


Le dernier potentiel est solution du probleme aux limites homogènes sur les plaques hautes et 
basses, à savoir : 
АТ (r,z)=0 re[0,0) ze|-z,z] T(+œ,z) fini 


To To C(2)=-T, Í dr SC 


+Т,(т,2) 


T (r,z) 


== a 2 2 T(r, SÉ 20 2 2 
r +Z Y tZ t? +Z 


Sinh(A(z, + z)) + Sinh(A(z 
Sinh(22z,) 


— +J +e] nep) mie): Re(v)>-1 


Cosh(A z) 
Cosh(A z,) 


T.(r,z) = | dÀ AC(A (ar) 072) T.(r,z) = Í dÀ AC(A (ar) 


Comme Je e "J (cx) = 

үр? + 

с. =z 1 

=> fae J (rà) === ———= 

ур +c р? + z? r^ + z? 
-zÀ 


1 е 
Par application du théorème de Hankel = | dr r _ === J, (rA) = — 
j JE ss ooo À 


= ja e ™J (cx) = 


r +z 
+œ -zoA —ZoÀ 
= E rJ, (Ar) _е eu eum t 
0 EE À À 
+оо ` Cosh(A z) Т, u Cosh(2 z) 
Т, =-T, | dÀ e *^ J (A T(r,z) = -== - T, | dà e J (Ar) 2—5 
¿ една) o ren prn [e ARES 
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Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation 
de Laplace : problëme de Neumann axial électrostatique dans de demi espace z>0, potentiel dü à 
un disque de rayon r, posé à la terre appliquant un courant continu 


Dans ce cas il s'agit d'un problème inhomogène de Neumann en z=0 : 


Cette fois sur le disque de rayon r0, la condition aux limites est inhomogëne de Robin : 
АТ(ғ,2)=0 r e[0,0) ze[0,+œ0] T(+œ0,z),T(r.+œo) fini 


T(r,z) = Í dA AJ(Ar)C(A)e * 
| ÔT(r;2) 0 


Oz 


C.L. = f(r)2 


C(A)= : [ar rJ (Ar) f (r) 
0 
Si la fonction limite est de la forme : 


ПОРЕ r «n 


à agdas P [dr (дг) e [d 22002) lie) 
f(r)=0 к>» о о 
An +œ 
_ T, 17 _ T, А _ Tolo À -4z 
c(A)- mo | dz zJ (z) = D u^ (An) > T(r,z) = 2 | 4А =з J (Ar), (Аде 


E душ Yu Í 4А ` J (A) I ne 
H 0 


Le profil sur l'axe z change de forme : 


+00 -Az 
Т(0,2) 1 fan J (Ar )e 


0 


T 7 


L'intégrale dans les tables de A.P.Prudnikov, Yu.A.Brychkov еї O.I.Marichev- < Integrals and Series, 
Volume 2 - Special Functions > page 183, section 2.12.8, formule n°4 donne par identification des 
termes à la valeur a=0, il vient : 

Junk" p e 


2 
z+ z^ en 


pez vel gen += [4 
0 


+œ -Az 
Т(0,2) 1 [O PR DA Z 1 5 1 е 
0 À z+ + r 2z 


Regardons la répartition de potentiel dans le cas limite oü r0->0 : 


T, H 


+œ —Az 
T(r,z) = P» faa 7 MONDE pan) s 25 
5,5 A 
sieja Eeer 
0 2 2 
SE Sal «y 
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Méthode de résolution des problèmes aux limites de l'équation de Laplace en coordonnées 
cylindriques 3D sur des secteurs d'angle cylindriques non bornés par les transformations 
intégrales de Kontorovich-Lebedev. Problèmes de Dirichlet inhomogène angulaire sur des 
secteurs d'angle illimité en r, illimité ou limité en z 


On aborde en quelques pages les problèmes aux limites cylindriques sur des domaines non bornés 
radialement et partiellement bornés sur l'axe z, en forme de secteur cylindrique, auxquels on 
applique des conditions aux limites inhomogènes angulaire. 


La plupart des exemples donnés ci-dessous sont tirés des deux ouvrages d'exemples et d'exercices 
de Physique Mathématique : N.N.Lebedev, Special Functions and their applications, 1965, section 
6.5 «The Dirichlet problem for the Wedge >x, N.N.Lebedev, I.P.Skalskaya et Y.S.Ufliand, 1965, 
« Problems of Mathematical Physics », section 5 « Integral Transform Involving Cylinder Functions 
of Imaginary order » pages 194 et suivantes. Les énoncés ont été parfois modifiés pour ne faire 
apparaítre que la formalisation mathématique du problème aux limites de Laplace, sans se 
préoccuper de la forme exacte des grandeurs physiques en jeu. 


Le probléme de Dirichlet inhomogéne angulaire est le suivant sur un solide cylindrique non borné 
en r et z, seules les conditions aux limites angulaire suffisent à le définir. 


C.L. T(r.0,z), = f} (r,z) T(r,0,z)|,_,, = f$ (r. z) 
Par la suite on définira d'autres problëmes sur des domaines non-bornés radialement avec des 
configurations géométriques plus complexes toujours sur des secteurs d'angle. 


La séparation des variables dans les trois coordonnées cylindriques donnent trois équations 

differentielles ordinaires, et les solutions de chacune, notamment les solutions radiales sont alors 

dans la catégorie des fonctions modifiées de Bessel d'ordre purement imaginaire, et dont la partie 

imaginaire de l'ordre est non nécessairement entière : 

d*T(r,0,z) x 1 ôT (r,0,z) d 1 eT(r,0,z) " eT(r,0,z) 
or? r Or г? 00? дг? 


EE с v/e(8) - 0 EE 0 


= 0 = T(r,0,z)= R(r)O(0)Z(z) 


dr? r dr 
Solution du type R(r)=AQJ„(iur)+B,Y, (iu r)= AJ, (u r)+ B K, (u r) 


Fonctions de Bessel modifiées d'ordre imaginaire 

@(0) = A,Cosh(v 0)+ B,Sinh(v 0) 

Z(z)= A,Cos(u 2)+ B,Sin(u 2) . 

Nous savons également que les fonctions de Bessel modifiée de premiére et deuxiéme espéce 
présentent une singularité en zéro par une oscillation rapide qui les rend impropre a étre utilisées 


dans des développement en série. Par ailleurs la fonction de premiére espéces I doit étre écartée 
car elle divergent lorsque r tend vers l'infini. La solution de base est donc la suivante : 


T(r,0,z)- K, (u r\4,Cosh(v 0)+ B,Sinh(v 0)(A.Cos(u z)+ B.Sin(u z)) 
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L'idée de Kontorovich et de son élève à l'époque N.N.Lebedev, est de tenter d'utiliser cette solution 
de base : 

T(r, 9 ,z)> Т, „(r,0,z)= K, (u r {4,Cosh(v 0)+ B,Sinh(v 0)XA,Cos(u z)+ B Sin(u z)) 
non dans une série, mais dans une intégrale sur le paramètre de valeur propre. A cet effet 
supposons que les fonctions limites sont paires en z, et admettent une représentation intégrale de 


Fourier-Cosinus : 

FEM DO el r ‚и)Сов(и z)e g^(r, u)- — 2 [az ss (r ,z)Cos(u z) 

Si la See est impaire alors, on utilise la transformée de Fourier sinus : 

Kat? )= [ди el r ‚ и)5їп(и z) gr. u)- — ^ (az j; Ío (r ‚2)5їп(и 2) 

Et si la fonction n'est ni paire ni impaire, оп sépare le problëme еп deux sous-problëmes avec 
respectivement les fonctions limites suivantes : 

ach EE pepj EZ 


On cherche alors la solution du problème sous la forme d'une double intégrale (sur le spectre 
radiale et angulaire, ici prenons une fonction limite paire) : 


T(r,0,z)= Jaa Cos(u gd K, (u rja (v.u) t ү E s) POH) SR = 


L'application des conditions aux limites de Dirichlet donne : 


T(r.,0, z) = = [du Cos(u z fav K, (u r)G, (v, u)- jau g? (r, u )Cos(u z) 
T(r, Ө, z) = = fau Cos( Hz Ju K, „(и r)G,(v, u) = (du g2(r, u)Cos(u z) 


gi (r.u)= fav K, (u т)б\(у, и) 


aile, ul [av Kp (u r)G:(v,u) 
0 


Dës lors la solution dépend des valeurs des fonctions G1 et G2 qui sont les coefficients de 
l'expansion intégrale par le поуаи de la fonction K d'ordre imaginaire. Cette expansion intégrale est 
possible grâce à une théorëme dü à ses deux co-auteurs, Kontorovich et Lebedev. Le théorëme 
permet de calculer les fonctions G1 et G2 par le biqis d'un transformation intégrale (et son 
inverse). La solution formelle du problème est alors la suivante (voir le théorëme de la 
transformation intégrale) : 


p =1,2 g,(r,u)= 2 јал f$ (r, z)Cos(u z) G, (v.u) = ту Sinh v) | ar $ gp (r, ик, (u r) 


r 


EE EE 


Résolution des problëmes aux limites de Laplace en cylindriques 3D par les transformations de Hankel et les 
transformations de Kontorovich-Lebedev- p 212 


Théoréme de Kontorovich-Lebedev avec les fonctions de Bessel modifiées d'ordre imaginaire 


également appelées fonctions de MacDonald 


Pour toute fonction f(x) sur l'intervalle [0,*e9) : 
1 - continue par morceau et dont la ашп est bornée sur tout intervalle [x1,x2] de [О,+оо] et 


2 - telle que les intégrales sont finies : ja үө, ER Jes et farw dn 
X 


2 
Cette derniére admet la représentation intégrale suivante, pour x50 : 


f (x) = Z [av v Sinh(z v) ell [at РО) = 2 
Au point de j E de la fonction la formule intégrale est la suivante : 
f(x+0)+ f(x-0) _ Km 
i A 


Ï dv v Sinh(z v) no [at f) = 


Par le changement de variable, et le changement de fonction : „ _ ur t=up Vx f (x) с g(u r) 


La formule к devient : 


Vx f(x) = g(u,r)= Z [ду v smile v)K, (u D [de fu Piu р = 
Ju pf (и p) = 2(u. p) g(u,r)= Zi v Sinh(z v)K,, (u r) Í dp et Pa) x p) 
En — l'expression de la А. g(u,r ES g(r ) 
= g(r)= À [av v Sinha v) n Í dp et аши K; (u р) 
p 
Si rh(r)= g(r )= rh(r)= 3 fav v Sinh(z v)K, (u n [dp H hlp )K,, (u p) 
= h(r)= 2 [av v Sinh(z y) кыши” faon h( pK, (u p) 


Et les SE de validité du EE intégral est assurée par la finitude des intégrales 
suivantes : 


nd н jet 


йй as nm ee 


ССВ ас cae fa ED cœ а (r) <œ 
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Les trois к intégrales : 
Ka (®© 


am 


g(r)- 2j fav v Sinh(z v )K,, (u r) Í dp ыы 


F) =Ê A dv v Sinh(z v) EO, | dt f (t) x8 


n(r)= le v Sina v) ED | dp Mel, (a p) 


peuvent se décomposer afin d'introduire les diverses formes des transformées et transformées 
inverses de Kontorovich-Lebedev des fonction f, g, h 


_2{ | K (x) Ku 7 
fe--— | dv v Sinh(z v) = | dt f (D) T fx) S fv) 


"f dt f (t) EE ) transformée de Kontorovich — Lebedev 


-() 


K ` (x) == | dv v Sinh(x v)f(v) = transformée inverse de Kontorovich — Lebedev 
| п х 


g(r)=Z [dv v Sim(z v)K„ (ur) [dp (о) EP а) = ge) 


Ks fap g(p)K,,(u p) 


transformée de Kontorovich — Lebedev 
= (2) \ 
S 2 +œ E 
K " (r) = = [4 Vv Sinh(x v)g (v )K „(и r) transformée inverse de Kontorovich — Lebedev 
0 


uc fav v Š 200 ke UD fap Mal, Qro) Moien RW) 


K, (v) = [do g(p)K,, (up) transformée de Kontorovich — Lebedev 

29] T | 

К r" Dis [4 v v Sinh(z v )h (E €9 transformée inverse de Kontorovich — Lebedev 
T $ r 
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Exemple : le développement intégral de la fonction suivante, relève d'une transformation de 


Kontorovich-Lebedev qui permet d'associer la fonction à sa transformée : 


ЖЕЕ, , Sinh( а 2 КОО)... gel 2 = Sinh(a у) 
= С Ga dv iv ( ) = d | i 
Ess gilt Sin(a) Vx Së л | =; Ѕіп(о) w) 


Par comparaison avec la première formules de la transformée inverse, on en déduit ceci : 


fG) = fem 2 | dy Sinh E 


Silo Vx 
T" Dm | K; œ E CTAR Sinh(o v) 
R; œ => [av v Sinks vr) ro^ => f (s) Jxe t f) e < T batter 
г m NT u Sinh(a v) 
= ja EE ()- [а < K; (0) =  Sin(a)sinh(z v) 


— ja e KE (t) 2 x Lim Sinha v) = — TY 
Á a>0 Sin(a)Sinh(xv) Sinh(xv) 
En prenant un paramètre “а a->i a : 
Sin(a у) 
Л 
Sinh(a )Sinh(z v) 
Prenons le б М еп N de Fourier-Cosinus de la fonction de MacDonald, selon son 
ordre v (voir la table des transformations de Fourier de F.Oberhettinger, « Tables of Fourier 
Transforms and Fourier Transforms of Distributions », édition 1990, section 1.20 Modified Bessel 
Functions of Variable Order, formule 20.5, page 100), on гесоппаїї là une formule de transformée 


intégrale de Kontorovich-Lebedev, et en utilisant l'une quelconque de ses formes, il vient : 


„ Sin(a v) +œ 
осна) _ 2 git dv y 27 = ` K, (x) = [йе "9", (t) = 
0 


K, 
BEE )- [arro ESO CR 
м, = [dv Cos(va) E) => 
Ух 7 ° Ух Р [dv v Sinh(z ve) e _ 2 far Cos(va) 6) 
f л? d Ïx л> Jx 
TEA л Cos(va) 
Sek v Sinh(z v) 
2 far e~ rla) p AD _ л Cos(va) far e CR (т) u л Cosh(va ) 
A t v Sinh(z v) E t v Sinh(z v) 
IC ek (t) л 
а з = 
SS | { t v Sinh(z v) 


Ces traitements sont non nécessairement valide voir plus bas 


Toutefois, les conditions d'application du théorëme de Kontorovich-Lebedev sur la finitude des 
intégrales ne sont pos vérifiées dans ce us : 


=x Cosh(a ) e * Cosh(a 
ro =C >ja 77, d E ja: = reef 1) › œ 


Mais en revanche pour la première remm les conditions sont remplis : 
1 1 


2 2 
J (x) aS f ar L roef) = [ax e^ SCH < œ 
0 Ух x 0 x 


Q 
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De тёте que SSES la fonction 


Тое 009 = ја ө, «| )- EE 2 — Log )<® 


Pourtant, on peut utiliser la dérivation par au paramètre pour retrouver des intégrales déjà 
calculées auparavant : 
+00 —t Cosh(a IK. +œ —t Cosh(a к. B 
[а £ ON л Cos(va) => -Sinh(a ell ac? „(О _ _ xv Sin(va) 
А t v Sinh(z v) t v Sinh(x v) 
л Sin(v a) 
Sinh(z v )Sinh(a) 


= pue (Сота) (f) = 


[ar P KA 2 Liri л Sin(va) | xv 
A a>0 Sinh(z v )Sinh(a ) Sinh(x v) 


Autre exemp le : D'aprës l'ouvrage A.Erdelyi, H.Bateman, Higher Transcendental Functions, volume 
2, page 55, formule (42) 


KY fe +b? — 2abCos(9))= 2 fav K (a)K,,(b)Cosh(v(x —0)) 
л 0 


En posant la valeur Ga, il vient : 


KY Ма +b? + 2ab]- Ki (a+b)= _ [dv K,,(a)K,, (b) 
л 0 


= Kë (a + x) = — [dv K,,(a)K,, (x) 
7 


En identifiant l'expression comme une transformée de Kontorovich-Lebedev, il vient en utilisant la 
premiére formule : 


= fa ro E 


K.(a-x) 2*7 K, (x) 
————— |dvK, (a) => 
Vx T | VM a x Se Ze => K. (x) 27e K, (x) 
Kz; ol Ta v Sinh(z v )f (v) EE ja K, (a) SE 
Ра == z K, (a) Re ES (a+t)K,, (t) == TK, (a) 
= dk 2 Sinh(z v) v EZ e t су Sinh(x v) 


Ko (OK, _ x K,, (0) 


d 
Ei Í f t v Sinh(z v) 
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Autre exemp le : D'après l'ouvrage A.Erdelyi, H.Bateman, Higher Transcendental Functions, volume 
2, page 55, formule (42) 


Si on revient à la premiëre expression : 


к [Va +x? – 2axCos(9))= 2 f av K,,(a)Cosh(v(x -0))K,, (x) 
л 0 


Par identification des termes de la transformée de Kontorovich-Lebedev, il vient еп utilisant la 
premiëre formule : 


EK (Ja? +x? — 2axCos(9)) _ 2 [av K,, (a)Cosh(v(n -0)) = (x) 
л 0 


Vx Ух 
R v)= [ало 60 
= 0 
Rj) 2. [av v зт reyes = 2 fav K,,(a)Cosh(v(z 0) K.) 
Can Ne T Ki, (a)Cosh(v(z -0)) _ `€ Ки (Ja +t? – 2atCos(9))K,„ (0) zK, (a)Cosh(v(z — 0)) 
AVE v Sinh(z v) = | t Е v Sinh(z v) 


Autre exemple : démontrons la formule intégrale : 


wg ( e e) соз" | —1 

Je Er, (her 2 

^ x v Sinh(x v) 

En prenant pour cela la première formule de la transformée inverse de Kontorovich-Lebedev, il 
vient : 

TV 


v 135 к) = : jo [con 23 E J кы) 


соң 


(i-e) ne fav v Sinh(z vyr 
0 


Jx" e 


S l-e "= Z [dv СЕВ) (х) Comme e = E [4v K,, (x) 
л 7 2 TY 


2 +оо 
Cela revient à démontrer que Mir | — [4v SEA (х)= 1 
л 0 


Cette dernière assertion est obtenu à l'aide de la formule intégrale (voir la table des 
transformations de Fourier de EOberhettinger, < Tables of Fourier Transforms and Fourier 


Transforms of Distributions », édition 1990, section 1.20 Modified Bessel Functions of Variable 
Order, formule 20.5, page 100), qui donne par la substitution a-»ia : 


gemi) _ Z [ау Cos(va)&, (x) 2 е" te) = 2 [dv Cos(via XK,, (x) 
л 0 л 0 
=> gon). 2 fav Cosh(va )K,, (x) oe => 1= 2 [dv SEA? (х) 
л 7 2 T 2 


Ce qui démontre l'assertion et par suite la valeur de l'intégrale à l'aide de la transformée de 
Kontorovich-Lebedev. 
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Autre exemple : démontrons la formule intégrale suivante 
Cost( Z | — Cos(va) 


[ыш 
0 


v Sinh(z v) 

Supposons que la formule s'écrit à l'aide de deux inconnus : 

v 
+оо [= eem) acos "У )+ B 
fa I, G)-—— 
^ x v Sinh(z v) 
Nous avons alors : 

TV 

ACosh| — |+ B 

(gra) 2 +оо | 2 | К. (х) 2 +оо лу К. (х) 
= dv v Sinh iv = dv | ACosh| — |+ B |—" 
Je sc) vv Sin (z v)z SAV) ке; Al v озң Z |+ e 


< 1— e we) = zi [4v КЕЗ A ei, (x) 
ds, 2 
Comme e *Ces(e) = 2 [dv Cos(va)K,(x) et 1= 2 [dv SEA (х) 
л 7 T 0 2 
-— ]- e *Cosh(a) = 2 [ау 63656 (х)= > [dv | cosh 8 )- cosa) к, (x) 
ж 2 л e 2 
=» А=1 et B= -Cos(va) 


On а donc établit les deux formules suivantes : 

T ( Le 9) Cosh( Z ) — Cos(va) 
[ах K, (x)= я 

0 " v Sinh(z v) 
Cost( Z | — Cosh(va) 


^ x | v Sinh(z v) 
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Exemple d'une transformation mixte Kontorovich-Lebedev et Mehler Fock 


Dans la table sur les transformations de Kontorovich-Lebedev et Mehler-Fock F.Oberhettinger, 
« Tables of Lebedev, Mehler and Generalized Mehler Transforms >, la transformée suivante de 
Kontorovich-Lebedev nous est donnée en page 10: 


[dv v Termes | K, ( (s)= [se MISE Te |— ? [av v Tanh(xv)P , (a)K, (x) 
d кыл, TX а 


P, (a) 
2р л? ну K (x) [z . 
=> |dv v Sinh - ET e 
л? | ке 2 Cosh(zv)] Vx s 


En inversant cette intégrale vue sous la forme d'une transformation de Kontorovich-Lebedev, il 
vient, une représentation intégrale de la fonction conique de Mehler de première espèce : 


f(x) =Ê al dv v Sinh(z v) Eo | dt f (D) Â. Ро) = fv) 


V = [а foy tee 


transformée de Kontorovich — Lebedev 


K " (x) = а Í dv v Sinh(x ore) transformée inverse de Kontorovich — Lebedev 
| л?! 
(а) 


х 
Р Р, (а) 
— л? Liv c pu c K. (t) л? bay 
Ici NAR Ne — = jt ax d -a x Div = 2 
SS SE 2 Cosh(xv) FN 2 ds | i \2 I Л 2 Cosh(xv) 


3 +00 
>P, (а) = V27 ?Cosh(nv) [at e"? K; (t) 


——+iv 
2 Í 0 [4 


Vue sous la forme d'une transformation intégrale de Mehler-Fock, il vient une représentation 
intégrale de la fonction de MacDonald : 


fec)» (dv P, (s)vTanh(av) [dy fP, (у) = SOS f» 
М, (v) = fv) =v Tanh(zv ) Í dx ҒОР, | (х) transformée de Mehler — Fock 


M; (x) = f(x) = fav f(v)P, (x) transformée inverse de Mehler — Fock 
H 5 SEH 


+оо 


Сотте fav v Тапй(лу)К\„(о)Р , (x)= Ce e * => f(v)=v Tanh(zv)K, (a) => f (x) = PA e 
— таня), (a) Йаа) ўа [Es еер, ()oK, 9-2} jeep, (9 


Et l'on retrouve cette formule en page 23 du même ouvrage, F.Oberhettinger, « ТЕ of jueces Mehler 
and Generalized Mehler Transforms », la transformée suivante de Kontorovich-Lebedev. L'examen de ces 
tables permet de déduire toute une série de représentation intégrale liée à la présence mixte des noyaux de 
ces deux transformations intégrales : Kontorovich-Lebedev et Mehler-Fock. 
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Revenons qu problëme de départ, en utilisant la deuxiëme formule des transformées de 
Kontorovich-Lebedev : 


ke gle)K,, (u p) 


š transformée de Kontorovich — Lebedev 


0 
K; (r) = = [ау v Sinh(z v)g(v)&,. (ur) transformée inverse de Kontorovich — Lebedev 
л 0 
gs (r.u)= [av С„(у,и)К„(иг) pour p=12 
0 


2 +œ 2 +œ 
= ejr. u)- = [dv v Sinh(z v)8(v)K„(u r)= [dv G,(v.u)K,,(a r) 
0 0 


elo, ul, (u p) 


=> G,(v,u)= EL Sinh(z vie v= Zv Sinh(z v) Í dp 
0 


La solution formelle du problème est donc : 
CL. T(r,0,2),, = 30,2) Т(,0,2),, = f26.2) 


gï (ram) =Ë Í dz f?(r.z)Cos(u z) (м) = 2 |а f(r.2)Cos(a z) 


+оо ө 


+œ 9 
Gv.) = Zv gute v) fap PR UP) Gi) Zv иң v) f ap Rm 
9 0 


T(r,0,z)= [du Cos(u z) [av К (и 2 7; 
0 0 
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II arrive souvent que les fonctions limites ne permettent pas de garantir la finitude de la premiëre 
l'intégrale pour l'application du théorëme de Kontorovich-Lebedev : 


[а ERA Gel 2) GEN 
0 


Dans ce cas on régularise le problème еп introduisant des fonctions modifiées : 
s; (2 и) z ei (r, u)- ei (0, uk, et l'on assume gue les conditions de finitude de l'intégrale sur la 


= 


fonction modifiée est assurée, soit : 


E s 
far n) ll Ech 
d r 


r 


On résout alors le problème avec ces fonctions modifiées : 


+œ 


A n š A 2 | 
g,(r.u)= [av Gv, и), (ur) G;(v.u)- <> V Sinh(r v) Í dp 
0 


0 


g,(p.u)K,, (u р) 


Il existe un développement connu de la fonction exponentielle en fonction de MacDonald : 
2 +оо 

e = — [dv KE. (x) 
л 0 


Pour le prouver оп utilise la transformée de Fourier-Cosinus sur ип ordre variable de la fonction de 
MacDonald : 


e Coh) — = [ау Cos(va)K, (x) > a =0=>e" = 2 [4v K, (x) 
л % 49 


Ce qui implique : 
T 2 +œ 
gp Ou)" =< g7(0, n) f dv K, ur) 
0 
Ce qui permet de revenir facilement aux termes G, de la maniére suivante : 
e, (r.u)- [av G; v. u)K, (u r) 
0 
g2(r.u)= [av G (v. u)K, (u r) 
0 
Bi 2 +œ 
g$ (0, u)e""" = < gp(0,u) f dv K; (ur) 
0 
А ЕТЕ 2 +œ 
ar. u)= g, н) a5 uk!" = aser n)- — gp (0, и) [dv K„(u r) 
0 


+00 , +оо 2 +œ 
=” [dv G (v. и)К „(и r)= [dv G (v.u)K,, (u r)- -25(0. 1) f dv К „(и r) 
0 0 0 


: 2 : 2 
=> v.ai)s G,(v. u)- eines G (v. ls G;(v. и)+ ein 
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La solution dans ce cas s'écrit maintenant : 
CL T(r.0,3)),_,=f^(r.z3) — T(r.0,3)),_, = C2) 


gr. u)- SÈT SE (r ,z)Cos(u z)= g2(0, u) u ) = - 2 |а 7; (0, z)Cos(u z) 


B. lr u)- 2 fae se Í, (r,z)Cos(u z)- e Í dz E (0, z)Cos(u J = 2 а (2 (r,z)- e" n (0, z))cos(u z) 


0 [U 


i + 
tel que | dr Sal Log 1 )< o р=1,2 
d r r 


А Jo, u)K, (u p) 2 
G;(v.u)= EL Sinh(z CH 81 et С,(у,и)= Gv, n) — 25 (0. u) 


T(r,0,z)= au Cos(u z) dv K, (u eu) Sinh(v (0, -0)) ot) Sinh(v v ur 
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Exemple : potentiel d'une charge q, placée devant un demi-plan conducteur à une distance a de 
l'origine. 


Le demi-plan conducteur parfait est dans le plan <x,o,z>, le bord du plan coincide avec l'axe z. 
On suppose que la charge est dans le тёте plan <x,o,z>, placé à la distance а de l'origine, soit en 
coordonnées cartésiennes q placé en (-a,0,0). 


Z 
2 y 
P Leg 
x " 
# 
Le potentiel crée par la charge q est le suivant : 
q q 
Т, (r,0,z) = = 
| Ale af Ais qr +а? + 2ах+у? +2? 
x = г Cos(0) 


4 
>T (r,0,z)= 
кн U Jr? + à? +2ar Cos(0)+ z? 


La solution est recherchée sous la forme d'une addition du potentiel crée par la charge q et d'un 
potentiel recherché : 


T(,9,2 7 1,6, 0,6, n LL OF -+ (r,0,z) 
r a ar Cos 2 


Т(т,Ө,2)]„=0 T(r,0,2)|,_,, = 0 
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Le potentiel T(r,0,z) répond alors au condition aux limites suivante : 


_ _ pfo = q SS d 
T, (r,0,z),, = T. (7,0. Si, = f; (r,z)=- 2 2 2 == 2 2 
Jr +a” +2ar + z (r +a) +Z 
La transformée de Fourier-Cosinus de la fonction limite est la suivante (voir la table des 
transformations de Fourier de F.Oberhettinger, < Tables of Fourier Transforms and Fourier 


Transforms of Distributions >, édition 1990, section 1.1 Algebraic Function, formule 1.32, раде 5): 


а) -24 а En) _ 22 кы, ш) 
Е. (r +a) +z T 


K¥ (x) fonction de MacDonald KA (x) = Lim К, (x) 
KY (x) =K;„(x) 
Dans ce cas la valeur en r=0 est la suivante : 


2q 
g0, u)=- Ky (ua) 
Et la fonction modifiée est la suivante : 


| 2 22 _,, 
gpr u)=- FK (ur a) е FRE (ша) 


ви) C Ky (ur a) e" Ky (u a) 


Posons nous la question : est-ce que l'intégrale suivante est finie : 


1 1 

2 9 2 M 

[а esr. u) Log] 24 | qy K, ta) Н 
° r SS r r 


r 


En utilisant par exemple Mathemqatica, on se convainc facilement que cette intégrale diverge : 
Clear[MacDonald]; 


MacDonald[0,x_]:=Limit[Re[BesselK[I v,x]],v->0]; 
Plot[MacDonald[0, x], {x, 0, 0.000011}] 


mu = 0.5; а= 2; 

Clear[IntMacDonald]; 

IntMacDonald[s ] := NIntegrate[MacDonald[0, mu (х + а)]/х Log[1/x], {х,=], 0.51]; 
Print[IntMacDonald[0]] ; donne 4.73288*10? 
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Posons nous maintenant la коре question : est-ce que l'intégrale suivante est finie : 


| lgi (r.e) Log} ETT (ку (u(r+a))+e *” KY (u 2770 


r r 


Avec le code Mathematica suivant : 

mu = 0.5; 

а= 2; 

Clear[IntMacDonaldModified]; 

IntMacDonaldModified[e ]:=NIntegrate[(MacDonald[0,mu (x*a)]-Exp[-mu x] MacDonald[0,mu а])/х Log[1/x], {x ғ, 
0.5}]; 

IntMacDonaldModified[0]--0.0667506 


Ces deux résultats justifient donc pleinement la régularisation opéréee par modification de la 
transformée de Fourier-cosinus de la fonction limite. 


Le coefficient G se calcule donc comme suit : 
elo, uk, (и р) 
p 


G HUS u)- Ev Sinh( (ту Ju 


* 4 7t к —е“*°КМ IK. 
es G'(v,u)= == Sinh(x v)[ap ( o (u(p+a))-e o (u a) „(и p) 
T ^ p j 
L'évaluation de l'intégrale est assez compliquée, faisant appel à plusieurs formules plus ou moins 
connue sur les intégrales définies des fonctions de MacDonald : 
ЕЕ fa (Kë (u(p a) - e" ^ K (u а))к, (и p) 
p 


© (Kyu a) e7 KP (и а))к, (х) 


x =U p => Í = |: = 


‚- Га KES нак) 
0 


— Ein да 09 
x 0 x 


Nous avons établit avec les quelques exemples données de transformations de Kontorovich- 
кеды (e I Е suivante s'évaluqit comme suit : 


[aK M(a+t)K, ,@) _ л K, (a) 
v Sinh(z v) 


II vient donc : 


+оо =Ë +00 M 
[a K, (O) _ 7e De même [ax É (x+ a)K,, (x) ` AE K (u a) 
^ t v Sinh(z v) " x v Sinh(z v) 
I = je Ko (x+u а)К„(х) vn (и а)[ах e"K, (x) 

0 x 0 x 

л M 
I =—  — IK, -K 

ES 2i „(и a)- KY (u a)) 


(K, (u a)- KY (u a)) 


. 4q Ln 49 | ç; Е 
=> Gi (v, u)= SR Sinh(x 2 = E Sinh(x KSE 
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Le résultat du calcul littéral de l'intégrale donne donc : 
G;(v,u)= “||; (и a)- K,, Qu a)| 

аби) = (и) 2.87 (0,и) = Gs и)-2® KY (u a) 
=G „(v, u) = -ZK, (u a) 

Ce qui donne la solution du problème : 


T,(r,0,z)= - A4 (au Cos(u z)fav К „(и 7) Ke Qt 0). [inh(v (2л — 0))+ Sinh(v Ө) 
л 0 0 


Sinh(2z у) 
4q T Y К, 
T,(r,0,z)= = Јан Cos(u Se K (ur) SE À [Sinh(v (2z -0))+ Sinh(v Ө) 


a=v(r-0) Ь=лу 
Sinh(v (2л — 0))+ Sinh(v 0)= Sinh(b + a)+ Sinh(b — а)= 2Sinh(b)Cosh(a) 
= 2Sinh(z v )Cosh(v(x -0)) 

e Cosh(v(z —0)) 


Сдвй\т v) K, (u г)к,, (и а) 


4 +оо +œ 
T,(r,0,z)= = [du Cos(u z)[4 
0 0 


On pourrait utiliser cette formule: K“ (u(r + a)) = 2 [av KE (ak, (u r) 
л 0 


Ainsi qu'il est établit dans l'ouvrage A.Erdelyi, H.Bateman, Higher Transcendental Functions, 
volume 2, page 55, formule (42) mettre la valeur Ga. 


KY (Ja? +b? - 2abCos(9))= 2 fav K,,(a)K,,(b)Cosh(v(x — 0)) 
л 0 


Mais en fait оп a mieux à faire, en intervertissant les intégrations de la solution : 


Colle —0)) к („ук Qi a) 


4 Фу +00 
156.» - S fau costa fav Ch 
0 0 


4q `F „ Cosh(v(x-0)) F 
БУ К, K, 
л? E Cosh(z v) [аи Cos(u z)K,, (u т)К „(и a) 


Regardons maintenant la transformée de Fourier-Cosinus du produit des fonctions de MacDonald, 


et l'on trouve dans la table des transformations de Fourier de F.Oberhettinger, < Tables of Fourier 
Transforms and Fourier Transforms of Distributions >x, édition 1990, section 1.17 Modified Bessel 
Functions of Arguments x, х? and 1/x, formule 17.23, page 90), la formule suivante : 

= л? а? +b? + у? 
17.23 ах Соѕ\х у)К (а х)К (Б y) 9 —————— P EU EC 
(17.23) fax Costs IK, (аков у= uer, ts 
viv xu br yz aa 


"m л? а? +? +2? 
du G K. K. = P = ласы = 
| и os(u z) „(u a) „(и r) ACosh(n v) Jar il 2ar | 
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4q ©. Cosh(v 
T,(r,0,z)= — „a Jav GH Í du Соз(и z)K,, (и r)K,,(u a) 


Cosh(v Gë 0)» É +r? +z d 


>F [av SCH лу) 2ar 


Maintenant c'est là ой les coïncidences font bien les choses, dans la section du présent document 
dédié qux coordonnées toroidales, dans l'exemple : < solution du problëme de Dirichlet sur le 
domaine extérieur d'une coupole ou coquille sphérique », on a démontré une formule extrêmement 
importante donnée dans l'ouvrage en russe < A.P.Prudnikov, Yu.A.Brychkov, O.I.Marichev, Integrals 
and Series - Volume 3 - More Special Functions» en page 181 formule 2.17.24.6 (mais la formule 
contient une légère erreur dans le deuxième terme à droite qui est corrigée dans ce qui suit), et 
que nous avons également démontré quelques pages plus loin : 


‘© Cosh(r b) 1 2 1+Cos(b) 
dr — — — P = 1-=A — 
| š Cosh° (лт) 5+ (e) 42 4c - Cos(b) | T zi c — Cos(b) ) 


= T,(r,0,z)=- 


ArcCotan s 
“© Cosh(r b) /2 1 л 1+ Cos(b) ll V2 c—Cos(b) 
> [dc =P, | (с)= Arctan = 
^. Cosh’(rr) Zu z dc-Cos(b)| 2 c — Cos(b) л c — Cos(b) 
ArcCos 1+Соз(ф) 
je Cosh(r 2P ( V2 1+с 
d Cosh° (лт) -yur T Jc - Cos(b) 
Appliquons cette formule à notre cas, il vient : 
** Cosh(r b) 1+ Cos(b 
ат —— — P A F = > 
| Сози (лт) СЕ i.e — Lal re E 
Tv b—m-0 ui Cos(z —0) - —Cos(0) 
2ar 
+оо 2 2 2 
[ау Cosh(v(z —0)) É Jr +2 )- | - En : T= Coste) 
^ Cosh? (x v) —+iv 2ar ү "n + Cos(0) л a +V +Z + Cos(9) 
2ar 2ar 


2ar gc Ban el з че! | 


Jada! ep +2 2arcos a^ +r? + z? «2arCos(0) 


Jar = arcra Bar " = Coste) | 


Г Ма? +r? ez +2arCos(6) a? +r? +2? +2arCos(0) 
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et la solution devient, aprës ces divers calculs : 


T(r.0,2)=>-4 ER (- : arci zar | m | 


Jar Ja? +r? ez +2arCos(0) а? +r? +2? +2arCos(0) 


Т,(,0,2)= 4 9 - Е arran zar | 1-Cos(9) | 


ш 2 2 2 
Ja? +r? +2? +2arCos а? +r? + z2 +2arCos(0) 


or 21-cos(o))=asm'(2) 


Et oh ! Miracle, le résultat final s'exprime d'une manière plutôt simple : 


д 2 arsi 2) 
T,(r,0,z)= — q 1— — Аксап 2 
Jr +a? + 2arCos(9)+ z? л Jr +a? + 2arCos(9)+ E 
q 
T(r,0,z)= +T(r,0,z 
) Jr + a? +2arCos(9)+ z? À ) 
2Varsin( 2) 


= T(r,0,z)= 24 Arctan 
луг? +a? +2arCos(9)+ z Jr? +a? +2arCos(9)+ z? 
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Exemple : potentiel d'une charge q, placée devant un demi-plan conducteur dans le plan <x,o,z> 


à un point quelconque dans le plan <x,o,y> externe qu conducteur, aux coordonnées polaire r=r, 
et 9= 0, 


Le demi-plan conducteur parfait est dans le plan <х,о,2>, le bord du plan coincide avec l'axe z. On 
suppose que la charge est dans le plan <х,о,у>, placée à la distance го de l'origine et d'angle O,, soit 
en coordonnées cartésiennes q placé en (roCos(8;),rOSin(8;),0). 


Z 
у 
2 p 
X 0 
$ оо q 
Le potentiel crée par la charge q est le suivant : 
Position de q (r,Cos(0, ) rSin(8, Lol 
T,(r,0,2) = R = а 
JG — ryCos(0, y + (y — rSin(0, y + z? Jx +y? + n - 27, (xCos(6, )+ ySin(0, )+ z? 

x = г Cos(0) 
y-rSin(0) => T,(r,0,z)= 4 


- Beer J + n =2rh (Cos(0)Cos(0,)+ Sin(0 )Sin(0, ))+ z? 


q 
=> Т (7,0,2) = 
: Jr? +n? - 2r nCos(0 —0,)+ z? 
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La solution est recherchée sous la forme d'une addition du potentiel créé par la charge q et d'un 
potentiel recherché : 


ИКЕ жох 0,)+ ИЕ 
0 — 0 -Vo)TZ 


T(r,9,2), =0 T(r.0,2), , = 


0=2л 


Le potentiel To(r,O,z) répond alors au condition aux limites suivantes : 


T,(r,O,z), | = T, (r,O,z) = filr, z)-- q 
5 Ke m =2л т -2rnCos(8,)4 z 


La transformée de Fourier-Cosinus de la fonction limite est la suivante (voir la table des 
transformations de Fourier de EOberhettinger, « Tables of Fourier Transforms and Fourier 


Transforms of Distributions », édition 1990, section 1.1 Algebraic Function, formule 1.32, page 5): 


ele, ul =— ed [а = = Cos(u z) = => 24 ку (ur + n —2r r,Cos(6,) ] 
to Jr + -2rrCos(0,)4 z T 


К (x) fonction de MacDonald KA (x) = Lim K;, (x) KY (x) = K,, Gi 


Dans ce cas la valeur en r=0 est la suivante : 
2 

g2(0, u) = . RF (ur) 

Et la fonction modifiée est la suivante : 


вн) 22 Ky (dr? +n? -2rrCos(0,))+ e K$ (ur) 
л 


ҮК кс et 2) 


Le coefficient G se calcule donc comme suit : 
8,(0,и)К,„ (и p) 


С (у, и)= 2. v Sinis v) јар 


. м vg, UR (d en Zenter "Eil, p) 
= С (у,и)= 5" Sinh(z v) Í do 


p 
0 
T, e we (O) л M ose AR (ир) M л 
d м = K ар *”* =K ERES MENS 
| i t v Sinh(x и Wa nf 4 p A n), Sinh(x v) 


4q | v Sinh(z “р Ко (Ау р? «n -2p raCos(0,) Ka Qu p) 


2 


=> G,v.u)-- 
л л A p 


-К (u n) 


^ Ky (d Up run -2upy rCos(0,) к, (up) += Ky (J р? tum -2pu rCos(6,) к, (р) 
dp = | dp 
Ó p d p 
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Nous avons établi dans les quelques exemples des transformées de Kontorovich-Lebedev : 
Comme 


ja Kò (Ja? +12 —2atCos(9) JK, (t) _ л К (a)Cosh(v(x -0)) 


pot 


d t v Sinh(z v) 
rz ку (Je + u? ro -2t u rCos(0,) ls, (4) z K, (u " )Cosh(v(z: ө, )) 
Е | | t x v Sinh(z v) 


Le résultat du calcul littéral de l'intégrale donne donc : 


G (v. u) => "2 [K (un, )costiv (s 0) - ке (u r.) 


А 2 * 4 
Gv u)- G; vsu) 7 g, (0,u)= Chal Zei (ur) 
4 
>G,(v.u)= - SK, (ur)Cosh(v(r -6,) 
Ce qui donne la solution du problëme : 


T,(r,0,z)= -44 [аи Cos(u z) [av K, (u r) K, (u r,)Cosh(v(x —0,)) 
0 0 


[Sinh(v (2x — Ө))+ Sinh(v Ө) 


Sinh(2x v) 
Ti(r.0,z)=-24 Í ди Cos(u 2) dv K, (a r) кш m eu [Sinh(v (2x —Ө))+ Sinh(v ө) 


a-v(r-0) b=xrv 
Sinh(v (2x — 0)) + Sinh(v 0)= Sinh(b + a)+ Sinh(b — a) = 2Sinh(b)Cosh(a) 
= 2Sinh(x v )Cosh(v(x —0)) 


T,(r,0,z)= ы Соз(и z)fa PE Оаа) K, (u r)K, (un) 


m A Cosh(z v) 
_ Cosh(v(2z —0 — 0,))+ Cosh(v(0 — 0, )) 
Cosh(v (т — 0))Cosh(v (x - 09, )) = 2 


| Ж Cosh(v(2z —0 —0,)) + Cosh(v (Ө —0,)) 
Á Cosh(z v) | 
A l'aide de la formule suivante dans F.Oberhettinger, < Tables of Fourier Transforms and Fourier 


Transforms of Distributions >, édition 1990, section 1.17 Modified Bessel Functions of Arguments x 


х? and 1/x, formule 17.23, page 90) : 


+оо > 5 К Я 
[аи Cos(u z)K, (u ro )K, (u r) = T Р, Ё +r ÈZ | 
0 


2 +œ 
=> T,(r,0,z)= + [аи Cos(u z) K, (и r)K,, (и n) 
0 


4Cosh(x v). йу. ZT 2nr 
I| vient 
2q t , Cosh(v(2z —0 —0,))+ Cosh(v(0 —0,)) t 
"dech [ау = vez B = (Á Df d Cos(u z)K, (u r)K, (u r.) 
0 0 


+œ М 5 3 
= T,(r,0,z) ROM GA fav Cosh(v(2z — 0 -0,)* Cosh(v(0 — 0, dp, | 5n +r? +z 
2 % Cosh? (z v) Jy ps 


2 
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La formule suivante dans « A.P.Prudnikov, Yu.A.Brychkov, O.l.Marichev, Integrals and Series - 
Volume 3 - More Special Functions» en page 181 formule 2.17.24.6 (attention erreur corrigée dans 
le présent document), donne dans notre cas : 


[ar Cosh( E à GE 1 122 Avetün 1+ Cos(b) 

Соѕ (z<) uem J2 Je — Cos(b) c— Cos(b) 

2 2 2 
vv b—2z-0-0, c TÉ Cos(2z —0 —0,) - Cos(0 +0,) 
ro r 

[av Cosh(v(2z -0-0 0)» n +1? + z? Е 
° Cosh? (л v) SC 2nr 
_ => Р m pro" Р m SAC +0,) 

2 D ке 0990.) à fo +” TZ _Cos(0+0,) 

2r,r 2nr 


2" 2 areta Bn | 2.3 БОШ | 


== 1— 
= gé n tr +2? 24 rCos(0 +06,) 


NHF 2 a Bc | ë 1+ Cos(0 +0,) ) 


1— 
А +r’ +z? – 25 — 2 ro +r? +z -2r, rCos(9+06,) 


2 2 2 
E. 0 ÈZ 
ту b—0-0, c 780 17 2 


2nr 
Tav Cosh(v(6 —Ө,))ь n ote um 
A Cosh’ (z v) -siv 2nr 
De — ` is 2 Here ; pu сөзе 
мү у ein ^ To +” tZ Cos(0—0,) 
2r,r 2nr 


2nr N^ les 1+ Cos(0 —0,) 
= SS ү q J 


Jan +r? + z? =2r, rCos n +° + z? -2nrCos(0 -0,) 


Jor É ac qe | 1+ Cos(0 —0,) ) 
(9 — 6,) de Ñ j 


= 2 2 
Jn? +r? + 22 An rCos n +r? *z'-2nrCos(0 -0,) 
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Et /а solution devient, aprës ces divers calculs : 


+ 2 2 2 
= T,(r,0,z)=- 4 [dv — аш у | h tr +z 
2 % Cosh (т v) QU 


_ q T7 , Созй(у(Ө -0,)) p r +1? + z? 
2Jnr 5 Соѕћ (z v) Eu 2nr 


= T.(r,0,z)=- 


2 | DES | I+ Cos(0 +0,) ) 


1— 
Se tr? + z? -2r rCos(0 + 0,) | m n +r? *z'-2nrCos(0 +0,) 
— 2 Б aar | 5 . 8) 
24 n +r? ez? —2r,rCos(0 — 6,) dd n tr^ +22 -2nrCos(0 —0,) 
Le résultat final s'exprime ainsi : 


T(r,0,z)= 4 +T,(r,0,z) 
үг? tn -2rrCos(0 — 0, )+ z? 
T(r,0,z)= 4 4 


2/r° +” —2r rCos(0 –0,)+ 2? Ë 24 n? tr +2? – 25 rCos(0 +0,) 


+ d Arc (2%, gd 1+ Cos(0 + 0,) | 


zn +r? + z? -2n rCos(0 + 0, n +r? *z'-2nrCos(0 +6,) 


2 жаш, | 2r. 2 1+ Cos(9 -0,) | 


+ 
n [ro tr? +z? -2r, rCos(0 -0,) n +r? + z? -2nrCos(0 —0,) 


En posant 0,-n, on retrouve bien la formule de l'exemple précédent : 


T(r,0,z)= q +T (r,0,z) 
Jr? +” -2rr,Cos(0 -0,)+ z? i 


0, =z => Cos(0 —z)= Cos(0 +z)= —Cos(0) 
q Е q 
24r? +r, +2r r,Cos(9)+ z? dre tr +2? + 2р r,Cos(0) 


T(r,0,z)= 


+ 4 Arc zr | = TECNO) | 
nn er + 22 +2rr,Cos(0) n +r? + z? +2rr,Cos(0) 


+ q areta Tor | ; EZ | 
nn er + 22 +2rr,Cos(0) n cr +2? -2rnCos(0) 


= T(r,0,z) = “4 Arc J37 a ur Cosi) | 


2 2 2 2 
zr, +r? +z’ +2rr,Cos(0 n tr + z? +2r r,Cos(0) 
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On peut écrire aussi la solution sous cette forme : 


T(r,0,z)= zi s and | 1+ Cos(8 -0,) |, 


т\т? An —2rr,Cos(0 — Ө, ПЕ en rCos(0 —0,) 


q Ë arera [Bu | 1+ Cos(0  6,) ) 


Д тү! +r +2? -2rnCos(0 + 0, ) n +r? +2? -2r,rCos(0 +0,) 


N 


л n +1? +7 
— et posons Cosh(a)- 9—————— 
2 2nr 


T(r,0,z)= 2 Ё клн 1+ Cos(9-0,) | 
т\т? +r, -2r r,Cos(9 -0,)+ z? (2 Cosh(a )- Cos(0 . 


_ q Ë _ arctan | 1+ Cos(0 + 0,) | 
т\л? +r? +2? -2rnCos(0 +0,)\ 2 Cosh(a )- Cos(0 +, ) 
ou bien 
T(r,0,z)= 4 Æ Al) 1+ Cos(8 — 6.) | + 
т\т? tn -2rrCos(0 —0,)+ z? 4 2 Cosh(a)-Cos(9 – Ө) 


4 "T arco | Cosh(a)=Cos(0 + 2) 


1 
Pour x > 0 Arctan(x)+ таал L.) = 
x 


т\л er + 22 -2rnCos 1+ Cos(0 +0,) 
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On уа tenter de retrouver la formule donnée dans l'ouvrage N.N.Lebedev, I.P.Skalskaya et 
Y.S.Ufliand, 1965, < Problems of Mathematical Physics », exercice 422, page 198. Pour cela 
introduisons les expressions suivantes : 


1 2 
Sachant que pour x > 0 Arctan(x)+ Aretan +) = 2 еї 2 Arctan(x) = атац X 7 | 
x =x 


SR — cod е = ) 2 Cosh( 2) — со 2) SR — cos 
x, = 0+0 zd 22 7512 0+0 Ce 0+0 
Cosi 2) + со +% | l-x, SR + co | TN SR + соз 0 | 
2 2 2 2 2 2 


D 
SE 
D 
S 
k 2 d 


2 
2 = 
bm Cosi Z ) + co 22% )- Со) со + d 
2 2 2 2 


(90) + Cos? EEN TE J — Cos ° T | 
ER M 2 2 2 2 
Term c (2*6 
2 
0 +0, a 0+0 
Cosh? Ze d d Cosh?*| €. |- Cos? о 
2х, E ji E Е 2 В ut. d 2 )  [|Cosh(a)- Cos(0 4-0, ) 
1— x° = ) m (2+4, 1+ Cos(0 + 0,) 
2 2 
Cosh Œ бөз] 
2 2x, Cosh(a)- Cos(0 = 0,) 
= а Ec 1+ Cos(0 — 0,) 
ACTE + ca ` x 0 


= Arctan Cosh(e )— Cos(@ + 6,) = 2Arctan(x,) et Arctan Cosh(e)— Cosi — 6.) = 2 Arctan(x,) 
1+ Cos(0 + 0,) 1+ Cos(0 — 0,) 


даа | 1+ Соз(Ө +0,) | E : i arcus | Cosh(a)— Cos(0 + 0,) 


ЕР Кш 2Arctan(x, ) 
Cosh(a )— Cos(0 + 0, ) 1+ Cos(0 + 0, ) 2 
= Arctan Ey Соз(Ө na Ө) =Z —2| Z _ Аксап ANE 2 Arctan Adm 
Cosh(a )— Cos(0 + 0,) 2 2 X5 Xo `Z 
Arctan E Соз(0 = 2 = 2 Arctan Е 
Соѕћ(о)– Cos(0 — 0, ) x 2 
=> 


Arctan ы Соз(Ө + 6.) = 2 Arctan AME: 
Cosh(æ)-Cos(0 + 0,) x, J 2 
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Posons : 


Cosh| “ |+ Cos бш, Cosh| €. |+ Соз еб, 
2 2 1 И 2 2 1 
+ 


= YW, = 
SR — Cos -0o | ži Со) Е Cos ° б, | CS 
2 2. 2 2 


Arctan 1+ Cos(9-0,) = 2 Arctan(w,)- 24 
Cosh(a )– Cos(0 — 0, ) 2 


1+ Cos(0 + 0) л 
A 2 =2A SS 
elek )- Cos(0 + 0, ) J кә) 2 


T(r,0,z)= L Ë + 2 Arctan(w )-=)+ 
т\т? tn —2r rCos(0 —0,)+ z? \2 E 


= q Ё - 2 Arctan(w )+ z) 
т\л? tr? +z? -2rrnCos(0 +0,) 2 E 


T(r,0,z) - 2gArctan(y,) E 24 Ё — Arctan(y)) 
т\т? + n -2r ryCos(0 –0,)+ 2° zn tr +2 -2r ryCos(0 + 0,) 2 


манат 1.) 
T(r 9 pa 2gArctan(y,) _ V; 
т\т? + n -2r ryCos(0 —0,)+ 2° т\л? +r7*+z*-2r ryCos(0 + 0,) 


aret / 
=> T(r,0,z)- 2q Arctan(v, ) _ V, 
Z ү + -2rnCos(0 —0,)+ 27 Jn e +2° —2r rCos(0 + 0, ) 


En quelque sorte, on retrouve l'expression dans N.N.Lebedev, I.P.Skalskaya et Y.S.Ufliand, 1965, 
« Problems of Mathematical Physics », exercice 422, page 198, mais il semble y avoir une petite 
erreur de signe dans la solution de l'exercice qui s'avére étre la suivante, si nos calculs sont justes : 


SR + Cos ° = | SR — od 2+6) 
2 2 _ 2 2 


= YW, = 
SR — Со 00) Со) + ca 242 | 
2 2 2 2 


=> T(r,0,z)= A Arctan(y,) Е Arctan(y., ) 
zi NE An — 2r r,Cos(8 —0,)+ z? Jr ir pz — 2r nCos(8 "my 
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Exemple : potentiel d'une charge q, placée extérieurement à un secteur d'angle 28, conducteur à 
une distance a de l'origine. 


Le secteur d'angle conducteur parfait est orienté avec le bord coïncidant avec l'axe z, il est centré 
sur l'axe x, avec un angle 20. On suppose que la charge est dans le plan <x,o,z>, placé 
extérieurement à la distance a de l'origine. Soit en coordonnées cartésiennes la charge q est placée 
en (-a,0,0). 


Y < 


Le potentiel créé par la charge q est le suivant : 


q q 
T (r,0,z) = = 
` J(x+a) + y? +? x? +а? +2ax+ y? + z? 


x = r Cos(0) T (r.8,2) q 
=> T, (r,0,z) = 
r = Jx + y? li Jr? + à? -2ar Cos(0)+ z? 


La solution est recherchée sous la forme d'une addition du potentiel crée par la charge q et d'un 
potentiel recherché : 


Tëlee O e, т = m =+T,Gr,0,2) 
r a ar Cos z 


T(r,8,2) o-a, =0 T(r,0,z) =0 


Ө=2л-9, 
9 > 9, 9, > 2л – 9, 
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Le potentiel T(r,0,z) répond alors aux conditions aux limites suivantes : 


T,(r,0,z), „ =T„(rs0,2), = f°(r,z)=- q 
d z š [к WW r? + a? +2ar Cos(9,)+ z? 


9 > 9, 9, > 2л — 9, | 

L'expression de la transformée de Fourier-Cosinus de la fonction limite est relativement identique 
en la forme à celle de l'exemple précédent (substitution го->а, 0o-» rt-0;), il vient donc : 

gr. и) = - ки Е +a? +2а rCos(8,)) e^ (0, и) = -22 ки (и а) 

Et la fonction modifiée est donc la suivante : 


g,r.u)- #4 к PA +а? +2a rCos(9,))+ e "" KY (u a) 


Le coefficient G se calcule donc de manière similaire à l'aide de la formule : 


Comme 
| Ж KY (Ja? +t? + 2atCos(9))K,, (t) _ z K, (a)Cosh(v Ө) 
à t v Sinh(z v) 


Ce qui donne : 


| 4 
С (v, u)= = [кн a)Cosh(v 0,)— KY (и а) 

. 2 | 4 
G (v, u)= буу, u)* — 8,00, и)= Gv u)- ei (и a) 


=G,(v,u)= E (u a)Cosh(v Ө,) 
л 
Се qui donne еп сһетіп vers Іа solution du problëme : 
99,  9,o2z2—9, > 9,9, -2zn-29, 9-9 — 9-0 9,-9-=>27-9 -9 
4q T Nn K, (u a)Cosh(v 0,) 
T r,0,z)2—— | du C. dv K. iv 
aler z) Ee | H os(u gd v „(и r) Sinh(v(2z —29,)) 
a-v(r-9,) b=v(z-0) Sinh(a+b)+ Sinh(a — b)= 2Sinh(a)Cosh(b) 
2Cosh(v 0, )Sinh(v(z — 9,))Cosh(v(z — 0)) 
Sinh(2v(x — 9, )) 
Cosh(v(z — 0))Cosh(v Ө,) 
K. K. 
Cosh(v(z — 9, )) „(и r) „(и а) 


Е Cosh(v(z -0+0, ))+ Cosh(v(z —-0-0, y 
Cosh(v (т — 0))Cosh(v 0, ) = 2 


2q f "^t = Cosh(v(z —0 4-0, )) - Cosh(v(z -0-0 
Feel [au Cos(u z) [av Di Ge o) 


0 


[Sinh(v (27 — 9, — Ө))+ Sinh(v (9 -0,))] 


K,, (u r)K, (ua) 


tiré [аи Cos(u z) [av 
[U 0 


K, (u r)K, (ua) 
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A l'aide de la formule suivante dans F.Oberhettinger, < Tables of Fourier Transforms and Fourier 
Transforms of Distributions >, édition 1990, section 1.17 Modified Bessel Functions of Arguments x, 
x? and 1/x, formule 17.23, page 90): 


[аи Cos(u z)K,, (un )K;,(u r) = 
0 


2 
л? i pe 


P 
4Cosh(x v) Ir r a 2nr 


Il vient 
2q Ф, Cosh(v(z —0 +0,))+ Cosh(v(z -0 —0,))' t 
Jr zl d du С K, K, 
(©, ,z) 2 | v л | H os(u z) „(и r) „(и а) 
пв) _ A fay SRE Pre Corte dro rs 
MEETS 25 r Cosh(z v )Cosh(v (т - 9, )) -——iv 2ar 


Et là on semble stopper dans l'élan des formules miracles, aucune n'est présente sous le boisseau. 


Ce qui ne nous empéche pas de donner la solution formelle du probléme : 


= Tr, 0, z) 1 m 
"ue +a? +2ar Cos(0)+ z? 


Cosh(v(z -0 +0,))+ Cosh(v(z —Ө -0,)) „ (= 


= =Í [е Cosh(z у )Cosh(v(z — 9, )) -Lriv 


Cas limite de l'angle 90-n/2 


Dans ce cas l'expression prend la forme 


=> T(r.0, z) zi = 
dne + a? +2ar Cos(0)4 z? 


q [av 2Cosh(v(z — 0))Cosh(v ө.) е +7? + = 


2Jar % Cosh(z v)Cosh(v 9,) -Lriv 2ar 
Ke 2 2 2 
=> T(r,0,z)= q ке о у=; Е + +2 | 
r? + а? +2ar Cos(0)+ z? Jar A Cosh(z v) -——iv 2ar 


En consultant l'ouvrage en russe < A.P.Prudnikov, Yu.A.Brychkov, O.I.Marichev, Integrals and Series 

- Volume 3 - More Special Functions» en page 182 formule 2.17.25.3, il est donné la formule 

suivante : 
(2.17.25.3) | ах 


u 


Cos (b x) (eis 5 #3 и al | [1+Cosh(b 
Cosh(x x yos 


(с+1)2(с+ Cosn(b)) 2 Р" c+ Cosh(b) 


On peut transformer cette formule en par le changement b-»ib, et donner ses deux déclinaisons 
pour la valeur u=0 : 


Cosh(b x) BER u PHA 1+ Cos(b 
2.17.25.3—bi dx == < =2 ? 1 Cos(b E ui e e 
| is) Í Cosh(z x) к. CAGE | c+ m 


[а Cos( (b x) sje 1 
u-0 teil Cosh(z x) = J2 [c + Cosh(b) 
P 0 
i 1 
E 


Ж s J2. [c + Cos(b) 
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L'application de la derniëre formule à notre cas donne immédiatement : 


+оо 2 2 2 
paces (c) SÉ туч ос Бшш 


1 
Cosh(z x) не SCH X2. c + Cos(b) K 2ar 


0 


T(r,0,z)= 4 _—@ Ta Corin cB (#2) 
Jr +a? + 2ar Cos(9)+ z? Jar А Cosh(z v) ri 2ar 
= 4 __4 1 
Jr? + а? +2ar Cos(0)+ z? Jar Va [ET tZ +соцт-ө) 
2ar 
=> T(r,0,z)- 4 4 


i J +a? +2ar Cos(0)+ z? Д Ма? tr? +z? -2arCos(0) 


Ce qui est exactement la solution trouvée par la méthode des images électriques, soit le champ 
d'un dipóle de charges inverses placées symétriquement au plan conducteur parfait. 
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Exemple : potentiel d'une charge q, placée intérieurement à un secteur d'angle 29, conducteur à 
une distance a de l'origine. 


Le secteur d'angle conducteur parfait est orienté avec le bord coïncidant avec l'axe z, il est centré 
sur l'axe x, avec un angle 29, On suppose que la charge est dans le plan «x,o,z», placé 
intérieurement à la distance a de l'origine. Soit en coordonnées cartésiennes la charge q est placée 
en (a,0,0). 


Y < 


Le potentiel créé par la charge q est le suivant : 


q q 
T,(r,0,z)= = 
` Je af ie x? +а? -2ax y? + z? 


x =r Cos(0) T (58,2) q 
=> r,0,Z) = 
r= x° + y? i Jr? +a? -2ar Соз(Ө)+ z? 


La solution est recherchée sous la forme d'une addition du potentiel crée par la charge q et d'un 
potentiel recherché : 


ш лге 2 = к = + T,(r,0,z) 
r а — zar Cos VA 


T(r,0,2)|,__, =0 E D =0 
9, > -9, 9, > 9, 
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Le potentiel Т.(г,0,2) répond alors aux conditions aux limites suivantes : 


T,(r,0,z), o =r, = f° (r,z)= — 4 
j lo ` z A г? + a? — 2ar Cos(9,)-- 2? 


9, => -9, 9, > 9, 
L'expression de la transformée de Fourier-Cosinus de la fonction limite est la suivante : 


Z 2 
ele, u) = =K ME * a? 2a rCos(8,)) g2(0, u) = SPERI (u a) 
Et la fonction modifiée est donc la suivante : 


ele, ul = а KZ PA +а? —2a „Соз(Ө,))+ e ""KY (u а)) 


Le coefficient G se calcule donc de manière similaire à l'aide de la formule : 


Comme 

P (Jaz + -2arcos(9)k,. © _ z K„(a)Cosh(v (z — 0) 
À t v Sinh(z v) 

Ce qui donne : 


G;(.u)- - ^4 [K (и a)Cosh(v (x -0,))-K3" (u а) 
G (v. и)= Gv. u)+É 8? (0.u)= Gv. u) Zei (u a) 


4 

=G (ии) = K, (u a)Cosh(v (z –0,)) 

Ce qui donne en chemin vers la solution du problème : 

9 —-9, 9,59, > 9,—9,-29, 9-9 ->9+0, 9,-9—0,—-9 

1, (6.2) - - 4 Tau Cos(u z) [dv K,,(u A E o. H 
qm ^ Sinh(2v9,) 

a=v9, b-v0 Sinh(a+b)+ Sinh(a — b) = 2Sinh(a)Cosh(b) 

2Cosh(v(x —0, ))Sinh(v 9, )Cosh(v Ө) 


K, (u ғ), (ua) 


Sinh(2v9,) 
_ 44 ф TC, Cosh(v(z -0,))Cosh(v Ө) 
=> T(r,9,z)=_-5 [au Cos(u aja Gate) K, (u r)K, (ua) 
2q `° + Cosh(v(z —0, +0))+ Cosh(v(z —0, - 0 
=> 10.6.2) -75 [du Cos(u z) [av = соз WS ) x, (a r)K, (ua) 
0 0 0 
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A l'aide de la formule suivante dans F.Oberhettinger, < Tables of Fourier Transforms and Fourier 
Transforms of Distributions >, édition 1990, section 1.17 Modified Bessel Functions of Arguments x, 
x? and 1/x, formule 17.23, page 90): 

[аи Cos(u z)K,, (u a)K,„ (u r)= 
0 


2 
л 


P É tr +z | 
ACosh(z у) Jar -yuv 2ar | 


Il vient 
24 Ф, Cosh(v(r-0,+0))+Cosh(v(r — 0, — 0)) 7 
(r.0,z)=—— 12 0 GER 0 [au Cos(u z)K,,(u r)K,,(u а) 
ST (r 0 z)-- [av Cosh(v vír- 8, + Ө))+ Cosh(v(z — Ө, -0)), a +r? + z? 
н Seel Cosh(z v)Cosh(v Al -5+iv 2ar 


Et là on semble encore stoppé dans l'élan des formules miracles ! Ce qui ne nous empêche pas de 
donner la solution formelle du problëme : 


=> T(r,0,z) zi 2 
Je +a ? -2ar Cos(9)+ z i 


q S S tt tes) (= 


2 Jar Ó Cosh(z v )Cosh(v 9,) 


2ar 


Cas limite de l'angle 90=x/2 


Dans ce cas l'expression prend la forme 


= Tr, 0, z) 7 = 
Je +a ? -2ar Cos(9)+ z i 


+оо 2 2 2 
q fav Ee É +r +Z J 


Ser Cosh(z v)Cosh(v 9,) 2ar 
=> T(r,0,z)= 4 jav Cote), (= 
J +a? —2ar Cos(9)+ z? ; 5 Cosh(z v) -54v DR 
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En consultant l'ouvrage en russe < A.P.Prudnikov, Yu.A.Brychkov, O.I.Marichev, Integrals and Series 


- Volume 3 - More Special Functions> en page 182 formule 2.17.25.3, il est donné la formule 
suivante : 


"t, Cos(bx) BER и Nu 1+ Cosh(b 
2.17.25.3 dx = > P" =2 2 1 Cosh(b pu] тте N S 
) | ` Cosh(z x) 4,0) GES QE d gd 


On peut transformer cette formule en par le changement b->ib, et donner ses deux déclinaisons 
pour la valeur u=0 : 


Cosh\b x) (b x) -u-l H Nu 1+Cos(b 
2.17.25.3— bi. =2 ? 1 Cos(b PSI E Een 
( is) IE Cosh(z x) (z x) Guck (+ Blies os( ) É d E 

[ах Cos(b x) yE 1 
я e Cosh(z x) => iv \/2../с + Cosh(b) 
и=0 Р?(х)=1= 
Га Cosh(b х) (е 1 
А Cosh(z х) -Lriv X2. Mc + Cos(b) 


L'application de la dernière formule à notre cas donne immédiatement : 


n Cosh(b x) E T 
DUI K b—0 

| ` Cosh(z х) -Lriv c)= ader cale) x > V > с => а 
T(r,0,z) q _ 4 [av Cosh(v 0) , É +r? = 

ar +а — 2 Jar ` Cosh(x v) -2+iv 2ar 

_ q x 1 

r? +a? -2ar Соз(Ө)+ z? SC BE 
ar 


= T(r,0, z) d x si 
gU + a? —2ar Cos(0)+ z? Ма? +r? +z2 «2arCos(0) 


Ce qui est exactement la solution trouvée par la méthode des images électriques, soit le champ 
d'un dipóle de charges inverses placées symétriquement au plan conducteur parfait. 
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Problème aux limites de Dirichlet inhomogëne angulaire sur un domaine bornée еп z et non 
bornée en r 


Etudions maintenant le problème suivant bornée en z, avec des conditions aux limites homogènes 
en z. 


CL. T(rQ,z),, = 30,2) Т(,0,2),, = f26.z) 
T(r,O,z)|, , = T(r,0, z), = 
Cette fois-ci les fonctions axiales présentent des valeurs propres discrètes répondant aux conditions 


homogènes axiales. Il convient donc de rechercher la solution sous forme d'un mixte de série et de 
représentation intégrale, comme suit : 


T(r,0,z)= 2 Sin(u, z)[av K, (u, ^J Gea Sinh(v (9, -0)) + С, (v, u ) Sinh(v (9 -0, Д 


Sinh(v (0, — 0, )) "^ Sinh(v (0, —6,)) 


А пл 
ой Uu, =" 


І 


2 


Chaque fonction limites se développe еп série de JA propres axiales 


Ío ? (r, z = eil r All ZZ je el r SN -jagi (r, 2 BE d 


n-l 


Il est alors aisée de calculer l'expression des mM С, comme suit : 


5 +оо K. 
K, (v) = [ар g(e)K,,(u p) transformée de Kontorovich — Lebedev 
p 
0 


K = (r) = À [4 vv Sinh(x v)g(v )K,, (ur) transformée inverse de Kontorovich — Lebedev 
л 0 
g,(r,u,)= | G,(v.u,)K, (u r) pour p -12 


=> е (г. ge Z fav v Sinh(x v)g (v ei (и, г) = [dv G (v.u,)K,, (u, r) 
0 


gp. u,)K,, (и, p) 


= G,(v.u,)- к. Sinh(x (у )= EL Sinh(x v) Í dp 
0 


La solution formelle du probléme est donc : 
C.L. ps trai —TG8z),, = 70.2) 


gr, и,) EV f (r J [= :] ой ц, = E p=] 
L 0 


z 2 


gp. nu, )K,,. Qu, р) 


G,(v.u,)= Lv Sinh(x "ae 


T(r.6,2)- SI 68|] Let Дев) Sinh(v (0,-0)) (р) =] 


= А Sinh(v (0, — 0,)) Sinh(v (0, — 0,)) 
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Il arrive également que les conditions de finitude de l'intégrale garantissant l'application du 
théoréme de Kontorovich-Lebedev ne soit pas remplie, et dans ce cas comme précédemment, on 
procède à une régularisation des coefficients de développement des fonctions limites. En réalisant 
les mémes types de calculs, la solution formelle s'écrit maintenant : 

CL T(6,2),-4' 02 Т0,0,2),, = 0,2) 


T(r,0,z)| _, =T(r,0,z)| , =0 


L L 
gr (r, и, ) = 2 [а £ (r, LE z] = go (0, и, ) = 2 [а d (0, as 97 d 


z 0 2 z 0 2 


L 
gran) Í dz EECH "2 d où H, = p=12 
z 0 


1 

; * 
tel que [а [7 (r. Zn 
` r 


2 


1 
Los E р=1,2 
r 


g,(p.u,)K,,(u, р) 


G'(v,u,)= Lv Sinh(x v) [dp 


0 


* 2 
et G (v,u,)=G' (v, u,)+ “g°(0,u,) 
Ж 


"E [ял үе Sinh(v (Ө, — 0)) Sinh(v (0 —0,)) 
T(r,0,z)= 2 ZÉ d es K; (u, LEE (0, 0 )) + Gv Hs) shy (Ө, 2) 


Z 
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Exemple : répartition de la température dans un secteur d'angle de rayon infini d'épaisseur lz, et 
d'angle 80, dont toutes les faces sont portées à la température nulle, tandis que la frontière 8- 


90 est à la température /(,) SE z] 


Le problème aux limites est donc le suivant, avec le développement de la solution : 


CL. T(rj&,z) ,-0 Т(,0,2),, = QUES d 
T(r,0,z), , =T(r,0,z)| _, =0 Ө,=0 0,=0, 
f (r.z2)=0= gi (r,u,)= 0 


Р (r,z)= ros" d f2(0, z) = pos z] 


2 


gilr, и,)= 2 | dz s Dum IR er GK d où ци, = < 


z L L 
= g;(r,u,)= Ze- gr 70) dz zi =< гр" z] 
gru )- 2 (/()-e *'/(0)р,„ =(1(r)=e*" O, 


G;(v,u,)= Ev Sinh(x v) Í do (/(o)- e" "ro. (u, p) 


* 


2 
G,(v,u,)= G;(v,u,)+ = 82 (0, t, ) 


| : +œ Sinh(v Ө 
T(r,0,z)= Sin ZZ =| e K, (u, r )G (v, и, nee) 


Continuons les développements, il vient : 
Lie)" 700) К, (и, P) 
р 


ch, CN )= 2 Sinh(z v) [de 
л 0 


G,(v.u, )- Giu, )«7 (0) 


K, (u, ro — 2 (f(e)-e* ге iv (а, 2, 
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Exemple : répartition de la température dans un secteur d'angle de rayon infini d'épaisseur lz, et 
d'angle 80, dont toutes les faces sont portées à la température nulle, tandis que la frontière 8- 
90 est à la température f(r,z) 


Ce sont les mêmes développements que l'exemple précédent, à ceci près que toutes les 
harmoniques axiales sont utilisées : 
0, = 0 0, = 0, 


C.L. T(r,0,2| =O TGr,0,2),., = 052) 


T(r,0,2)|., = 
2 À 2 À 
g5 ° (r, u,) len r ЕСЕ | >g? (0, u,) T. =; jes f(0, 2 is SCH J 
в0и,)= 0) "ein où и, Tt 


Posons f, (r)= [& RUE z] Ful 0)= f dz /( fir ЕЕ z] 


+œ 


9 r 0 
G;(v.u,)- Ze Sinh(z v) | dp loune" < g2(0.u,)IK,, (u, p) 
0 


G;(v,u,)= E Sinh(x v) Í do (t) enert, (OK, (и, p) 


» 0 p 
=H. 
T(r,0,z) => KS OI SES E v Sinh(x v) Í do Ux ice n ole, (д, p) ( 
Sinh(v Ө) | | 


à Sinh(v 0, ) 


<| r) Sinh(v Ө) E 


l. jSinh(v Ө,) 


(0)+ v Sinh(z v) fap 


AUS 
л A p 
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Sil'on suppose que la fonction limite est indépendante du rayon, alors, il vient : 


E r) Sinh(v Ө) 


1 Sinh(v Ө,) 


2 


ml, m L d |+ v Sinh( he v) ар (1- e "r P )K K, (u, 2) 
0 p 
K. пл, Sinh(v 0) С 
п ка 1 ` ) Sinh(v 0,) 


[r eN бее” ten 


Comme nous avons vu dans les premiers exemples de transformations de Kontorovich-Lebedev 
n don — p ёп P +оо | pP 
que l'intégrale [4e ( < К Ай p) 2 [ар (1 е kk. (p) | 
p 
0 0 


prenait la valeur suivante : 


Í do (1-e"")K.(o) _ „ H RE Ja Lee a 


А р v Sinh(z v) T A p 
II vient : 


4 Wem пл Sinh(v Ө) TV 
T(r,0,z)= ne Sin z 27 DE | ; r] Sinh(v 2 Cost á ) 


n=. = 


On peut encore simplifier ce calcul, en utilisant la représentation intégrale suivante ( voir l'ouvrage 
A.Erdelyi, H.Bateman, Higher Transcendental Functions, volume 2, page 82, section < Modified 
Hankel functions >, formule (25) ) : 


K, (x Е 2 -J- [а Cos(x Sinh(t ))Cosh(v t) si x>O et -1< Re(v)<1 


v Эіу = K,, (x)= ВЕЕ: ЗЕЕ. E Cos(x Sinh(t ))Со$ (у t) 


ol) 


= K,,(u,x)= =L ja Cos(u,x Sinh(t))Cos(v t) 
Cosh ЛУ k 
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Il vient : 
yS Cosi ЛУ 2 Ө) " Со "У 2 Ө) 
dv K, (u,r)= | dv dt Cos(u,r Sinh(t))Cos(v t) 
| Sinh(v 0,) | Sinh(v 0,) SE 3E | 
= [ау < Зико); Se Cos(u„r Sinh(t ))Cos( v t) -[а Cos(u„r Sinh(t)) [dv Sinh(v ө)Соз(у ) 
d Sinh(v 0 ^ Sinh(v 0,) 
0 
mu . Si T2 | 
Сотте [ dv Зану Are t) " 0, 
9 KÉ SA SR + Ska 
9, 9, 
A Ө А 
Se Gol ЛУ 7 У | Sinh 0) ы sl [сеш Sinh(t)) 
d K, = |d ` 
= | "= Sinh(v 6,) Juck: |а 


9 сөң Zt) + cos 72) 
6, Ө, 
La solution du problème s'écrit alors : 


0з) LT СЭЙ m zi ze са", Sinh(t 2 


? SR + ek 
% 9, 
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Exemple : répartition de la température dans le quadrant 0<x<œ,0<y<œ,0<z<1z, toutes les faces 
sont portées à la température nulle, tandis que la x-0 est à la température TO 


Cette configuration géométrique correspond à un secteur d'angle droit. Il vient de la solution 
précédente : 

л 
0 = 0 0, = 2 


C.L. T(r,0,z2 sf X T(r.0,23)),_x =T, Zëss ffr, = 


g2(r. u.) ЕЕ kA 110 (C1) )o g (r,u,,)= 0 


AT, 
82 (7o Honni) = 82 (0, Honi) = t ` 


(2п+1)= 
2n lr 


l 


z 


2l T, e, 
= d ius (r) = dish (0) = бег OU Mini F 


z fer r] Sinh(v 0) 
iv L | лу 
29 | Sien 
8 n=+œ 1 n 
Ө,2)= 2 : d 
(r,0,z) x2 > (Ол +1) | i = (1 ce tam P)K (2n+1) 
ig v Sinh( (z v) [ар L 
0 p 
us e (nl | Sinh(v 0) ` 
ve ed Seck, a ч 1, Sin С 
(r,0,z) E @ F fav 2 
qp == n 


ИЕ? E v ж 9 fap (1- к 2 


Nous avons vu dans les premiers exemples de transformations de Kontorovich-Lebedev que 


l'intégrale Í do (i-e ^), (Mansi P) _ [ae (1-e )к,, (p) | 
0 0 p 


p 


prenait la valeur suivante : 


=T 
p v Sinh(x v) л p 


et) Соң Z )-1 role) ee") (2) A 
[ар iv 31+ [ар iv P) = SE 
0 

Le résultat final attendu est donc le suivant : 


RS : À 5 si Beck, fa Cosi ЛУ 5 `" | Lr | bast З 


= ati) а al 7” 
2 


Résolution des problëmes aux limites de Laplace en cylindriques 3D par les transformations de Hankel et les 
transformations de Kontorovich-Lebedev- p 212 


En utilisant la représentation intégrale suivante : 
1 


Cosi Z zr P 
2 


=> K;, (u,x) SZ 


K, (x) = [а Cos(x Sinh(t ))Cos(v t) 


ja Cos(u,x Sinh(t))Cos(v t) 


Il vient : 


| Са "У У ut v 0) e Cosh( ZY `" buy Ө) 
dv 


) Sinh e а — d SEN C zx 
= fav З Í dt Cos(u,r Sinh(t))Cos(v t) = Í dt Cos(u,r Sinh(r)) [av Simo NUE) 


0 КЕ 0 Sinh E 3j 
2 2 


je Cos(u,r Sinh(t))Cos(v t) 


DR fav Sinh(v 0)Cos(v t) ü Sin(2 0) 
£ ЕЗ Cosh(2t)+ Cos(20) 
» Ts Cosi "У `" тб, 2 (i À= fa Sin(20)Cos(u,r Sinh(t)) 
0 Sinh Z) iv Wn ^ Cosh(2 DÉI Cos(2 0) 


Le résultat peut donc s'écrire sous une forme plus simple : 


ed Beck, 

ST, "к L Ka Cos(u,r Sinh(t)) 
T(r.0.2)= eu 2 а ^ 

(n6,2)- s > (2n+1) e Cosh(2t)+ Cos(20) 


ue ed Cre z] Jm Co tie r sinh) 


z 


0 (2n+1) e? Cosh(2t)+ Cos(20) 
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fav Sinh(v Ө)Со s (v t) 


Démontrons la formule qui permet de calculer l'intégrale : 
Sinn ^ ) 


Il vient : 


(mp 
Sin] ZE 
[ar Sinh(pt)Cos(rt) _ x d | 


Ô Sinh(qt) ТЕЗ "E Tre Í a oct. a 


Sinh(pt Em гт) Sinh(pr) ren (irc) Sinh(( p +ir)t )+ Sinh((p —ir)c) 
| - far =; [ar 
d Sinh(qc) Sinh(qc) 


Sinh( (ат) 


Сотте [ar Sinh(at.) _ ч SE d b > |а 
d Sinh(bt) 2b 
T 


far seit Z run? PE fa Sp) sae 
Ja Son) я {nf test), (nasi) MEN s| 


7 Е : Neo (zz Sei E к, 


Cos(a)Cos(b)= —(Cos(a + b)+ Cos(a — b)) Sin(a)Cos(b)+ Cos(a)Sin(b) = Sin(a + b) 


= Е 
SBS 


Soit en réalisant les substitutions suivantes, le résultat attendu : 
3 T 
Ska 


[dr E q g 
0 їпй\дДТ q ЕЕ 
4 4 
rov p=0 r=t ER > fav Sinh(0 v )Cos(t v) _ Sin(2 0) 
278 зт Y Cosh(2 t)+ Cos(20) 
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Et de тёте: 


Ф Sinh(pr)Cos(rr) _ ed 
r 
[ат SCH = ` — 4 p e [0,4[ 
9 SS (ат 1 cosi ==) -cof 22) 
q q 
0 

D sm е! 

T>V p=0 r=t q в = [а Sinh(0v)Cos(t v) о, 
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Problëmes de Dirichlet inhomoqëne anqulaire, Dirichlet homogëne radial 


Le problëme de Dirichlet inhomogëne angulaire décrit précédemment était non borné en r et 
borné ou non-borné еп z. Introduisons une condition aux limites radiale homogène : 
AT(r,0,z)= 0 C.L. TO. = 3 (r,z) T(r,0,z),, = (т>) те, = 0 
La solution de base issue de la séparation des variables est cette fois Іа suivante : 
T(r,9,z)- (K, (и DAT L.)- K;, (u l, e (и r)XA,Cosh(v 0)+ B,Sinh(v H ESCH DÉI B.Sin(u z)) 
Pour une fonction ni paire ni impaire, on sépare le problème en deux sous-problëmes avec 
respectivement les fonctions limites suivantes : 

° (r, + f° Le MUS — f? (r.- 
fi (r. z)= Í ( z) L ( z) L (r, z)= Ío ( z) L ( z) 
On cherche la solution du problème sous la forme d'une double intégrale (sur le spectre radiale et 
angulaire) en introduisant un noyau intégrale qui s'annule en r=lr : 


L, (u r,ul,)- I„ (u r)K, Lui Li (p L)K,. (u r) 


T(r,0, z) = Í du Cos(u 2) Í dr L. (u r, ul, E (т, и) не: = +G, (z, m Ne EU. ) | 


L'application des conditions aux limites de Dirichlet donne : 


р=1,2 g9(r,u)= SL ? (r, z)Cos(u z) 


T(r,0,z),_ a = Í du Cos Hz fee nt Hr, ul.)G (r, u) = [au g? (r, u)Cos(u z) 
0 

T(r,0,z), T: = [ди соң uz GA ur, u, )G,(z, и)= [аи g5 (т, u)Cos(u z) 
0 


= g,(r,4)= [dr L.(ur.ul,)G, (c. u) 
0 


Dés lors la solution dépend des valeurs des fonctions G1 et G2 qui sont les coefficients de 
l'expansion intégrale par le noyau L d'ordre imaginaire. Comme nous exposons les résultats de 
l'article, voici la solution formelle du probléme : 
a & 
Stat 2 ENSE. g(ur)L (ur. ut) 
0 


F 


0 


2 T RE т) 
= — a0 ul L. 1 
e(r) 2 [ FR L. ( L) g(r, D „(ити » 


Comme 


Si A r ,u)= [ar G, T. u)L „(ur,ul,) pour p =1,2 


2 Pu Djar HE ul.) 


p 


Alors G, (z, u)= 


LA 
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Et la solution formelle du problëme est la suivante : 
L. (u r Hl, )- I. (u r)K, (ul, )- 1. (u LI. (и r) 


р=12 go(r,u)= 2 fa f?(r,z)Cos(u z) 


, Sin TT (ur. ul.) 
G ( a )=2 L;. 
Í a а 


T(r,0,z)= Í du Cos(u z) Í dr L. (u in 
0 0 


Sinh(z (0, — 0)) Sinh(z (0 — 0 )) 
Sinh(r(6, 0) с (ex) Sinh(z(0, — 0, d 


Cette expansion intégrale est possible дгйсе à un article qui étend la transformation intégrale de 
Kontorovich-Lebedev, dans une article de 1965 des auteurs Ya. S. Uflyand, E. A. Yushkova : < 
Solution of the Dirichlet problem for a finite wedge by means of special integral transforms using 
Bessel functions > 


Transformation intégrale de < Ya. S. Uflyand, E. A. Yushkova_> pour le problëme de Dirichlet 


Pour toute fonction f(x) sur l'intervalle [0,a), telle que H a)=0, si : 
1 - continue par morceau et dont la variation est bornée sur tout intervalle [x1,x2] de [0, a] et 


2 - telle que l'intégrale est finie : je rog[ LJ <œ 
x x 


Alors cette dernière admet la représentation intégrale suivante, pour x € [0,a] : 
L, BE . = 1. (x)K,, (a)-1;, GC D 


f(x) =< | ат т Sinh(z т) Lal dt /@) 
=Í EG] 


Sch a) 


En terme de transformée et ое inverse modifiée de Kontorovich-Lebedev, il vient : 


L, (x,a:) te a) 
f(x)== |атт Sinh(z т) T dt f(t) 
= | Deg er 


fG,a) = [а pe Pake) transformée 


=> 
f(x) = E3 far T Sinh(x г) оа) J,a) transformée inverse 
л 0 


Ila) 


L'article indique que lorsque la limite o tend vers l'infini, alors on retrouve la transformation de 
Kontorovitch-Lebedev : 


L. (х,а) = Lim iT 


| Lol @—› 


2: (a) = +00 > Lim 


а —00 


> sŒ) =Z; [de z simis) da ro 
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Pour un problème aux limites sur le secteur d'angle borné en z, alors la solution s'écrit : 
Ij Io lr) 10 (x) х= А, а = Al, 


Vno Endurance ct 2] 


п=1 


М пл 
ой use 


l 


z 


Chaque fonction limite se développe en série de fonctions propres axiales : 


So (rz => 2 F Se? = r, 4, ) -2[e se (r,z СЕ d 


п=1 


II est alors aisée de pai l | о, des coefficient S comme suit : 


)= ja £36 a) z(e, u) -far g(u,r Jeter. pa) 
Comme Si a D saed 
in (тт 
7 gír, l,l, )L. (ur, ul.) 
-— AURAI 
Si g,(r,4,) = [ат (r,u,)L.(a,r,u,l.) pour p =12 


0 


2 _ Sinh(z z) ' Í gp £o mL nu) 
: Tic (uL 0 

La solution formelle du probléme est donc : 

C.L. T(r,0,z)|,_, = fff (r. z) T(r,0,2)|,_,, = fy (r. z) 


Alors G, (z, ijs 


Л r 


T(r,0, cl, = o =T(r.0,2)|.., =0 
g(r.u,)- Zei me E: ой бе - p =1,2 


L, 1 
б.(с,н„)=-2т SHE E) pu, apr urs mt) 
Ii. (u,1,) 


Рае Sinh(r (0, — Ө)) Sinh(r (0 — 0,)) 
T(r.0,2)= X, Sin(u„z Í dt (шу, ui) с simne 0) 241) Sial 9, а] 


r 
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Problèmes de Dirichlet inhomogène angulaire, Robin homogène radial 


Dans l'article suivant daté de 1966 de E. A. Yushkova : « Certain singular boundary value problems 
for the Bessel equation and their applications in mathematical physics », l'auteur introduit une 
représentation intégrale qui donne la solution du problème de Dirichlet inhomogène sur un secteur 
d'angle de rayon limité, avec une condition de Robin homogène radiale. 


+BT(r,0,z) =0 


r=l, 


AT(r,0,2)=0 CL. T(r.0,2),.,=f?(r.2) Sech, m f (а) El 
1 =0, p 


La construction de la solution du probléme se raméne à la représentation intégrale d'un forme f(x) 
sur un intervalle bornée [0, a] répondant à une condition de Robin en x= a de la forme : 


Transformation intégrale de « E. A. Yushkova » pour le probléme mixte de Dirichlet et Robin 


Pour toute fonction f(x) sur l'intervalle [0,y), telle que f'(y) - hf(y) —0, si : 
1 - continue par morceau et dont la variation est bornée sur tout intervalle [x1,x2] de [0, a] et 
2 - telle que l'intégrale est finie : EEN < œ 

Ïx x 


Alors cette dernière admet la représentation intégrale suivante, pour x € [0,a] : 


н„(х,у)= K; Gor, (y )+ AL n) - Li GR, GP) hK, (0) 
Í dt rs E n 


joe [ат z Sinh(z 2 Hos y) 
cT à I 


1„'(у)+ nr, (yy 
пы, ны ONE ar) hL, (y) 


En terme de transformée et "at inverse modifiée de Kontorovich-Lebedev, il vient : 


а j H, (х,у) Y н.) 
f(x) = Li dt c Sinh(x t Е t (p) H, 
4 | L(y)eht. DI: d 


fG.y)= ja f pt) transformée 


F (x)= SE [ат T Sinh(z т) Н (x.7) 5 fGa,y) transformée inverse 
л 0 


1'(у)+ (у) 
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Lorsque l'on prend pour paramëtre h=6/a, soit : 


н (ко) Ke 9) Pru) Jr ко) Ê K.) 


р M a H, (х,у) ан.) y) 
Ек pu ren QI Ç Lat y: Sr dt f (p) нз?! 


Posons . 
Rs y)- e Н„(х,у)= K, exo 1, G) BT )- 1. (x)(@ К„'(у)+ B K) 
a f'(y)*Bf()-0 
= то R, (х,у) к) 
f(x) = =a [4 T Sinh(x Be (у )- Br, (у y: [a РО) ——— 
Dans le cas 6=0, a-1, condition homogène de Neumann, il vient : 
R,,(x,y)= Kx). Gr)- 1, G)K o) 
f'(y)-0 


т fe) = Z Jo E 


SÉ fp ETS R ke y) 
M 


Dans le cas 6-1, а=0, condition homogène de Dirichlet, il vient la solution décrite précédemment à 
savoir : 


f@)=O К„(х„у)=К„(х)„(у)—1„(х)К„ (у) 


rey= [ar r Simi 9) Gef t f Аф y) 


it | 


L'article indique que lorsque la limite y tend vers l'infini, alors on retrouve la transformation de 
Kontorovitch-Lebedev : 


LG) 
Lim P Gi = hk, (z) 
уэ» acht 
>f@=4 far t Sinh(z en " ; - [а = 
Ü іт ir Y Eh 0 it 
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Application de la transformation intégrale < E. A. Yushkova > aux problèmes aux limites sur un 
secteur d'angle de rayon limité et de hauteur limitée ou infinie : 


La solution formelle du problème illimité en z: 
" OT(r,0,z) 


CL бб, „ей ez ^ T(n82),,-402 — 


Ry (x, y) = Ka (xa LC) B 1.) - (xa Kv) вк, (у) 


g(r.u)- - ја у Cf, z)Cos(u z) у= ul, 


+BT(r,0,z) =0 


r=l, 


pour p =1,2 


2 


2 et, (v) Br. (y) 0 d 


T(r,0,z)= (du Cos(u zar К, (ur, DL (z, uu) sure zm £6; u) Sinh(z (0 — 0, )) | 


1, 0 
G (ss) 2 ч Sinh(x т) [а gir. ш), (ит, ul) 


La solution formelle du problëme limité еп z est donc : 
CL T(r,0,z2)|, = 700,2) Т0,0,2),, = 0,2)  TG,0,2j,,=TG.,0,z),., = 
BCE -0 
Or r=l, 
R,,(x,y)= K, «(x Xa 1„'(у)+ 81,„(у))—1„()(е К„'(у)+ B К„(у)) 


e r ri) Bes (r, 2 СЕ z] où L, ЕЛ. p=1,2 
0 


+BT(r,0,z) 


£ 


E Sinh(x г) 
@,(,и,)= si lat; (y )- Br, (у y 


T(r,0,z) 2 -È Sin(u,z ја R. (ur, DD (г, и.) (Ө, –0)) Sepe SHE а 0, E 


: R, l 
dr elen.) SEH » y = ul, 
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Quelques propriétés des fonctions de Lommel et fonctions de Lommel modifiées 


Fonctions de Lommel de première et deuxième espèce : 


Solution de l'équation différentielle de Bessel inhomogène 223"+z23'+(2? -y? y =z"! > | 


U+V+3 и-у+3. z2 
2 ` 2 ” 4 


ch 
s (z)= ap 


s, (2) 
(2) 


(u +1) -v? 
else än? DÉ . ц“ =) SE = 2 GL co н 22 2 
rz dax 
Ce J, (z) J, (z) 
Ѕіп(лу) [== dar 
2 2 
ati fat) 
24 2 - 2 |с — Ue) Cos z L dë 
Sin лу) 2 2 


e Suv (z)= Suy (z) i 


enemy (een) 


Liens avec les fonctions de Struve H et K : 


"oa Jaco ea e) 


K, (z)= H, ()- Y 


Deuxième Ordre v purement imaginaire, у=їт 


NUN {у Ven (уй, eer гу A 
ЕР {у HKO 


Si le premier ordre u est réel et l'argument réel alors la fonction de premier espèce est réel, ainsi 


que la fonction de 


deuxième espèce : 


Gär? 
2 


2 


8 ni (x) = u 


1+i l-i Co. 
Sui [z)= s. a (z) SE tir (z = 


2 


А ; 2 
ТА Аа аали 2 
bw E : Ë eR S 
(u+1 «e? ga 

r Sa (E 

T 2 

S, (z)= s )- 27 — J; J 

uil) suce) iSinh(zc) r ge ch {== 

2 2 


lant) d 


iSinh[zr) 


| СЕ 2 = © oos н = jz e) 
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Lorsque le premier ordre u est un entier, il vient : 


l=n 


Sor (о) pr Viel) Sona) sanse () CU" tomo IK (2141) 
Si (e) е) ROE) ` Stell suc) CY sr (of f> +a) 


Comportement asymptotique Z -> co 


Développement en série de puissance ascendante autour de 0 


2 
Re ШИ: VERVE уса | Йй 
ON кешн s 2 2 Jeya 
Ж (и+1) -v? (0+1)? =v? F ret dE d д2" 
I=+oo 121 I=+co Ï 2l m-l 
ablenne Eent YC) sinn! a(ua)=] nmi 2) 
1-0 ((u + 2m 1) E) £x alay) m-l 
т=1 
22 [u +1) eeh v+3 ) 
Et également а (и v)=1 a (и v)= : 
бм 22 (eu) (ex 
r Г 
2 2 
Développement еп série de puissance descendante autour de œ 
1=+о m-l 
S, (z)= z" DO) а (ру)? ag(- u,v)-21 a(-u,v)= (Qm -1- и) -v?) 
1=0 m-l 


l-u-v l-u-v 
ГИ) = =< +7 [D +1 
=À x 2 H 2 ү 

1 


Attention le développement еп série de puissance descendante est purement formel car exprimé 
sous forme d'une fonction hypergéométrique 3F0, il est clairement divergent, sauf si la série se 
termine, soit une forme polynomiale. Et pour cela, il convient que les valeurs (3-u+v)/2 ou (3-u-v)/2 
soit un entier négatif, ou encore que (u-v) ou (u+v) soit un entier impair strictement positif. Pour 
ces seuls cas on a le développement en série licite et convergent autour de = : 


Si u+v=2p+1 ou и-у=2р+1 


2 Ï = i-4- 4 
= 5,,(2)= z" E =, = Y 1; 2 
2 2 Z 


Développement des fonctions de Lommel de première espèce en fonction de Bessel 


Lia 1=+00 
j _ quel (1+и+1)Г@+и+1) р aen SS (1+ u +1) (+ u +1 
5 (0) Kee GEES ech) <= Sin -v?) 


) Jann (z) 


< (+1) 
1=+00 so (z)= 25 21 11V 8 Jami (z) 
ve s.j (z) 27 (21+ uer - ul) Wee (2141) e) 


ПЕТ nare 1 +1) +r?) leg 


А = 8 Ј 
Sii (z) Zr.) a(z) 
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Comportement asymptotigue uniforme de la fonction de Lommel de première et de deuxième espèce 
s lorsque v -> °° 


En appliquant la rëgle (с) Alen suy (2) 0, (Ci) (voir plus loin le comportement 


asymptotique uniforme des fonctions modifiées de Lommel), la construction est la suivante : 


— quu š 
8 „y (z) =-i t (-iz) A E _ à 
uH i235, G1| ЊУ, и+1 ко H,v,—i— 
u+ Z V z v 


e z = ul 2 Ду 
ul о У, = Suy (a)= b ( i) 2 > 21 MENU 2, 21 
2 v y^ = v y^ s v 


tV (z ` => 2l 
1=0 


Comme G, (u,v,6 lest uniquement fonction de £? G, La CL б, (uv?) 


Soit la construction suivante 
= d = d? ~ 
(ue UG (u.v,ç)+ (2 +3) ae Gen )+ É aci G, (u,v,¿) 


1-62 


G(u,v,6)= 


= Z 
HE 1240 G HV, 
2 4 


Suv (z) = ү? 21 
1=0 


isa- k- ei 


Avec Gv )=—- = Gi (u.v. c) (1 - e 


(1+ uy 462(7+2u u?) 26 "(175+34u 38u? би? +3џ*) Ach 39и+2и°+5и? vil &®%(и—1) (и —3)* 


G,(uv,¢)= 


0-2 


Réduit аи seul premier terme, l'expression est relativement simple, tant pour ип paramètre v réel 
ou purement imaginaire : 


zu+1 z"*1 


Su,v (Z) = — = Syair (Z) = ———— 
Š YË = ы S £ + z? 


Par ailleurs le terme additionnelle pour définir la fonction de Lommel de deuxième espèce а 
également une expression asymptotique < relativement > simple . 


s, ) s, (ef Er n EE 


+1+v +1-v 
e ill 277 (es We Ji (stent ette 


Sin( (ту) 


La fonction батта se développe comme suit : 


E Е 
T(z)e J2z z Pa)" gd id 
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Dans le cas d'un parametre v réel les fonctions de Bessel ont le développement suivant : 


Ce qui donne : 
1 1 
z V=> MW _ 2, 
S  (z)» rt dei ?y `e "gt 
2 у 
Et dans le cas d'un paramètre v purement imaginaire : 


sin[—] Cos |=] 
Ке [2 :[2]] - ———— Im[3i :[2]] 


_ 1+и+1т 2 
S, : 254 2⁄1 r| ———— | хг — 
wi e [7] saa ezntr[ 7 ы! Cosh[ =] Sinh [22] 


1+и-1т 
2 


Et les deux développements des parties réelle et imaginaire de la fonction de Bessel réduit аи 
premier terme sont les suivants : 


үк” + 22 7 i z 
£ [=s 2] = —— +log|[_ == ==, т] = Nt +Z + z Log] _ === | 
т 2, 52 2.2 
T+ T +z T+ T tZ 
лт 2 1 
е2 T 4 л 
Ке [75 :[2]] = É d Cos [z £z, z1- =] 
2хлт \T +Z 4 


e 


e? z? 4 x 
Im[J;.[Z]] = | | sin[z £[z, SM 
Ê 4 


ү2лт 


Ce qui donne une expression approchée du deuxième terme pour la fonction de Lommel de 
deuxiéme espéce : 


Ѕіп[ 22 ] Cos [22 | 42 т? 4 хи л 42 хи x 
- ^ > Re [Ji .[zZ]] = InDi [21] = | - =) Sin[=>-rt&[z, z] +>] = Sin[ — - « etz, с] +>] 
Cosh[**] sinh [22] Vat \t?+z 2 44 aat 2 4 
Ou encore : 
“été z] = z Log[ =] > A ` ч Ee -*] е мы ыш, d E t 
2t a pm 2+ 4 P = spa 2c 4 


ra etur] gari In[.[z2]1] _ V2 


хи ez л 
- = sin[ -z Log[ ]+ ] 
2 Cosh [22] 2 Ѕіп [22] afre 2 Zeil 4 


E: sin[ 


Comme d'autre part on peut exprimer une expression asymptotique du terme en fonctions Gamma, 
il vient finalement l'expression asymptotique réduite au premier terme : 


; S 1-и+1т (и-1) 
Comme: r[z]s 42x z" 2 e% et Log| ] = їх +21 
1-и-1т t 
II vient : 
1-u-ic at 1«u-ic rt 
|r[ => ] P 2 d 2 ST ев |r[=> > ] |? =2** xe 2 d 
2 2 
+1 лт 1 mu zz) 
> Su,i t [Z] = +V2x e 2 z^? Sin[=> -z etz, SM 
т?+7 2 т? + z? 
z"1 лт _1 mu л 741 
> Su,i t [Z] = A Wine 2 z^? sin[= -с1ов|—] +-| = 
vez 2 2 4 т2 +22 
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Fonctions modifiées de Lommel de première et deuxième espèce : 
Fonctions modifiées de Lommel à valeur complexe (voir article de 1936 P.Szymanski, < On the 


Integral Representations of the Lommel Functions ») et fonctions modifiées réguliéres en 0 


by (z) SS SN (iz) 


Tav (z) = AUS, (iz) 


Solution de l'équation différentielle de Bessel modifiée inhomogéne zy e? t v2)y Lg | 


.Htv 3 u-v43 z? 
( ) lu (i ) TE (i ) ( ) ns Д 2 ' 2 ` 3 ан z” ( ) r1 ( 2 1} 2) 
t ,(z)=-i "s ,(Gz)=e s, uz) t, ,(z)=z -z avec aj(u,v)- u +2m v 
S e xd: (+ уйл а 


—iT— ежа = 
PE "(2 Hir (e Me Du a ca [s 8" ute) 


= Pme Š I (z)=e 
T, (z) =t„„(z) prete et d 1,(2)+ 2 (zy 2 oett] = ; 


О a Te) (2 TEE (ete) 


Si les ordres u et v sont réels alors les fonctions régulières en 0 définies ainsi sont réelles. 

Liens avec les fonctions modifiées de Struve L et M : 

ы@)=-е THE) Yu) Asien) M) n )- 1,0) 

Е гі И Ch, вије гу я „=н, la Je, ^ js (2) 

E E r = air " ik (ўз И == o) cendo oct {у 2 323 (iz) 

E E Jor, e 3 («2 K, (нт) iCos(nv)} ele 2 k, (озии) = 1 , (2) 
elen 2 "elen s E) c eost) 

nu G2" ar v eO- 1. 0-2 eost) 

> T)ar v+ Ms 0-2K, Osina) ŽK, (саан) 


л 
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Il s'ensuit que les parties régulières en 0 des fonctions modifiées de Lommel sont proportionnelles 
aux fonctions modifiées de Struve : 


M, (z)- L,(z)- 1,(z) 
gemi) t, (z)- 277 Jar v+ (о) 
Н pu 


туша (т) =, (z)— т arv +s 3 ,G)= SCDE d 19 


Fonctions modifiées de Lommel à valeur réelle, (voir également article de 1964, C.N.Rollinger 
« Lommel functions with imaginary argument >) 


Développons la partie réelle des fonctions introduites précédemment : 


Solution de l'équation différentielle de Bessel modifiée inhomogéne 2? у"+2у' +(z? +v 2) =2" ә | 


Re Liu Ci es Ds Die >. 


Travaillons sur la partie réelle de T, (z) 


T, (z)=t,,() eere eme Bey Posez) 
— Sch 01212920) [рг = däi 2 Conf r E dite 


2 


+I+v +1-v 
T, (z)=t„,(z)+2*" i = ie БА v5 vd zs 1, (z)- Co Ey Te 


Sin(nv) 


— Кејт, , (z)}= bry (2)+2*"7 
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Comportement asymptotique de la partie réelle de la fonction modifiée de Lommel de deuxième 
espèce z ->œ 


LV UN 
ri 1-и+у l-u-v 4 
Suy (z) ~= z! E 2 , 2 BE £) 
a ei rus l-u+v l-u-v 4 1 l-u+v 1-u-v 4 
= Т, @@ (z) —i' SCH 12” SA 2 E 2 1; 2) = S E 2 5 2 WE 
I=+co т=1 
= Ref, (z)}= Tay (2) = Zu Y at: uk") ag(=u,v)=1 a(-u,v)= H»-i-4» _v2) 
1-0 т=1 


Comme l'usage de la fonction hypergéométrique 3F0 est à prendre avec des pincettes, il у а 
identification des développements à l'infini dans les deux cas cités auparavant : soit (u-v) soit (u+v) 
est un entier impair strictement positif, et l'on a le développement en série licite et convergent 
autour de = : 


Si u+v=2p+1 ои и-у=2р+1 


EN i-es À 
Comme T (z)- =z” E IT A ` d T -) 
2 


Si u-v=2p+1 т, (2) = Relr, e)» 1, 27 (EE Jr (s) oor p (LAA LE À) 
2 


weini T, (z)- веј, (з (z)- 277 al Eee enn ( ER ev LAY |; 3 
2. 
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Comportement asymptotique uniforme de la fonction modifiée de Lommel de première espèce 
lorque v -> °° 


Dans l'ouvrage de 1974 de F.W.J.Olver, W.Rheinbolt < Asymptotics and Special Functions > en pages 
388-390 formules 11.01, 11.02, 11.03, 11.04 et suivantes (prendre m=u+1 et z pour є=2/№), les 
auteurs donnent en exemple la construction d'un développement asymptotique uniforme de la 
fonction de Lommel modifiée. Après identification par tâtonnement il s'agit de bien de la fonction 
de première espèce. En l'occurence la construction est la suivante : 


ud ¿2 


Ge REES 


d? 
aci Get) а) pour 1>0 


, 4 4 G 
1+6? dÇ 


1=0 v 
gH Lu 4972 2 4 2 
Avec. Go(u v.c) - е 5 G(u,v,c)- £^! (+в) +2& 6 2 c*(1- u) 
1+6 T 
G, (u,v ozer enit +2u и?) act (us -34u -38u? бш? + 3u* | 4c (49 - 39 äu +5? "d Su ND 3} 
YT 7 
(+67) 


En posant Gul ç *1G (u,v,¿) 


u+l 


2 e d d? ~ 

= G(u,v,6)= çH G(uv,c)- == ( б, \(и,у,&)+(2и +3) ac Qe Ce God) 
~ ~ 2 ~ 

(is) Quee ас Gi (161.6) 6? а Gens) 


G(u,v,¢)= 1+? 


и) (L+ u c6 и?) Cf 

u le) 
E (Le n) ac 2n u?) 26 *(175+34u 38u? + 2 ch 39и+2и° +5и? ut) &*(u-1y (u-3) 
d 1+? 


Avec u ETT 


Par déduction, il est possible de trouver le comportement asymptotique de la fonction de Lommel 
de première espèce (voir plus haut). 
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Fonctions modifiées de Lommel de paramètre v=it purement imaginaire 


La fonction modifiée de première espèce est à valeur réelle lorsque l'argument est réelle : 


F p uU tir 3 и-їт+3 z 
14 2| ^ ? Â l=+0 21 m-l 
u+l = 


eler EE avec aen) Hm +e) 


(и + 1} +r? 1=0 dj (и,т) m-l 


l uir (z) SE 


Quant à la fonction modifiée de deuxième espèce elle n'est pas à valeur réelle lorsque l'argument 
est réelle : 


. E j TET m 
T, (2)=t, l) (аа "lu Zen oda?) 


2 2 


T 
i OEIT E ç 
AEN e e Conf rt) ref, (SK, (G) 
T T 


А Sin z ET соз п Lead 7 ET 
л 2 2 2 
= Belt, (z)}+ к.а) E чых E ES iSinh(zr) 
T 2 2 
к) Ç MT (соз x ET 
л 2 2 
De plus si: E =< 6 E =] St Ë = ode = = iSinh(zr) 


Et 2со{л EH bd = = Cos(zu)+ Cosh(rr) 


us t = a - Refr, Li. Н) к, (ejr: Costau)+Coshlae) p („ 


Lii (z)- 


T. .(z)= AC + E it e + - it fret, Gi. Sin(zu) „ (osi Соз(ти)+ Ок), 


T T 


La partie imaginaire de la fonction modifiée de Lommel de deuxième espèce est donc 
proportionnelle à la fonction de Bessel K d'ordre purement imaginaire. La partie réelle de la 
fonction modifiée de Lommel de deuxième espèce est plus intéressante : 


кеў ЕЕ) О ды су, EEk, e) 


2 2 
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Des transformations intégrales et leur inversion qui fonctionnent plus ou moins bien 
numériquement 
Le but ici est d'illustrer certaines généralisations de la transformation intégrale de Kontorovitch- 


Lebedev, notamment avec les fonctions de Bessel J d'ordre imaginaire, ou bien encore les fonctions 
de Whittaker W. 


Certaines formules de Titchmarsh, 1962, dans « Eigenfunction expansions associated with second- 
order differential equations. Part 1 », page 95, formules 4.14.1 et 4.14.2 (il suffit de changer 
g(x)=xf(x) et VA=r, и=ух, у=уу) se présentent comme suit : 


DI ле SC 4nJy, (yx С уу) 2 [ar с г; Ce JE mS = dee 
1 n=+00 
(2) f(x)=- > (4n + 2an (yx С ro} 3L SES Gre (yx J mU; (уу) rO) 


La transformation se présente donc avec un terme de spectre continu et un terme de spectre discret, 
notons les comme suit: f (x)= wx) + f (x) i 


terme discret terme continu 


La représentation de ces transformations intégrales sous cette forme présente un meilleure stabilité 
numérique. 


Par contre lorsque les formules précédentes s'écrivent sous la forme : 


@ LE anre) [rs dro) far Ren l| dv Ren. Gron) 
@ л) (m ser) jr dro) Z tn c mU fdv mo) 


Alors l'instabilité numérique dans le calcul est plus manifeste en général. 
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Plus généralement, S.B. Yakubovich 1997 < Index integral transformations of Titchmarsh type > et 
1998, < The Titchmarsh integral transformation by the index of a Bessel function >, écrit le couple 
de transformations (avec une valeur y=1) pour une valeur du paramètre u€[0,1] : 


O a forro Zë DE wl 1e S fae este Ela Air 
o S es ean (EE) M ea > WEE r Sinh(zr) ee], et 


] T 7 
es [ar t Sinh[zz) 


0 


i za Re(J,,(x)) |= at) Im(, (x) rer) 


Cosi Z) SR 
2 2. 


r 1-и+їт r 1+и+їт 
2 2 


т5їпї(лт KS Re(J, (x))- us Im(J, (x)) 
jida Snai 2 " Cosh| = Sinh 
f9- 27 [dr esee) u д | s, dete) је x EC glo) 
rSinh(zr . Re(J, (x))- iu) In(4, (x) 
rise Z ar comite (Iz s nate) L far + 5) 


2x ли TT ли ЛТ stc) 
0 Cos? Cosh? 4 Sin? Sinh? 
2 2 2 2 


по) Jar estes dir, gate) 27 ut) AS it) mU, (s) s) 


A x3 — Cos(zu)+ Cosh(zr) SR d sm Z) 


L'utilisation de la fonctions de Lommel de deuxiéme espéce introduit de fortes instabilités 
numériques que l'on a pas dans les formules originales de Titchmarsh. Exprimons la formule de 
Yakubovich sous une forme permettant de revenir aux formules de Titchmarsh, il vient : 


(=) 
2 


BE "la Se) [ол S up un 


H 


f()= 


+œ 
Í dt t Sinh(zz) 
0 


27 
2 
T xX 


2 +œ si Z) 
SËCH йу ДЕ R. 
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Pour éviter là encore des instabilités numériques j'ai remarqué que les expressions suivantes sont 
plus commodes : 


-u +œ 


f(x)= = far rSinh(zr ) 
T X 


d 2 
(=) š 
2 


e | Sin © Cos| É d Sin] 2А. Cos| H 
x j^ ; Li CH NS dÉ ш(/„(х)) jo fb») M va Im 


Mais attention ce n'est pas vraiment sous cette forme que j'ai testé sa validité car elle présente un 
aspect divergent sur le terme comportant la fonction de Lommel de premiére espéce. Pour 
« régulariser » ce terme, il sera transformé en un terme discret par l'utilisation du développement 
de la fonction de Lommel de première espèce et l'intervention de la transformée de Laplace de la 


fonction f(x). 


Pour u=0, on trouve l'expression (et son équivalent plus stable): 


f) ar ess) [de rm (42 fae zitt fär rom, O) 


e id [ar ess) fab Om o), fae | EE imb (|a ONE mon 


Pour u- 1, on trouve l'expression (et son équivalent plus stable): 


+œ 


1 Sinh(rr) 


л) аг: @ (2) sao aide 092 fae net Блокер, 0 
9 Cosh| =Â 0 0 0 
2 
= È [ar а D [dr Окер, |4: < aR (s) f dy ORe) 
idee far ss) dr Луке, far deet et de ro or ве) 


En développant la fonction de Lommel de première espèce en série dans les premiers termes des 
membres droit des expressions, on retrouve les formules de Titchmarsh : 


п=+со 


2 + 
(1) f(x x Ку 4nJ,,( yx СЯ y)/(y)+ me | 


n-l 


TT 
Sinh(zc 


Re, (их) [а Ref, (УО) 


n=+00 +оо 


(2) л) У (n 2), )[ Фуу »)/б)+ 2 [dr SE yim Gi f Im[, (yy) 


X Жї Zx = Sinh(art 
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En effet utilisons deux vésultats, sous certaines conditions de convergence d'intégrales (voir les 
articles de S. Yakubovich) : 


Rd -izLog( p+ p=) 
Transformée de Laplace de J, (x) > [dy e” J. (y)= É = 
0 


p +1 


KE Cosle Log(p+ Ji +1) 

[dy e” вел, (v))= - 

0 Jp +1 

eoru) Lena) _ вн (гов (р + fr) 
2i p^ +1 ур? +1 

Identités restreintes aux fonctions Rel J * (y), et ШАК Е (y), et J, (y) 


кде Í dp Col Lerle + Fu), 


[dy e” im, (}= 


p' +1 


або) = [v fre, (y) 


Si 


са 


ji THC Som, (0) ou (о) = jo =ош 
QU d zj JON, ou glv “| Kp) 
g(v) K DIAS g(v) [dp GH 


ой h( p)est "l'inverse" de la transformée de Laplace de f > f (x) - [ар e пр) 
0 


Les deux développements еп série des fonctions de Lommel de première espèce sont les suivants : 


1=+00 2141 l1=+œ 1? 
sk (x)=2 У: g) Jara (x) 9 (x) =8 Y. == Jais) 


eL +r — AP +? 
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Injectons ces deux développements dans les deux expressions obtenues, et l'on retrouve bien les 
formules de Titchmarsh : 


Intégrale | dr 


—al +œ М 
ro ete gue 
0 


()= п A E fae e TU — [ar aa ms fo stet) 
()= f()- A Seny T la т Ze [ат = Jp ` UE sm mU (e) ® ro mb. G9 


+оо Si i 
k -te ah) q = - |р? +1 [pas d 
0 


яси +1 jJ. Eos)" 
+ 


Í 


P (p4 TES у"? — mU. (х)}[ ay fom, G2) 


* Sable) 


()= A) = 330r јар À 


(21+1 
: +œ EM p È Jp +1 ES = 
Transformée de Laplace de J. (x) => [ау grs | = 
0 ур +1 


OS il 2 Barat, аа) [dp pen far mb Co f 7 шалу) 


X 120 


+œ 


OS s) E Q9 ауа)» far pur 


ls ri)» X fae co x Jal) fa iescht) +2 far EICH 


ex лт) 


OSs- Yrs) [aeos [dp cobi y, e EEGEN 


See gel. [s ente or) +2 far = E Rel E) dv лоҝе, (v) 


AP +r? Sinh(xr) 


far Соѕ(ат) me^" Nx | » [ат CosfeLoglp [p +1) zo | a) 
! A 


4P +r? 


— KP) (p. J)” +2 far Re) | dv FORO) 


= o үр +1 Sinh(z) 
4e Э 21 
Transformée ае Laplace de J, (x 19 e T J, (у )= P+ ы +1 
p +1 
2 +2 T +00 
EE KE уат je e P” J (y ) + 5 as Six X NG f (y)RetJ;.(v)j 


(2) = sa) Fuul Vo ro rope far roe (3) [dv лодар, О) 


Sinh(zz 
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Avec un paramètre u €[0,1] quelconque il faudrait exprimer le terme discret w(x) : 


ок JL cSinh(zc) (tt). efe ШЕ; SE E 
S ini ZE 
E 2 je (Et. eo 3 КТЕ ЗЕ E)r Ее) 


sous Іа forme d'une somme de fonction de Bessel d'ordre réel іпјесіопѕ le développement de Іа 
fonction de Lommel de première espèce en fonctions de Bessel et appliquons une transformation de 
Laplace sur les parties réelles et imaginaires des fonctions de Bessel : 


>м] Jes S U+ a+ Ur + +1) 
шїт 


J > 
Cu IGI и +1) +r?) ara (x) 


HF] RA fo Se вер, 


d (21+u+1) +r s 


l-u+it i 

2 (A+ u+ + ul) n E^ 
w(x)- > n Жан) [de т$їпї(лт) g Т; S P 
T x To ! 0 tuti +T E 0 


l=+% +00 
wx) = — ША SD т$їпї(лт) 


HX. o I! 


ТЕ M dk = 
2 | 2 


[ay fom) 


=i 


Р sim Z - JE Z ансор pr) 
ies pr E Ae + Co = eof = ово À fr) 


+œ 


{ 


„йд: (A+ u Ee, ts 


h 
TXT H (p) 


: CoslrLoglp + |р? +1 L 


d Qreuslf «e 


eur] 
2 


С az) 
muye 2 
+Cos dt t Sinlr Lo | +4 ^4) 
ej (21+u+1) +c? TERME 


4 «(201+ ur ul 
le ЎА) (of 


T X ш 


2 


oyp +l 


On note ainsi les intégrales : 


2: 
Co. 
2 


(21+u+1f «c 


Sin Loglp+ p +i) 


4 (U+ u+ (+ u+] 
= wx)= ( H |Г@+и+1) 


re I Jaa N, 


Et l'on constate que numériquement le terme w(x) comporte < disons moins d'instabilité 
numérique. 
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Pour les valeurs u=0 et u=1, cela donne : 


n ns с 


(1+1) +r? 


«Je S dp Mp) dr 27 Cos[eLoglp fr) 


21+2Ÿ «c? 
(21 «2) 


Е EN 7 жа ар xe 2л сов[гТовр+.[рї +1] 
np) Хт ик, [ае AP +r? 


дие l'on doit comparer ауес les développements suivants : 


_ 24% _ 4 8S T dp "d лт Side Loglp Je? +1) 
u=0= w(x)= ES ey Sech da о) Лаб) Zell ach l ues (1+1) +r? 


итә wo) =! Euys) [ay r0) EP H [a клы 


x T o Vp? +1 б AP +r? 


Soit l'évaluation numérique doit conduire aux identités terme à terme suivantes pour les deux cas 
u=0 et u=1 : 


i +œ 27 dp W: кисы ы) 
и=0= kaa ол" [а (21+1) + 


Í dp ip) [ac TOME +1) 


р? +1 d AP +r? 


Je résume ainsi une implémentation possible (peut-être plus stable numériquement) de la formule 
généralisée de Titchmarsh lorsque l'on connait < l'inverse > de la transformée de Laplace de la 


fonction f(x) : 
MES rot 


E rhel «pu Е Jes än ) [drr Co Z EE pa) 


TE "Daf EE eue) e 
oü р) езі "l'inverse" de la transformée de Laplace de f > f(x)- [Ф ер) 
0 
(A+ ERU RENTES! 
ё zy ` ; у. Дд 
mx 1-0 

T | ek cof z) къ ni 3 cof z 

таеш 2 ар) Zäit ill M 
ЛТ TT 


rT x ï 2[Cos[mu)+Coshízr) Co z) 
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Instabilité numérique de la transformation de Yakubovich-Titchmarsh 


La formule de Yakubovitch-Tichmarsh présente un gros inconvénient d'instabilité numérique dans 
on implémentation en effet en la forme : 


Lo [ T ш КУЙ. e Р 
0 [e z | Sinh| = ) | 


[1-и+їт 
dese 


Tandis que sous sa forme développée en deux termes, l'instabilité numérique est réglée moyennant 
le développement en série de la fonction de Lommel et un passage aux transformées de Laplace 
dans les intégrandes (voir ci-dessus) : 


e 


о. 2 far rSinh(zr) S,.(x)g(r) 
zx n 


2 ad" Cos z) 
_ 2" Pas Lait, (үү ы. > NAZ 
f) far Sine) (ZH) s, (у S ST d ini b) 


E Í REEL) za Re „ы ы 1. (3) j A fb ТЕ a. eil GI 


ell [e A _ Im 
mx 2(Coslzu)+Cosh(zr) SR иң ZE) ô Cosi Z) ad? 


Pour se convaincre de cette instabilité potentielle reprenons la première ou deuxième expression, et 
voyons avec le développement asymptotique au premier terme ce qui se passe avec les fonctions de 
Lommel soit de première, soit de deuxième espèce : 


xtd P 
S uv (x) d E = È Suir (x) © rax 
1 е: ad? Cos Z) 
Spa) „(+2 | КЕ!) Ref, (x)}- im, (x)} 
| i 2 лт 
Cost Z | 


Posons  tE(x,r)=Vt? x? +t Log] ————= 
бт) à т+уүт +x? 


лт 1 лт 


Rob, Ge сао) elt l E ale? 
sin Z) "S co 


D ——7— Rei; — mU, (x)! =. - ZÉ ТЕ, 3 
SE SU Gg 2 " 


2 


лт 


zl "le 274 


= nte 31 dn +H+ 3 
2 
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ce qui donne une estimation des intégrales : 


f [x] 22 far t Sinh[x с] | [5255] p [2 Su cExI [ f[y] STET aa tyll UST PUR tyll 

== nit vto ETES HDN- сүл] ze Dr 

e 2 2 

f[x] = w[x] + f[x] 
w[x] = 2% far t sannta) D Eed p sy. i c [x] Í f[yl Sy [77] Re[2; Di? кы кш sss, [у11 

= ЕТЕ ) [> Su,it У f[y се =n [=] ir[y Sinh[ =] ir[y 
боа e fa LL емы l Tupaia S| Z] Lao, туп! Jena ЗАА]. ара, КЕ Бык ett 

x Cosh[z и] + Cosh[z т] Cosh [22] Sin nh [Z2] Co: sh [T] Sin nh [ = ] 
Siu d E e Louie UR „езу үл IE ax | ар н S Le ашы 
A x Cosh[z и] + Cosh[z т] Cosh[Z*] SES sinh [22] SES 4 dii Cosh [22] TEM Sinh[ ==] ы 


э f[x, N] „ 2 far san ZZ et a+ Z] fa et lsi [> - © ety =] + 2] 
жге" еш 2 ? 4 T 2 ? 4 
N LJ 


sat Ë 1- i ^ sin[ zz] Cos [Z] 
w[x, N] = = fe саһ] |r [=] |? sui c [x] [etn —  Re[3: о[у] - => 1т[34+[у11 
n? х 2 T Cosh [ =] Sinh[ Z] 
N e 2 2 
[x, N] 42 xH Fa s Ui f: fly] Si [= el 1 2 || 
= w[x = t —— y fiy] Sin] — - t Ely, t] + = 
S хх т? + x? 2 ~ 4 
e 


On voit que les intégrandes du deuxiéme terme sont des fonctions bornées, et qu'il y a de bonne 
chance sans le démontrer formellement que les intégrales convergent rapidement. Mais pour le 
premier terme, l'intégrande en т est une fonction exponentielle croissante de ce dernier rendant 
quasiment impossible un calcul numérique sans un processus de régularisation adéquat. C'est 
d'ailleurs celui que l'on réalise lorsque l'on substitue à la fonction de Lommel de première espèce 
sa représentation en série et que l'on change d'espace d'intégration avec les transformées de 
Laplace. 


Avec l'expression modifiée du premier terme en prenant les variables d'intégration de Laplace, nous 
avons : 


| ‘su z) | 
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x x Eo ! i ° 4p | {= кү| 


pon е 
(21+ u +1) +r? 
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(z 


| 
ou (ros p+ o? +1) 


Еп restreignant l'intégration sur des valeurs suffisamment grande de т de telle manière à appliquer 
les formules asymptotiques, il vient : 


SEND . A2 
Г ЕЕЕ Sinh| ZE |x 2"ac" Г и 
2 2 2 


2 !22(2/+1+ u) * xp) F? 
N)x DIU +14 d, d 
=> Р(х, ) == > n ( Vus] р x | T 


prse Z | eloelo Hp ) 


(21+1+u) +c? 


Expression qui а l'infime avantage de présenter des intégrandes bornées, et dont Іа dernière 
intégrales est elle méme une fonction de t probablement bornée lorsque le paramètre u€[0, 1]. 
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La transformation d'Al-Musallam, 2002, < A Modified and a finite Index Weber Transform > : 


L'identité de Hardy, 1926, R.G.Cooke « The Inversion Formulae of Hardy and Titchmarsh » 


+œ 
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Une identité de Hardy équivalente : 
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Passage de l'identité de Hardy à des transformations intégrales connues 


La transformée de Hankel classique : 


‚булугу (X) =2""T(v)J,(x) 7(y)= os (ху) 


Ly +2 


u=v-l+emu+v=2v-l+e et ` E e 
Lim —— = 2 Tend? 5l dy xy Lim? v- н O) 3€, 


SS SÉ 


Ge 
= Jade, (x)/ (y) 


La transformée de Neumann-Struve : 


Lued 


TEE as Ge) À 
= fe 70, (x) 
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Une transformation due à Jet Wimp. 1964 < A class of integral transforms », formules 4.10 et 
4.11 page 37 : 


g(r)= dÉ -u ir) 


f(x) = ER e т$їпй(2лт)й/ „ү (x)e(c) f(x)- ER | ат s Sin(2ze)t( 4 - и ir) 


2 +00 +00 


Гат, Go ` [at Јат, „Сл 


Helle) 


Lorsque u=0, on retrouve la transformation de Kontorovitch-Lebedev : 
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Une transformation de Noyau une fonction hypergéométrique diversement décrite dans plusieurs 


articles 
Les articles les plus précis et « concordant » quant à la forme sont les suivants : 


Dans l'article de 2001 de Y.A.Neretin « Index hypergeometric transform and imitation of analysis of 
Berezin kernels » , ainsi que dans son article de 2012 « Index hypergeometric integral transform », 
elle est décrite avec divers facteur multiplicatif qui une fois corrigé et testé numériquement donne 
ceci : 


g(r)= | [ax fem ear ete) 


T(b c) 5 b+c 
KEE NEE 
On peut donner une forme similaire : 
Comme |r(2it)f = AUGE 29 F s G 
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Dans l'article de 2006 de S.Yakubovich « ON THE PLANCHEREL THEOREM FOR THE 
OLEVSKII TRANSFORM » , les deux expressions 1.1 et 2.2 utilisées sont les suivantes : 
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Ces deux systèmes de transformation sont bien équivalents. En effet : 
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E.C.Titchmarsh en pages 101-102 de son ouvrage < Eigenfunction expansions associated with 
second-order differential equations. Part 1» datant de 1949 dans sa première version, donne la 
transformation suivante : 
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Avec les transformations suivantes ce systëme est bien équivalent à celui donné par Neretin : 


r= JA э T =2dr Bed y=b+c 
J 
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Voici également deux publications faisant état d'une transformation intégrale dítes de Fourier- 


Jacobi, ou Jacobi tout court : 1974, T.Koornwinder, < A new proof of a Paley-Wiener type theorem 
1977, Y.Hasegawa, « On the integrability of Fourier-Jacobi 


for the Jacobi transform », 
transforms » présentant la transformation suivante et son inversion : 
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Là encore c'est la même transformation décrite par Olevskiy, puis Neretin et d'autres chercheurs 
cités dans ce document : 
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Avec les substitutions suivantes, le système est bien équivalent а celui de Neretin : 
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Deux autres formes de cette transformation sont encore présentes dans le livre des auteurs 
S.B.Yakubovich et Yu.FLuchko « The hypergeometric approach to integral transforms and 
convolutions », formule 7.1 page 121, ou il est écrit : 


eO- ja re ln 


C 
c-a-1 


x Nn | | А | с-и+їт,с-и—їт 
f(x) = zr | T 5їпһ(2лт )е(т)Г(и +їт)Г(и —їт)Г(с -u +їт)Г(с— u H = =x) 


Résolution des problëmes aux limites de Laplace en cylindriques 3D par les transformations de Hankel et les 
transformations de Kontorovich-Lebedev- p 212 


Et dans le livre de S.B.Yakubovich < Index Transforms >, pages 200 et suivantes, les expressions 
suivantes s'avérent étre équivalentes : 
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Ces deux systèmes sont également équivalents à celui de Neretin, en utilisant les propriétés des 
transformations linéaires des fonctions hypergéométriques, à savoir : 
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Une formule due à Jet Wimp, 1964 « A class of integral transforms >, formules 4.13 et 4.14 page 
38: 
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Encore une fois се systëme est équivalent à celui donné par Neretin , comme il suit : 
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Il reste à effectuer quelques derniéres substitutions pour retrouver le systéme de transformations de 
Neretin : 
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On voit donc que beaucoup d'articles présentent des transformations en apparence différentes mais 

qui en réalité sont les mémes. On peut également conclure que le systëme donné par Neretin 

apparait donc comme étant dans sa forme la plus simple. Par ailleurs notons que la fonction 

| Se ире Si est apparentée d ипе fonction de Jacobi (définition hypergéométrigue et поп 
b+c 

elliptique), la raison pour laquelle la transformation est également appelée « Transformation 

intégrale de Jacobi ou de Fourier-Jacobi » par certains auteurs. 
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Forme simplifiée de la Transformation d'Olevskyi b=c 


Avec Б=с, il vient : 
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